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Abstract 



A halting protocol is generally irreversible in the classical computation, 
but surprisingly it is usually also irreversible and non-unitary in the uni- 
versal quantum computational models. The inherent incompatibility within 
the universal quantum computational models between the irreversible and 
non-unitary halting protocol and the unitary quantum computational pro- 
cess that obeys the Schrodinger equation in quantum physics has been known 
since the early set-up of the universal quantum computational models. The 
irreversibility and non-unitarity of the halting protocol is closely related to 
the inherent irreversibility and non-unitarity of quantum measurement in 
quantum mechanics. The unitary dynamics quantum mechanically implies 
that the halting protocol in the universal quantum computational models 
should be made reversible and unitary so as to eliminate the inherent incom- 



patibility. It has been shown in Ref. [24] (Arxiv: quant-ph/0507236) that 
a universal quantum computer could be powerful enough to solve efficiently 
the quantum search problem in the cyclic-group state subspaces, and the re- 
versible and unitary halting protocol is the key component to construct the 
efficient quantum search processes based on the unitary quantum dynam- 
ics, while the state-locking pulse field is the key component to generate the 
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reversible and unitary halting protocol. In this paper the reversible and uni- 
tary halting protocol and the generalized state-locking pulse field have been 
extensively investigated theoretically. The basic principles to construct the 
state-locking pulse field and design the reversible and unitary halting pro- 
tocol are described and studied in detail. A generalized state-locking pulse 
field is generally dependent upon the time and space variables. It could be 
a sequence of time- and space-dependent electromagnetic pulse fields and 
could also contain the time- and space-dependent potential fields. Thus, the 
reversible and unitary halting protocol built up out of the state-locking pulse 
field generally consists of a sequence of time- and space-dependent unitary 
evolution processes. It is shown how the quantum control process is con- 
structed to simulate efficiently the reversible and unitary halting protocol. 
An improved subspace-reduction quantum program and circuit based on the 
reversible and unitary halting protocol, which is much simpler than that one 
in the previous paper [24], is proposed as the key component to construct 
further an efficient quantum search process. A simple atomic physical system 
which is an atomic ion or a neutral atom in the double-well potential field is 
proposed to show how the state-locking pulse field is generated and how to 
implement the reversible and unitary halting protocol. 

1. Introduction 

A halting protocol (or a halting operation briefly) of a computational 
model such as the Turing machine is one of the key components in computa- 
tion. It is well known that the halting operation generally is not a reversible 
operation in the classical computation and often is related to measurement 
of computational results. Though in the reversible computational model [1] 
and probably in the quantum Turing machine [2] the halting protocol could 
be made reversible, it is surprised that the halting protocol usually is not 
yet reversible and unitary in the universal quantum computational mod- 
els including the universal quantum Turing machine [3, 9] and the universal 
quantum circuit model [4]. When the universal quantum Turing machine was 
proposed in the early day [3] , the halting protocol was also introduced as one 
of the important components of the universal quantum computational model. 
According as the universal quantum Turing machine, in addition to the com- 
putational quantum system that is used to perform quantum computation 
there is also an extra quantum bit used to instruct what time a quantum 
computational process is halted on the universal quantum Turing machine. 
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This extra quantum bit is named the halting quantum bit or briefly the halt- 
ing qubit. The halting quantum bit should be observed periodically from the 
outside in a non-perturbation manner so that the quantum measurement on 
the halting quantum bit does not disturb the unitary evolution process of 
the quantum system during the quantum computational process. Once the 
halting qubit is found to be in the halting state, e.g., the state |1), the quan- 
tum computational process terminates. Therefore, the halting operation is 
performed to stop the quantum computational process only after a sequence 
of non-perturbation quantum measurements on the halting qubit were car- 
ried out during the quantum computational process. However, the halting 
protocol could not be compatible with unitary quantum computational pro- 
cesses within the universal quantum computational models. It is well known 
that a quantum computational process on a universal quantum computer 
obeys the quantum physical laws [3, 4], that is, the time evolution process 
of the quantum system of the quantum computer obeys the Schrodinger 
equation in quantum mechanics during the quantum computational process. 
It is also known that the quantum parallel principle [3] that the quantum 
computational process is performed on a superposition of the basis states 
of the quantum system is the characteristic feature of the universal quan- 
tum computation models. The halting protocol is essentially different from 
a unitary quantum computational process in the universal quantum com- 
putational models in that the halting protocol involved in the non-unitary 
quantum measurement is generally irreversible and non-unitary. This basic 
and inherent incompatibility within the universal quantum computational 
models between the facts that the quantum computational process obeys the 
unitary quantum dynamics and that the halting protocol is irreversible and 
non-unitary due to the non-unitary quantum measurement really originates 
from the basic quantum mechanical laws. It is well known that there is an 
inherent incompatibility in quantum mechanics [5] between the non-unitary 
quantum measurement and the unitary time evolution process of a closed 
quantum system that obeys the Schrodinger equation. 

For many years since the early proposals of the universal quantum com- 
putational models [3, 4] the computational power of the universal quantum 
computational models has been investigated extensively and continuously. In 
the past two decades it has been shown that there are several possible candi- 
dates to fuel the quantum computational power, which include the superposi- 
tion and the quantum parallel principle [3, 4], the quantum entanglement [29] 
and the multiple-quantum coherence [17, 20], the quantum coherence inter- 
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ference [18, 20], and the unitary dynamics quantum mechanically associated 
with the symmetry and structure of a quantum system [17, 18, 20, 22, 23, 24]. 
A large number of quantum algorithms [13, 14] including the Shor's prime 
factorization and discrete logarithm algorithms [6] and the Grover's quan- 
tum search algorithm [7] which outperform their classical counterparts have 
been discovered and developed in the past two decades. Most of these quan- 
tum algorithms are based on the universal quantum circuit model [4]. The 
computational results of these quantum algorithms generally are output at 
the end of the computational processes by the proper quantum measurement. 
Thus, these quantum algorithms are not related to the halting protocol of the 
universal quantum computational models. It has been believed extensively 
that power of the quantum computation that outperforms the classical com- 
putation could come mainly from the quantum parallel principle [3]. But it 
was suspected whether or not an arbitrary recursive function mathematically 
could be computed more efficiently on a universal quantum computer than 
a classical computer with the help of the quantum parallel principle [8] . The 
quantum parallel principle requires that a quantum computational process 
take a superposition state as its initial state. However, an incompatibility 
within the universal quantum computational models arises when a universal 
quantum computer computes a general recursive function in mathematics by 
starting at a superposition state. This incompatibility is really due to the 
conflict between the halting protocol and the quantum parallel principle [8]. 
It results in the question whether or not a universal quantum computer is 
capable of computing more efficiently an arbitrary recursive function when 
the quantum parallel principle is employed. A quantum computational pro- 
cess allows its initial state to be a superposition state of the quantum system 
of a universal quantum computer [3, 4, 9]. A superposition state may be 
expressed as a linear combination of the conventional computational bases 
of the quantum system. Each computational base of the superposition state 
could be thought of as one input state of a quantum algorithm running on the 
universal quantum computer if the superposition state is taken as the initial 
state of the quantum computational process of the quantum algorithm. This 
implies that the quantum algorithm can be performed in a parallel form on 
the universal quantum computer by taking at the same time all these differ- 
ent computational bases of the superposition state as its input states. From 
the point of view of the classical computational model a different input state 
corresponds to a different classical computational process. Then one could 
imagine that in effect the quantum computational process really performs 
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simultaneously many different 'classical' computational processes with differ- 
ent input states when the initial state of the quantum computational process 
is a superposition state. Although these different 'classical' computational 
processes start at the same initial time, they could end at different times in 
computation, respectively. As a typical example, this situation will occur if 
the quantum algorithm is used to compute a general recursive function in 
mathematics and the end state of the computation is fixed. Because there 
is not the same end time for all these different 'classical' computational pro- 
cesses a conflict arises when the quantum computer decides what time the 
quantum computational process is halted. Thus, from the viewpoint of the 
mathematical logic principle of computational model the halting protocol is 
incompatible with the quantum parallel principle. 

This conflict is apparently between the halting protocol and the quan- 
tum parallel principle, but it really related inherently to the non-unitarity of 
quantum measurement which is well known in quantum mechanics [5]. Ac- 
cording as the halting protocol the halting quantum bit usually is set to the 
halting state \nh) = |0) at the initial time to of the quantum computational 
process. Suppose that the initial superposition state of the whole quantum 
system including the halting qubit of the quantum computer is given by 

!*(*„)) = ^(^K)!^)} = X>(*o)|o>|</^o)>- 

According as the halting protocol, when some 'classical' computational pro- 
cesses arrive in their end states during the quantum computational process 
their halting state \nu) becomes the state |1). Since there is not the same end 
time for all these different 'classical' computational processes, there exists 
some time t in the quantum computational process such that some 'classical' 
computational processes arrive in their end states ({|<ft (£))}) and hence their 
halting state becomes the state |1), while the rest have not yet arrived in 
their end states and their halting state remains unchanged and is still the 
state |0). Then at the time t the quantum state of the whole quantum system 
is written as 

!*(*)> = £>*(*)|o>b fc (f)> + 5>(*)|i>bi(*)>- 

k l+k 

The state \^{t)) is clearly a superposition state involved in the halting state 
\rih) and the computational states {\<Pj(t))} of the quantum system. Again 
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according as the halting protocol [3], the quantum measurement is carried out 
on the halting quantum bit periodically from the outside during the quantum 
computational process. Then the quantum measurement will change the 
state \^(t)) and spoil the quantum computational process [8]. Actually, 
according to the quantum mechanics [5] the quantum-state collapse on the 
superposition state \^(t)) occurs inevitably when the quantum measurement 
on the halting qubit is performed on the superposition state \^/(t)). It is 
well known that a quantum-state collapse process is a non-unitary process in 
quantum mechanics [5b]. Thus, the halting protocol generally is irreversible 
and non-unitary in the universal quantum computational models in which 
the quantum parallel principle is a basic principle. 

There are a number of works [3, 8, 9, 10, 11, 12] to discuss in detail 
the conflict and to propose schemes to avoid it for the universal quantum 
computational models. However, so far there is not any satisfactory and 
universal scheme to avoid the conflict when the initial state for a quantum 
computational process is a superposition state. On the other hand, the halt- 
ing protocol may be made reversible in the reversible computational model 
[54] due to that any initial state of computation is generally a single basis 
state and hence there is not such a conflict in the reversible computational 
model. It has also been shown [3, 8, 9, 10, 11, 12] that this conflict could 
be avoided when the quantum computational process consists of a single 
'classical' computational process or several different 'classical' computational 
processes which arrive in their end states at the same time. This is because 
in this case the state \^{t)) above is either the first term with the halting 
state |0) or the second term with the halting state |1) but not a superposi- 
tion of the two terms. Since an initial basis state corresponds one-to-one to a 
'classical' computational process in the quantum computational process one 
may also say that the conflict could be avoided if the initial state of the quan- 
tum computational process is limited to a single basis state. Although from 
the viewpoint of the mathematical logic principle of computational model 
the halting protocol becomes reasonable if the initial state of the quantum 
computational process is limited to a single basis state, it could not be gen- 
erally reversible and unitary as the halting protocol contains the non-unitary 
quantum measurement. The original halting protocol [3] uses the quantum 
measurement to achieve the halting operation, that is, if the halting qubit 
is found in the halting state |1) through the periodic quantum measurement 
the quantum computational process is stopped by the brute-force method 
from the outside. Therefore, even in the case that the initial state of the 
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quantum computational process is restricted to be a basis state the original 
halting protocol generally is not reversible and unitary. The original halting 
protocol [3] also implies that there must be a conditional operation such that 
when the computational end state appears during the computational process 
this conditional operation converts immediately the initial halting state |0) 
into the state |1). Some improved halting protocols [11, 12] use explicitly the 
conditional halting operation to replace the original halting operation. It was 
also proposed in the halting protocol [11] that when the quantum computa- 
tional process is stopped due to that the halting state |0) is changed to the 
state |1) it needs to start an extra unitary evolution process at the same time 
so that while both the computational end state and the halting state |1) are 
kept unchanged after the halting operation, the total evolution process can 
be unitary for the whole quantum system including the computational quan- 
tum system, the halting qubit, and the auxiliary qubits (the halting protocol 
[11] needs to use an extra auxiliary register). In the reversible computa- 
tional model the reversible halting protocol could be achieved by executing 
successively a computational process, the conditional (logic) halting oper- 
ation, and the inverse computational process [54]. However, the quantum 
measurement must be given up in these halting protocols [3, 9, 10, 11, 12] if 
one wants to make these halting protocols reversible and unitary thoroughly 
as the quantum measurement generally could lead to information loss of a 
quantum system even when the initial state is a single basis state. There- 
fore, the irreversibility and non-unitarity of these halting protocols in the 
universal quantum computational models are traced ultimately back to the 
irreversibility and non-unitarity of quantum measurement in quantum me- 
chanics. 

The halting protocol usually is not reversible and unitary in the universal 
quantum computational models due to the non-unitary quantum measure- 
ment, but for most of the present quantum algorithms based on these uni- 
versal quantum computational models the computational results usually are 
not significantly affected by the non-unitarity of quantum measurement even 
though the non-unitarity could lead to loss of partial information of the com- 
putational results. There are also the incompatibilities within the universal 
quantum computational models between the quantum parallel principle and 
the halting protocol from viewpoint of the mathematical logic principle and 
between the unitary quantum computational process and the irreversible and 
non-unitary halting protocol from the viewpoint of the quantum physics, but 
up to now a number of powerful quantum algorithms [13, 14] that outperform 
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their classical counterparts have been found within the universal quantum 
computational models and the halting protocol has not yet any significant 
effect on these powerful quantum algorithms. Therefore, on one hand, the 
quantum parallel principle is paid much attention, on the other hand, the 
halting protocol has a negligible influence on the quantum computational sci- 
ence in the past years. It has been explored in theory to eliminate the incom- 
patibilities within the universal quantum computational models, but people 
are not clear whether or not the incompatibilities could lead to an essential 
impact on the power of quantum computation and do not yet know whether 
or not the reversibility and unitarity of the halting protocol is important to 
discover and develop new and efficient quantum algorithms. However, these 
powerful quantum algorithms [13, 14] can only treat successfully few special 
mathematical problems. Most of them are based on the unitary quantum 
circuit model [4] in which the irreversible and non-unitary halting proto- 
col neither is used nor has any essential effect on the computational power 
of these quantum algorithms. It is also known that the incompatibilities 
could lead to that a universal quantum computer could not be more efficient 
than a classical computer in computing an arbitrary recursive function in 
mathematics. On the other hand, it is well known that a quantum computa- 
tional process obeys the unitary quantum dynamics and is compatible with 
the mathematical logic principles used by the computational process. These 
facts show that it is worthwhile to investigate further how the halting pro- 
tocol can be made reversible and unitary so that the incompatibilities could 
be eliminated in the universal quantum computation models. 

The unitary dynamical principle of quantum mechanics plays a key role 
in the scalable quantum computation in a mixed-state quantum ensemble. 
This basic principle states simply that both a closed quantum system and 
its quantum ensemble obey the same unitary dynamics quantum mechani- 
cally [17, 18, 20, 22, 23, 24]. The unitary quantum dynamics associated with 
the symmetry and structure of a quantum system has been used to discover 
and develop new quantum search processes in a quantum ensemble [17] and a 
pure-state quantum system [23, 24]. Quantum search processes are extremely 
important in quantum computation as the unsorted quantum search process 
has an extensive application in computational science and can be used to 
solve the NP-hard problems. It has been shown that the square speedup for 
the standard quantum search algorithm [7] is optimal in a pure-state quan- 
tum system [15]. The first attempt to break through the quantum-searching 
square speedup limit was carried out in a nuclear magnetic resonance (NMR) 
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spin ensemble [16, 17a]. Though these quantum search processes [16, 17] are 
not scalable in a spin ensemble, their speed is really exceed greatly the al- 
lowed value of any quantum search algorithm with the square speedup limit 
in a range of a few quantum bits. However, these quantum search processes 
do not achieve a real breakthrough of the square speedup limit because their 
output NMR signal intensities decay exponentially as the qubit number of the 
spin system increases and hence beyond a few quantum bits their quantum- 
searching speed falls off rapidly. So far a scalable quantum search algorithm 
working in a quantum ensemble has not yet been found. Both the unitary dy- 
namics quantum mechanically and the quantum coherence interference play 
a key role in achieving the fast and scalable quantum computation in a spin 
ensemble [17, 18, 20]. It has been shown that many oracle-based quantum 
algorithms including the parity-determination algorithm are subjected to the 
polynomial speedup bounds in a pure-state quantum system [19]. Again the 
unitary quantum dynamics has been shown to play an important role in 
achieving a much fast computational speed to solve the parity-determination 
problem in a spin ensemble which is beyond the polynomial speedup bounds 
on the oracle-based quantum algorithms in a range of a few quantum bits 
[18]. Up to now, the scalable problem has not yet been solved for the quan- 
tum parity-determination algorithm in a quantum ensemble. These results 
obtained in a spin ensemble encourage one to explore further the potential 
ways to break through the quantum-searching square speedup limit and the 
polynomial speedup bounds upon the oracle-based quantum algorithms. The 
unitary quantum dynamics associated with the symmetry and structure of 
a spin system has also showed that the prime factorization for a large com- 
posite integer may be implemented in a scalable form in a spin ensemble 
[20]. The efficient factoring algorithm was first discovered in a pure-state 
quantum system [6]. In order to achieve the scalable quantum computation 
for the prime factorization in a spin ensemble it is necessary to exploit the 
symmetric property and structure of the spin system of the spin ensemble to 
help the unitary quantum dynamics to solve the prime factorization prob- 
lem. Here both the time-reversal symmetry and the rotation symmetry in 
spin space [5a] of the spin system play a key role in the scalable factoring 
algorithm in a spin ensemble [20]. The factoring algorithm in the spin en- 
semble usually is divided into four time periods [20] , which is similar to the 
conventional NMR experimental counterparts [21]. The first period is to gen- 
erate the multiple-quantum coherences of the spin ensemble. In this period 
information of the order of the modular exponential function is loaded on the 
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multiple-quantum coherences. The second is the time evolution process to 
carry out the frequency labeling for these different order multiple-quantum 
coherences generated in the first period. The third is the time-reversal pro- 
cess of the first period. The final is the quantum measurement to output 
the NMR multiple-quantum coherence signal that carries the information 
of the order of the modular exponential function. Here the output NMR 
multiple-quantum spectrum usually is used to determine the order of the 
modular exponential function. The multiple-quantum spectrum may be ob- 
tained by Fourier transforming the output NMR multiple-quantum coherence 
signal that is measured indirectly in experment. The time-reversal symme- 
try [5a] ensures that the dephased NMR multiple-quantum coherences in 
the first period can be refocused in the third period and hence the multiple- 
quantum coherences become inphase so that the multiple-quantum coherence 
interference can lead to coherent enhancement of the output NMR multiple- 
quantum coherence signal in the factoring algorithm. Therefore, it becomes 
possible that the output NMR multiple-quantum coherence signal does not 
decrease exponentially as the qubit number of the spin system increases. The 
time-reversal symmetry has been used extensively to obtain highly sensitive 
and inphase NMR multiple-quantum spectra in high-resolution nuclear mag- 
netic resonance experiments in solid [21]. On the other hand, the rotation 
symmetry in spin space [5a] of a spin system is the basis of the multiple- 
quantum operator algebra spaces of the Liouville operator space of the spin 
ensemble [22]. Though the total number of independent NMR multiple- 
quantum transitions in an n— qubit spin system increases exponentially as 
the qubit number of the spin system, a large number of the independent 
multiple-quantum transitions are really degenerative or nearly degenerative 
in transition frequency due to the rotation symmetry in spin space, and the 
multiple-quantum transitions with significantly different transition frequen- 
cies are really very few in the spin system. For example, there may be only 
(n+ 1) different order multiple-quantum transitions in the n— qubit spin sys- 
tem and each has its own transition frequency, while the total number of the 
independent multiple-quantum transitions of the spin system is (4 n — 2 n )/2. 
In theory the number of the multiple-quantum-transition spectral lines in 
the NMR multiple-quantum spectrum may be equal to the number of the 
independent multiple-quantum transitions of the spin system. Note that the 
total spectral intensity of the NMR multiple-quantum-transition spectrum 
of the spin system is not more than the total magnetization (~ n2 n ) of the 
spin system in the thermal equilibrium state [55], while the total number 
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(4 n — 2 n ) /2 of the NMR multiple-quantum-transition spectral lines increases 
exponentially (~ 4") as the qubit number. The intensity of each multiple- 
quantum-transition spectral line therefore weakens exponentially (~ n/2 n ) 
as the qubit number increases, although all these multiple-quantum spectral 
lines are inphase due to the time-reversal symmetry Then due to noise in 
the detected NMR signal each of these multiple-quantum-transition spectral 
lines will become unobservable even when the qubit number is moderate if 
most of these multiple-quantum-transition spectral lines have different reso- 
nance frequencies (i.e., transition frequencies). Fortunately, due to the rota- 
tion symmetry in spin space there are very few multiple-quantum transitions 
with significantly different transition frequencies in the spin system and hence 
very few observable multiple-quantum spectral lines with significantly differ- 
ent resonance frequencies in the multiple-quantum spectrum. A large number 
of the inphase multiple-quantum spectral lines overlap with each other in the 
multiple-quantum spectrum due to that they have the same resonance fre- 
quency. As a result, all these inphase multiple-quantum spectral lines with 
the same resonance frequency become really a single multiple-quantum spec- 
tral line, and the intensity of the single multiple-quantum spectral line is re- 
ally the sum of all the intensities of these inphase multiple-quantum spectral 
lines. This intensity becomes so large that the noise in the detected NMR 
signal can not have any significant effect on the single multiple-quantum 
spectral line. As shown in Ref. [20], if now there are only the (n + 1) dif- 
ferent order multiple-quantum transitions in the n— qubit spin system, each 
with a significantly different transition frequency, then each of the (n + 1) 
different order multiple-quantum transitions is composed of the degenerative 
multiple-quantum transitions which have the number ~ (4 n /2 — 2™~ 1 )/(n-|-l) 
on average and its multiple-quantum-transition spectral intensity is propor- 
tional to the number (4 n /2-2 n-1 )/(n+ 1) on average. Here the total spectral 
intensity of the (4 n /2 — 2™~ 1 ) inphase multiple-quantum-transition spectral 
lines plus the spectral intensity of the longitudinal magnetization and spin 
order components is really equal to the total magnetization of the spin system 
[55] which can be observable in the conventional NMR experiments even for a 
very large qubit number. Obviously, the average intensity for each multiple- 
quantum-transition spectral line is inversely proportional to the qubit num- 
ber and some of these (n + 1) multiple-quantum spectral lines do not weaken 
exponentially as the qubit number increases. Thus, both the time-reversal 
symmetry and the rotation symmetry in spin space ensure that the output 
NMR multiple-quantum-transition spectral intensities can be efficiently de- 
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tected in the spin ensemble in the factoring algorithm. On the other hand, 
so far quantum entanglement has not yet been proven strictly to play a key 
role in speeding up a quantum computation, whereas the locally efficient and 
scalable factoring algorithm in the spin ensemble [20] shows that quantum 
entanglement could not be a key component to make quantum computation 
much more powerful than classical computation. Therefore, the unitary dy- 
namics quantum mechanically associated with the symmetry and structure of 
a quantum system could be the key component to power the quantum com- 
putation and hence is a general guidance to discover and develop new and 
efficient quantum algorithms in quantum computational science. The unitary 
dynamical principle quantum mechanically implies that the irreversible and 
non-unitary halting protocol should be modified to be reversible and unitary 
so that the incompatibility between the irreversible and non-unitary halting 
protocol and the unitary quantum computational process can be eliminated 
within the universal quantum computational models. 

The unitary dynamics quantum mechanically associated with the sym- 
metry and structure of a quantum system also plays a key role in discovering 
the scalable and efficient quantum algorithms in a pure-state quantum sys- 
tem. It is well known that there are 2 n basis states in the Hilbert space of 
an n— qubit spin system. These basis states may be chosen as the eigen- 
states of the component operator J z of the total spin angular momentum 
of the spin system. But due to the rotation symmetry in spin space of the 
n— qubit spin system many of these 2 n basis states are degenerative and have 
the same eigenvalue. Thus, the whole Hilbert space span by these 2 n eigen- 
states of the z— component operator J z may be divided into (n + 1) different 
state subspaces according to the rotation symmetry in spin space of the spin 
system [5a, 5c]. The quantum search space which is just the whole Hilbert 
space of the n— qubit spin system for the standard quantum search algorithm 
therefore is reduced to the largest state subspace among these (n + 1) state 
subspaces. Hence the standard quantum search algorithm is improved [23], 
since the largest state subspace is still much smaller than the whole Hilbert 
space. Moreover, it has been shown [23] that any unknown quantum state 
can be efficiently transferred to a larger state subspace from a state subspace 
in the Hilbert space of the n— qubit spin system, while the inverse process 
of the state transfer is generally harder to be carried out. This general rule 
not only is useful for solving the quantum search problem but also helpful 
for understanding deeply non-equilibrium processes of a quantum ensemble 
from the viewpoint of the unitary quantum dynamics instead of the conven- 
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tional probability theory. As shown in Ref. [24], this general rule is also 
closely related to that there exists a computational-power difference between 
a unitary evolution process and its inverse process in a quantum system. Of 
course, this inverse unitary process may exist or may not in the quantum 
system. A direct extension of the idea that the symmetric property and 
structure of a quantum system could help the unitary quantum dynamics to 
solve the quantum search problem is to exploit further the symmetric prop- 
erty and structure of a general group such as a cyclic group in a quantum 
system to help the unitary quantum dynamics to discover and develop new 
and efficient quantum search algorithms [24]. With the help of the reversible 
and unitary halting protocol based on the state-locking pulse field and the 
property and structure of a cyclic group in a quantum system it has been 
shown by the unitary quantum dynamics that a universal quantum computer 
could be enough powerful to solve efficiently the quantum search problem in 
a cyclic-group state space [24]. Here the important point to arrive at this 
conclusion is that the halting protocol of the universal quantum computa- 
tional models is available and may be made reversible and unitary for the 
quantum search process if the initial state of the quantum search process is 
limited to a single computational basis state. 

All these conventional halting protocols [3, 9, 10, 11, 12] of the univer- 
sal quantum computational models generally can not be used to construct 
an efficient quantum search algorithm based on the unitary quantum dy- 
namics [24]. This is because they are either irreversible and non- unitary 
or dependent sensitively upon initial states of the quantum computational 
process under study. For example, in the reversible computational model 
the step number of computational process in the reversible halting proto- 
col [54] is dependent sensitively upon the initial state if the output state is 
fixed in the reversible halting protocol, while in the reversible and unitary 
halting protocol [11] the output state is dependent sensitively upon the ini- 
tial state if the step number of computational process is fixed. Thus, both 
the reversible and unitary halting protocols are not suitable to construct an 
efficient quantum search process based on the unitary quantum dynamics. 
Only the specific reversible and unitary halting protocols that are based on 
the state-locking pulse field [24] could be useful for solving efficiently the 
quantum search problem in the cyclic group state space. This is because 
neither the output state nor the step number of computational process in 
such a reversible and unitary halting protocol is dependent sensitively upon 
any initial state of the computational process and this is the key point for the 
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reversible and unitary halting protocol to be useful for solving efficiently the 
quantum search problem. While the reversible and unitary halting protocol 
is the key component of the quantum search process to solve efficiently the 
quantum search problem in the cyclic group state space, the state-locking 
pulse field plays a key role in constructing such a reversible and unitary 
halting protocol. This is because the state-locking pulse field could make 
both the output state and the step number of computational process in the 
reversible and unitary halting protocol insensitive to any unknown marked 
state of the quantum search problem. The computational complexity for the 
quantum search process in the cyclic group state space [24] could be mainly 
dependent on the performance of the state-locking pulse field used in the 
quantum search process. The reversible and unitary halting protocol not 
only plays an important role in solving efficiently a quantum search problem 
in the cyclic group state space but also has an extensive application in quan- 
tum computation. This is one of the key components to realize a universal 
quantum computer to replace fully a classical computer in future. Therefore, 
it is necessary to design a good-performance unitary quantum control unit 
(or circuit) to simulate faithfully and efficiently the reversible and unitary 
halting protocol of the universal quantum computational models. A unitary 
quantum control unit that consists of a trigger pulse, a state-locking pulse 
field, and a control state subspace has been proposed to simulate faithfully 
and efficiently the reversible and unitary halting protocol [24]. The key com- 
ponent of the unitary quantum control unit is the state-locking pulse field. A 
state-locking pulse field is able to keep a desired state almost unchanged in a 
unitary form for a long time in a quantum computational process, and this is 
the reason why the output state of the reversible and unitary halting protocol 
based on the state-locking pulse field does not depend sensitively upon any 
initial states. How to design a good-performance state-locking pulse field is 
a challenge and also an important research subject in quantum computation 
in future. In general, a general state-locking pulse field may be dependent 
upon the time variable and the space variables and even the quantum field 
variables. A state-locking pulse field generally could be a sequence of time- 
and space-dependent electromagnetic pulse fields such as the laser pulses and 
could also contain any time- and space-dependent potential fields. 

In this paper a new quantum program and circuit is constructed explicitly 
to reduce the quantum search space which may be generally a cyclic-group 
state space or the Hilbert space of an n— qubit quantum system. This quan- 
tum program and circuit is much simpler than that one in the previous paper 
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[24] , where in order to show clearly that an ideal universal quantum computer 
could be powerful enough to solve efficiently the quantum search problem in 
the multiplicative-cyclic-group state space the quantum program and circuit 
used to reduce the quantum search space is designed in a more complex 
form. It can also be used to construct further a quantum search process 
to solve efficiently the unsorted quantum search problem in a general cyclic 
group state space or the Hilbert space of an n— qubit quantum system. The 
quantum program and circuit may be divided into two almost independent 
units: the unitary quantum computational unit and the unitary quantum 
control unit. The quantum control unit simulates efficiently the reversible 
and unitary halting protocol based on the state-locking pulse field, while the 
quantum computational unit is responsible for the reduction of the quantum 
search space. In the paper the quantum program and circuit is first analyzed 
completely. Then the basic properties of an ideal state-locking pulse field are 
described in detail, and it is shown how the quantum control unit simulates 
faithfully and efficiently the reversible and unitary halting protocol. The 
basic principles to construct a general state-locking pulse field and simulate 
efficiently the reversible and unitary halting protocol are suggested and then 
explained in detail. It is proposed in the paper that a simple atomic physical 
system which consists of an atomic ion or a neutral atom in the double-well 
potential field is used to realize the reversible and unitary halting protocol. 
In the atomic physical system a generalized state-locking pulse field used to 
build up the reversible and unitary halting protocol is also constructed ex- 
plicitly. 

2. The reversible and unitary halting protocol and the state- 
locking pulse field 

As an important application the reversible and unitary halting protocol 
may be used to build up the reversible and unitary quantum program and 
circuit to compute some mathematical functions. Here by solving a simple 
problem given below one may illustrate how to use the reversible and unitary 
halting protocol to solve a general mathematical problem in quantum com- 
putation. Suppose that given a periodic function f(x) = f(x + xt) with the 
period xt and the integer variable x — 0, 1, ...,xt — 1, there is a computa- 
tional circuit Uf to compute the functional value f(x + l) from the functional 
value f(x) for any integer x — 0, 1, xt~ 1, where the functional value f(x) 
is a distinct integer for every distinct integer x for < x < x T . The func- 
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tional operation Uf could be expressed as Uff(x) = f(x + 1) [13, 14] for 
x = 0, 1, ...,Xt — 1. Now given an unknown functional value f(x ) with the 
unknown integer xq G {0, 1, ...,xt — 1}, one wants the unknown functional 
value f(xo) to be changed to the desired functional value f(xf) with the 
known integer Xf G {0, 1, ...,xt — 1}, e.g., f(xf) = or 1. A simple scheme 
to solve this simple problem is that one computes one-by-one the functional 
value J'(xq + x) for the integer x — 0,1, Xt — 1 by the computational cir- 
cuit Uf and checks the functional value for each computing step, and when 
the functional value f(xo + x) is found to be equal to the desired functional 
value f(xf) the computational process is halted. Then the final result of the 
computational process is clearly the desired functional value f(xf). Though 
this problem is very simple from the viewpoint of the computational science, 
it is surprising that the scheme that can solve efficiently this simple problem 
in a reversible and unitary form could also be used to solve efficiently the 
quantum search problem. Note that given any initial value Xo and the final 
value Xf there is the unique integer x within the range < x < xt such 
that f(x + x) — f(xf). Now a classical computational program Q d for the 
scheme could be written down as 

n h = 

f(x) = f(x ) 
For i — 1 to xt 
If f(x) = f(xf) then n h = 1 

while nh — 1, halting 
else f(x) — > f(x + 1) end if 
end for 

where rih = or 1 is the halting bit that is used to indicate when the program 
terminates. As shown in the program, when the halting bit value — 1 the 
program terminates. The program consists of Xt cycles with the cyclic index 
% — 1, 2, xt- Evidently, this simple program outputs the desired functional 
value f(xf) no matter what the initial functional value f(xo) is with the 
possible integer xo = 0, 1, xt — 1. Here the unknown functional value f(xo) 
may be stored in memory on the computer in advance or is input from outside 
the program. The program first checks whether the initial functional value 
f(x ) equals f(xf). If f(xf) = f(x ), the halting bit value n h is changed to 1 
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from the initial value and then program terminates due to — 1, and the 
output result is f(xf); otherwise the program computes the functional value 
f(xo + 1) and checks again whether the functional value f(x + 1) is f(xf) 
or not. This computing process is repeated until the computed functional 
value f(xo + i) is found to be equal to f(xf), here (xo + i)modxT = Xf. 
Then the halting bit value — is changed to 1 so that the program is 
halted. Thus, the output functional value for the computational program Q c i 
is always f(xf) no matter what the initial functional value f(xo) is. This 
is the important property of the classical halting protocol. This property is 
also very important for the reversible and unitary halting protocol. 

If the classical computational program Q d would be reversible and uni- 
tary, then it could be suitably used to build up an efficient quantum search 
process to solve the quantum search problem in a cyclic-group state space, 
and such a quantum search process would be much simpler than that one 
in the previous paper [24]. For the quantum search problem in the multi- 
plicative cyclic-group state space S(C p -i) [24] the periodic function f(x) 
in the program Q c i may be chosen as the modular exponential function 
fk( x ) = (g Mk ) x modp with the order (or period) xt = and the integer 
variable x — 0, 1, ...,rrik — 1 for k — 1,2, ...,r. Then the corresponding re- 
versible functional operation U / may be defined as the unitary cyclic-group 
operation U g M h : fk(x) — > fk(x + 1). Here p is a prime and g the gen- 
erator of the multiplicative cyclic group C p _i, and the group operation of 
a multiplicative cyclic group is the modular multiplication operation. The 
multiplicative cyclic group C p _i has the order p — 1 = p^p^ 2 ■■■P% r , where 
Pi, p 2 , p r are distinct primes and the exponents ai, a 2 , a r > 0. Here 
the integer = p°£ and p — 1 = m^Mfc for k = 1, 2, r. Any pair of the 
integers m 8 and rrij are coprime to each other for 1 < % < j < r, and for 
convenience usually set mi < m 2 < ... < m r . Obviously, every functional 
value is distinct for the modular exponential function fk(x) which satisfies 
fk(x) > 1 for x — 0, 1, ...,rrik — 1. The functional states {\fk(x))} really form 
the multiplicative-cyclic-group state subspace S(rrik) = {\{g Mk ) Sk modp), 
5^ = 0, l,...,m k — 1} of the factor cyclic subgroup C «-k of the multiplica- 
tive cyclic group C p _i = {g s modp} = C p <n x C p «2 x ... x C p «r, whereas the 
multiplicative-cyclic-group state space S(C p -i) with dimension p—1 is the di- 
rect product of the multiplicative-cyclic-group state subspaces {S(mk)} with 
dimensions {m k } : 

5(C7 p _i) = S( mi ) (g) S(m 2 ) (g) ... (g) S(m r ). 
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If the quantum search problem is solved in the additive-cyclic-group state 
space S(Z p _i) with dimension p — 1 = m 1 m 2 ...m r , then the periodic func- 
tion f(x) may be taken as the modular function fk{x) = xmodm^ with the 
period Xt = rn^ and x = 0,l,...,mk — 1 for k = 1,2, ...,r. Every modu- 
lar functional value fk(x) is distinct and fk(x) = 0, 1, ...,m,k — 1. Obviously 
the modular functional states {| /&(£))} form really the additive-cyclic-group 
state subspace S(Z mk ) = {|sfemodmfc), = 0, l,...,m.k — 1}. The additive 
cyclic group Z p _\ = {0, 1, — 2} may be decomposed into the direct sum of 
the factor additive cyclic subgroups {Z mk } : Z p ^ = Z mi Z m2 ... Z mr , 
here the group operation of an additive cyclic group is the modular addition 
operation. The additive-cyclic-group state space S(Z p _ 1 ) then may be de- 
composed into the direct product of the factor additive-cyclic-group state 
subspaces {S(Z mk )} with dimensions {m^} : 

= S(Z mi ) (g) S(Z m2 ) (g) ... (g) S(Z mr ). 

Though the 2 n — dimensional Hilbert space S(Z2") of an n— qubit quantum 
system may be thought of as an additive-cyclic-group state space under the 
modular addition operation (mod2 n ), for the Hilbert space S(Z 2 ™) there is 
not a direct-sum decomposition like S(Z p _i) above with r > 1. However, 
the Hilbert space S(Z 2 n) may be decomposed into the direct product of n 
2— dimensional additive-cyclic-group state subspaces {S , (Z 2 )}, 

S(Z 2n ) = S(Z 2 ) (g) S(Z 2 ) (g) ... (g) S(Z 2 ). 

If now the periodic function f(x) is chosen as the modular function fk{x) = 
a; mod 2 with x — 0, 1 for k — 1, 2, n, then the quantum search problem in 
the 2™— dimensional Hilbert space of the n— qubit quantum system may also 
be solved just like that one in the additive-cyclic-group state space S(Z p -i) 
or in the multiplicative-cyclic-group state space S(C p -i) [24]. The symmetric 
properties and structures of both the cyclic groups C p _i and Z p _\ have been 
suggested to help the unitary quantum dynamics to solve the quantum search 
problem in the cyclic-group state spaces S(C p -i) and S(Z p -i), respectively 
[24]. 

Obviously, the classical computational program Q c \ is not reversible due 
to the irreversible operations in the program and especially due to the fact 
that the halting operation within the program is not reversible. However, 
in order that it can be used to construct a quantum search process to solve 
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the quantum search problem in the cyclic-group state space the whole pro- 
gram Qd including the halting operation must be made reversible and uni- 
tary. Generally, a classical irreversible computational program may be made 
reversible in the frame of the reversible computational model [1] with the 
help of the reversible mathematical logic gates and especially the reversible 
operation of a general function mathematically [1], the universal quantum 
gates and especially the conditional quantum gates [25], and other unitary 
operations quantum mechanically [26]. One can also construct directly a 
unitary quantum circuit equivalent to a classical computational program on 
the universal quantum circuit model [4]. It has been shown that there is 
an equivalent quantum circuit to a given reversible quantum program on 
the universal quantum Turing machine [27]. However, these conventional 
methods by which an irreversible classical computational program is made 
reversible and unitary could not be always suitable for transforming the ir- 
reversible halting protocol to the reversible and unitary one in the universal 
quantum computational models. First, there could not be a universal halting 
protocol in the universal quantum computational models. Second, although 
there could be a halting protocol in the universal quantum computational 
models when the input state of a quantum computational process is limited 
to a single basis state, this halting protocol could not be thought to be re- 
versible in the sense that an information loss occurs in the halting protocol 
due to the non-unitary quantum measurement. For example, by the quan- 
tum measurement one could know what time the halting state |0) is changed 
to the state |1) and hence obtains the information of the instant of time at 
which the quantum system arrives in the end state of the quantum compu- 
tational process, implying that the quantum computational process loses the 
information of the instant of time. However, even if the current quantum 
measurement could not be non-unitary according to quantum mechanics [5b] 
as the measured base now is a single basis state of the measurement opera- 
tor and such an information loss could not have a significant effect on some 
quantum computational processes, this information is not available accord- 
ing as the mathematical logic principles of the quantum search problem and 
can not be yet allowed to use if the halting protocol is used in solving the 
quantum search problem based on the unitary quantum dynamics. This is 
because a quantum computational process obeys the quantum physical laws 
and is compatible with the used mathematical logic principles. Therefore, 
the classical program Q d may be made reversible and unitary according as 
these conventional methods except the irreversible halting operation within 
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the program. Of course, in the case that any initial state is restricted to be 
a single basis state the halting operation may also be made reversible by the 
conventional methods [11, 54] for the conventional computational tasks other 
than the current tasks. In the reversible computational model the reversible 
halting operation in the program Q c i may be achieved in such a way [54] 
that after the computational process from the initial functional value f(x ) 
to the output functional value f{xf) is done the conditional halting oper- 
ation is executed and then the inverse computational process is performed 
from the functional value f(xf) back to the original functional value f(x ), 
and the cyclic process consisting of the computational process and inverse 
computational process repeats incessantly. Evidently, this extra inverse com- 
putational process is dependent upon the cyclic index i of the program as 
a different initial functional value f(xo) is changed to the functional value 
f(xf) in a different cycle in the program. One could also use the halting pro- 
tocol [11] to achieve the reversible halting operation for the program. When 
the program arrives at the output functional state (corresponding to f(xf), 
see below) the conditional halting operation is executed and then an extra 
unitary evolution process starts, that is, the state of the auxiliary register 
starts to evolve [11]. If now the step number of the program Q c i is fixed (for 
example the step number may be set to xt), then a different initial func- 
tional state (corresponding to f(x )) will result in a different output state of 
the auxiliary register. Thus, the total output state of the program is really 
dependent upon the initial functional state. One therefore concludes that if 
any of the two reversible halting protocols [11, 54] is used in the program, 
then the program is dependent sensitively upon the initial functional state. 
Such a reversible program is not suited as a component of the quantum 
search process based on the unitary quantum dynamics. It is necessary to 
use the specific method to make the halting operation within the program 
Qd reversible and unitary if one wants to use the program to construct a 
quantum search process. This specific scheme to make the halting protocol 
reversible and unitary may be stated below. As the first point of the scheme, 
in order that the halting protocol is reasonable from the viewpoint of the 
mathematical logic principle any initial state is limited to a single basis state 
for the halting protocol. The cost for this point is that the quantum parallel 
principle could become less important. As the second point, the halting pro- 
tocol should not contain any quantum measurement so as to keep away from 
any irreversible and non-unitary process and avoid any information loss of the 
quantum system. Generally, the unitary quantum dynamics avoids using any 
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non-unitary quantum measurement as its quantum operation within quan- 
tum computational processes. This is quite different from the early proposals 
that the quantum measurement may also be used as a quantum operation to 
build up quantum circuits and algorithms [13, 28]. If the halting operation 
within the classical program Q c i is achieved by the brute-force method, then 
the end state (i.e., the output state) of the classical program is independent 
of any initial functional value J(xq). But the brute-force method to stop the 
program from the outside could cause the whole program irreversible. There- 
fore, for the third point of the scheme the brute-force method is replaced with 
the unitary operation conditionally depending upon the halting state to stop 
the program [11, 12, 54]. However, all these three points above in the scheme 
can not ensure that the end state of the program which could include the 
computational state, halting-qubit state, and auxiliary state is independent 
of any initial functional state, that is, all these three points can not lead to 
the important property of the halting operation that the output state is inde- 
pendent of any initial state. It is of crucial importance to make the end state 
of the program independent or almost independent of any initial functional 
state. This is because whether or not the quantum search process based on 
the reversible and unitary halting protocol is efficient is mainly dependent 
upon this point. Thus, as the fourth point, the end state of the program is 
locked by the state-locking pulse field so that it is not dependent sensitively 
upon any initial functional state. This scheme which consists of the four 
points above has been used in the previous paper [24]. However, this scheme 
is not intuitive to understand the reversibility and unitarity of the reversible 
and unitary halting protocol mainly due to that it is hard to understand how 
the state-locking pulse field is capable of keeping a quantum state almost 
unchanged in a unitary form for a long time in a quantum computational 
process. 

It was proposed in the previous paper [24] that a unitary quantum con- 
trol unit which consists of a trigger pulse (Pt), a state- locking pulse field 
(Psl), and a two-state control state subspace is used to simulate efficiently 
the reversible and unitary halting protocol. This quantum control unit does 
not use any non-unitary quantum measurement as its component. The state- 
locking pulse field which is the key component of the unitary quantum control 
unit is used to keep the desired state almost unchanged in a unitary form 
for a long time in the quantum computational process under study. Now 
by using a similar unitary quantum control unit to simulate faithfully and 
efficiently the halting operation within the classical program Q d and with 
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the help of the reversible mathematical logic operations [1], the universal 
quantum gates and especially the conditional quantum gates [25], and any 
other unitary quantum operations, a reversible and unitary quantum compu- 
tational program Q c and its equivalent quantum circuit which correspond to 
the classical program Q c i can be constructed explicitly in the frames of the 
reversible computational model [1] and the universal quantum circuit model 



[4], respectively, and they may be represented intuitively by 

State-Locking Pulse Field (P SL ) -ON (PI) 

K> = |0> (P2) 

\h) = |0> (P3) 

\fr(x)) = \f r (x )) (P4) 

For % — 1 to m r (P5) 

If \f r ( x )} = |i) then U c h : \b h ) -> \b h + 1) end if (P6) 
While |/ r (x)> = |1), Do U c h : \n h )\f r (x)) = |0)|1) - |0)|0), 

P t : \n h )\f r (x)) = |0)|0) - | Cl )|0), P SL : |ci) - |c 2 ) (P7) 

If \b h ) = |0) then U f : \f r (x)} -> \f r (x + 1)) end if (P8) 

end for (P9) 

State-Locking Pulse Field (P SL ) ■ OFF (P10) 



Note that the halting bit n h and the function f(x) in the classical program 
Qd have already been replaced with the halting state \nh) and the func- 
tional state \f(x)) in the quantum program Q c , respectively. In the quantum 
program Q c the functional state \f(x)) is set to the modular exponentia- 
tion state \f r (x)) with the period xt = m r of the multiplicative-cyclic-group 
state subspace S(m r ) = {\f r (x)) = \{g Mr ) x modp); x = 0, l,...,m r — 1} 
and the desired functional state \f r (%f)) se t to the state |1). Note that the 
multiplicative-cyclic-group state subspace S(m r ) does not contain the state 
|0) because the functional value f r (x) > 1 for any integer x — 0,1, ...,m r — 1. 
Owing to f r (x) 7^ here the state \f r (x)) = |0) means that the state in the 
register of the functional state \f r (x)) takes the state |0) rather than that the 
functional value f r (x) = 0. The unitary functional operation Uf is defined by 
Uf\f r (x)) = |/ r (a; + l)) for \f r (x)) G S(m r ), here Uf is really the cyclic-group 
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unitary operation U g M r of the factor cyclic subgroup C p *r of the multiplicative 
cyclic group C p -\. If the state \f r (x)) = |0), then Uf\f r (x)) = Uf\0) = |0) as 
the state |0) does not belong to the multiplicative-cyclic-group state subspace 
S(m r ). More generally the functional operation Uf satisfies Uf\g(x)) = \g(x)) 
if the state \g(x)) is not in the state subspace S(m r ). Though here the quan- 
tum program Q c and its equivalent quantum circuit are designed for the 
multiplicative-cyclic-group state subspaces {S'(mfc), k = 1,2, ...,r} (the in- 
dex r may be replaced with k in the quantum program Q c ) , similar quantum 
programs and circuits can also be constructed for the additive-cyclic-group 
state subspace S(Z mk ) = {\fk(x)) = \xmodrrik); x = 0,l,...,m,k — 1} for 
k — 1,2, ...,r and for a general periodic-function state space. The quantum 
program Q c contains mainly m r cycles with the cyclic index % — 1 to m r and 
the state-locking pulse field {Psl)- 

The quantum program Q c consists of ten statements, which are denoted 
conveniently as the statement PI, statement P2, statement P10, respec- 
tively. In particular, the statement P7 consists of the conditional unitary 
operation U%, the trigger pulse P t , and the state-locking pulse field Psl- The 
quantum program Q c is more complex than the classical one Q c i mainly due 
to the statement P7. The statement P7 contains the unitary quantum con- 
trol unit that simulates efficiently the reversible and unitary halting protocol. 
The unitary quantum control unit consists of the trigger pulse P t , the state- 
locking pulse field Psl, and the control state subspace S(C) = {|ci), |c 2 )}. 
Any one of these two states |ci) and \c2) of the control state subspace is 
different from the initial halting state \rih) = |0). Actually the three states, 
i.e., the halting state \rih) = |0) and these two states |ci) and \c2), should be 
orthogonal to each other and belong to the same register which is named the 
halting register here. These three states are also called the halting-register 
states or briefly the halting states. It will be seen in next section that the halt- 
ing register may be generated from the simple physical system of an atomic 
ion or a neutral atom in the double-well potential field. The physical system 
of the halting register may also be called the quantum control system whose 
Hilbert space contains the control state subspace S(C) and probably other 
relevant states. Both the trigger pulse P t and the state-locking pulse field 
Psl are applied only to the physical system of the halting register. Actually, 
in the quantum program Q c the state- locking pulse field Psl is applied only 
to the control state subspace S(C). The trigger pulse P t acts on the initial 
halting state \nh) = |0) to convert conditionally it into the state \c\) of the 
control state subspace only if the state |/ r (x)) in the register of the functional 
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state is the state |0). Thus, here the unitary transformation U t during the 
trigger pulse P t could be defined simply by 

C/ t :|n h )|/ r (a;)) = |0)|0)^|c 1 )|0). 

Both the trigger pulse and the state-locking pulse field are very important 
in the unitary quantum control unit and will be discussed separately in de- 
tail later. In order that the halting quantum bit {\nh),Uh = 0,1} is sep- 
arated from other qubits of the computational state subspace in the quan- 
tum program an extra quantum bit named the branch-control quantum bit 
{\bh),bh = 0,1} is used to control directly the functional operation of the 
function f r (x) in place of the halting quantum bit. Therefore, there is the 
unitary conditional functional operation Uf for the functional state \f r (x)) 
defined by 

TTC lh V| f _ / |0>|/ r (x+l)> if b h = 

U f \b h )\f r (x)) - ^ MMx)) iibh ^ Q 

This definition shows that the conditional functional operation Uj is applied 
only to both the functional state \f r (x)) and the branch-control state \bh)- 
If the branch-control state \bh) is the state |0) the functional operation Uf 
changes the functional state \f r (x)) to another functional state \f r (x + 1)), 
otherwise the operation U / does not change the functional state. Therefore, 
the branch-control state \bh) can control conditionally the action of the func- 
tional operation U / upon the functional state. When the branch-control state 
\bh) is changed from the state |0) to the state |1) the functional operation Uf 
is halted to act on the functional state \f r (x)) even though the conditional 
functional operation Uf continues to apply to the functional state. In the 
quantum program the halting quantum bit {1^^)} is designed to control the 
branch-control quantum bit {|&^)} indirectly and hence controls ultimately 
the action of the functional operation Uf on the functional state. The condi- 
tional unitary transformation U£ in the quantum program could be defined 
simply by 

mh)lf(r)) _f \(b h + l)mod2)\f r (x)), if|/ r (x)> = |l> 

U b \WrW) ~ \ \ bh )\f r ( x) ), if |/ r (x)>^|l> 

and the conditional unitary transformation U% : \n h )\f r (x)} = |0)|1) <-> |0)|0) 
is defined explicitly by 

|0)|0), if \f r {x)) = \l) and K> = |0> 
U c h \n h )\f r {x)) = { |0)|1), if |/ r (a;)) = |0) and \n h ) = |0) 
\ n h)\fr( x ))i otherwise 
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The detailed analysis for the quantum program Q c is given below. The 
state-locking pulse field Psl is first switched on at the beginning of the quan- 
tum program (see the statement PI in the program) and could be switched off 
(or partly switched off) after the quantum program finished (the statement 
P10). It is mainly used to manipulate the states \c±) and |c 2 ) and lock the 
state |c 2 ) of the control state subspace S(C). Both the initial halting state 
\nh) (statement P2) and the initial branch-control state \bh) (the statement 
P3) are simply set to the state |0). The initial functional state \f r {%o)) (the 
statement P4) could be unknown and may be stored in the memory of the 
quantum computer in advance or is input from outside the quantum program. 
There are m r possible initial functional states {\f r (x )), x o — 0, 1, m r — 1} 
at most. The quantum program first checks whether or not the initial func- 
tional state \ f r (xo)} is the desired functional state \f r (xf)) = |1) (the state- 
ment P6). Then there are two possible cases to be considered, that is, ei- 
ther |/ r (xo)) = \fr(xf)} or \f r (x )} ^ \f r (xf)}. Consider the first case that 
\fr(%o)} — \fr(xf)) — Since the state \f r (x )) = \f r (xf)) the branch- 
control state \bh) = |0) is first changed to the state |1) by the conditional 
unitary operation (the statement P6). Then the desired state \f r (xf)) is 
changed conditionally to the state |0) (the statement P7) by the conditional 
unitary operation U% due to that the halting state \nh) now is the state |0), 
the initial halting state \n h ) = |0) then is changed to the state |ci) of the 
control state subspace by the trigger pulse P t and then the state |ci) to an- 
other orthogonal state |c 2 ) of the control state subspace, since then the state 
|c 2 ) is kept unchanged by the state-locking pulse field Psl (the statement 
P7). When the quantum program executes the statement P8 the conditional 
functional operation Uj will not have a net effect on the functional state 
\f r (x)) according as the definition of the operation Vj because the branch- 
control state \bh) now is the state |1) and the state \f r (x)) = |0), although 
now the operation Uj is still applied to the whole quantum system of the 
quantum computer. This shows that the functional operation Uf acting on 
the functional state \f r (x)) is really halted after the initial branch-control 
state |0) is changed to the state |1). Now the quantum program finished the 
first cycle with the index % — 1. It then returns and executes the statement 
P5 of the second cycle with the index % — 2. When the quantum program 
executes the statement P6 in the second cycle, the current branch-control 
state \bh) = |1) keeps unchanged under the action of the operation U% as the 
current state \f r (x)) = |0). Since the current state \nh)\f r (x)) = | c 2 ) 1 0) , that 
is, the state \rih) ^ |0), the unitary operation U% does not have a real effect 
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on the quantum system when the quantum program executes the statement 
P7, and since now the halting-register state \rih) is the state \c2), that is, the 
state \nh) is neither |0) nor |ci), the trigger pulse Pt does not yet have a real 
effect on the quantum system. The key point in the quantum program is 
that the state |c 2 ) of the control state subspace has been locked by the state- 
locking pulse field P$l since the quantum program executes the statement 
P7 in the first cycle (i = 1). Hence the state \c2) is not yet changed when the 
statement P7 is executed in the second cycle. Actually, the state \c2) may be 
kept unchanged till the end of the quantum program after the statement P7 
was executed in the first cycle. If now the program continues to execute the 
rest statements and even run till the end of the program, then all these states 
\ n h) — \ c 2)) \bh) — and \f r (x)) = |0) of the whole quantum system of the 
quantum computer are still kept unchanged due to the fact that the halting- 
register state \rih) is kept in the state \c2) by the state-locking pulse field. 
Note that the conditional functional operation Uj is applied continuously to 
the quantum system of the quantum computer even after the branch-control 
state \bh) = |0) is changed to the state |1), which leads to that the functional 
operation Uf acting on the functional state \f r (x)} is halted. Of course, 
in this case the conditional functional operation Uj has not a net effect on 
the functional state \f r (x)). This process is repeated from the second cycle 
{i = 2) to the end {i = m r ) of the quantum program. The analysis above 
shows that if the functional state \f r {x)) is the desired state \f(xf)), then 
when the quantum program executes the statement P6 the branch-control 
state \bh) = |0) is changed to the state |1), and then on the statement P7 the 
initial halting state \rih) = |0) is changed to the state |ci) and then further to 
the state \c2) which is ultimately kept unchanged by the state- locking pulse 
field, since then the state \n h )\bh)\f r (%)) = | c 2 ) j 1) j 0) of the whole quantum 
system is kept unchanged till the end of the program. Therefore, the quan- 
tum program outputs the final state \nh)\bh)\f r (x)) = | C2 } 1 1 ) 1 } (the state 
\fr(x)) = jO) is easily changed to the desired state \f r (xf)) = Next con- 
sider the second case: the state \f r (xo)) 7^ \fr{xj))- If the initial functional 
state \ f r (x )) 7^ \ f r (xj)) — |1), then the initial branch-control state \bh) = |0) 
is not changed to the state |1) by the unitary operation U£ when the program 
executes the statement P6. Again due to that the state |/ r (^o)) 7^ |1) and 
|0) the unitary operation does not really act on the state \rih = 0)\f r (xo)) 
and both the trigger pulse P t and the state-locking pulse field Psl do not 
yet act on the initial halting state \rih) = |0) when the program executes the 
statement P7. Since the branch-control state \bh) = |0) the functional state 
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\f r (xo)) is changed to the state \f r (xo + 1)) by the conditional functional 
operation after executing the statement P8. Now the quantum program 
returns to execute the statement P6 of the second cycle after the cyclic index 
% — 1 is changed to % = 2 on the statement P5. Again the program first 
checks whether or not the functional state \f r (x + 1)) is the desired state 
\f r (xf)). Just like before there are also two possible cases, the first case is 
\f r (x +l)) = \f r (xf)) and the second \ f r (xo+l)) ^ \ f r (xf)). As shown above, 
for the first case \f r (xo + 1)) = \fr(xf)) the program will output the final 
state \rih)\bh)\f r (x)) = |c2)|l)|0). For the second case \f r (xo + 1)) ^ \fr(xf)) 
the functional state \f r (x + 1)) will be changed to the state |/ r (xo + 2)) 
at the end of the second cycle (i = 2) of the program (the statement P8). 
This process is repeated till the k— th cycle (m r > k > 1) when the func- 
tional state \f r (xo + k)) = \f r (xf)) at the end of the A;— th cycle. Here the 
index k is unique for < k < m r and k = corresponds to the earlier case 
\fr(xo)) = \fr(xf)). Now for the (k + 1)— th cycle the initial branch- control 
state \b h ) = |0) is first changed to the state |1) (the statement P6), following 
the statement P6 the functional state \f r (x + k)) = \f r (xf)) is changed to 
the state |0), then the initial halting state \rih) = |0) to the state |ci) by the 
trigger pulse P t and further to the state \c2) by the state-locking pulse field 
(the statement P7), and since then the state \c?) is kept unchanged by the 
state-locking pulse field till the end of the program. Therefore, the quantum 
program outputs finally the state \n h )\bh)\f r {x)) = |c 2 )|l)|0). This shows 
that after executing quantum program Q c the initial state \rih) \bh) \f r (x)) = 
|0)|0)|/ r (xo)) is always transferred to the output state |c2)|l)|0) no matter 
what the initial functional state |/ r (^o)) is with xq = 0, l,...,m r — 1. Note 
that the quantum program Q c is reversible and unitary because all these op- 
erations of the quantum program are reversible and unitary. One therefore 
concludes that by the unitary quantum program Q c different initial states 
{|0)|0)l/r(^o))} are transferred to the same output state |c2)|l)|0) and hence 
the output state of the quantum program is not dependent sensitively upon 
any initial states. However, the first part of the conclusion is apparently in 
conflict with the fact that different input states can not be completely trans- 
ferred to the same output state by a given unitary transformation. Therefore, 
the first part of the conclusion is expressed exactly as that different initial 
states {|0)|0)|/ r (:ro))} are transferred to the same output state |c 2 )|l)|0) in 
probabilities approaching infinitely 100% in theory by the unitary quantum 
program Q c . In theory there is only one initial state \0)\0)\f r (xo)) that may 
be completely transferred to the output state |c 2 )|l)|0) by the unitary quan- 



27 



turn program. As shown in next sections, due to the fact that the quantum 
program Q c is reversible and unitary the output state |c2)|l)|0) can be re- 
ally obtained from the initial state |0)|0)|/ r (xo)) only in a probability close 
to 100% rather than in the probability 100% for any initial functional state 
\f r (xo)) in a real physical system. 

The state-locking pulse field Psl plays a key role in the quantum control 
process that simulates efficiently the reversible and unitary halting protocol 
in the quantum program Q c . It is the state-locking pulse field that keeps 
the state \c2) of the control state subspace unchanged till the end of the 
quantum program after the functional state \f r {x)) is changed to the desired 
functional state \f r (xf)), while keeping the state |c 2 ) unchanged for a long 
time is the key step to achieve the reversible and unitary halting protocol. 
As pointed out before, the statement P7 of the quantum program which is 
involved in the state- locking pulse field Psl is mainly used to simulate the 
reversible and unitary halting protocol. The statement P7 really forms a 
unitary quantum control process (or unit). This process (or unit) is almost 
independent of the quantum computational process (or unit) to compute 
the desired functional state \f r (xf)) in the quantum program, but it really 
controls the quantum computational process (or unit). A quantum program 
and its quantum circuit may be generally divided into two parts, one part 
is the quantum computational unit (or process) and another the quantum 
control unit (or process). As an example, for the quantum program Q c the 
quantum computational unit is used to compute the desired functional state 
\f r {xf)), while the quantum control unit is used to perform the reversible 
and unitary halting protocol. As pointed out before, the unitary quantum 
control unit consists of the trigger pulse P t) the state- locking pulse field 
Psl, and the control state subspace S(C). Generally, the control state sub- 
space S(C) is different from the computational state subspace such as the 
multiplicative-cyclic-group state subspace S(m T ) in the quantum program, 
but they all belong to the Hilbert space of the whole quantum system of 
the quantum computer. The simplest control state subspace is a two-state 
subspace such as S{C) = {|ci), |c 2 )}, but generally a control state subspace 
is not restricted to be a two-state subspace in the quantum program. The 
trigger pulse P t could be used for communication between the control state 
subspace and the computational state subspace. It generally triggers a time- 
and space-dependent unitary evolution process in the quantum control sys- 
tem with the control state subspace S(C). The unitary transformation for 
the trigger pulse may be defined explicitly. For example, in the quantum 
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program Q c the initial halting state \rih) = |0) is converted into the state \ci) 
and vice versa by the trigger pulse P t only if the state \f(x)) is the state |0). 
Then the unitary transformation for the trigger pulse Pt is defined simply 
by Ut, as shown before. Here the state \f r (x)) = |0) in the register of the 
functional state could be thought to be related to the computational state 
subspace, while the state |ci) is of the control state subspace. A different 
definition for the unitary transformation during the trigger pulse P t can be 
seen in next section. The unitary quantum control process that simulates the 
reversible and unitary halting protocol may be simply described below. Here 
it is first pointed out that the state- locking pulse field Psl is applied only 
to the quantum control system with the control state subspace S(C) in the 
quantum program. The state-locking pulse field is first switched on to apply 
to the quantum control system of the quantum computer at the beginning of 
the quantum program (see the statement PI), but it usually does not take a 
real action on the quantum control system at the beginning time. However, 
when the functional state \f r (x)) is changed to the desired state \f r (%f)) and 
then the initial halting state \nh) = |0) changed conditionally to the state 
|ci), the state-locking pulse field Psl which has been applying to the control 
state subspace since the beginning of the quantum program starts to take a 
real action on the states of the control state subspace. The state \c±) is first 
sent to the state |c 2 ) in the control state subspace under the state-locking 
pulse field. This process usually is a time- and space-dependent unitary evo- 
lution process. Then the state \c2) is kept unchanged by the state-locking 
pulse field to the end of the quantum program and circuit. Because now 
the branch-control state \bh) leaves the initial one and is kept in the state 
|1) unchanged due to that the state |c 2 ) is kept unchanged by the state- 
locking pulse field the computational process is halted conditionally and the 
reversible and unitary halting operation therefore is achieved. According to 
this picture that the quantum control process simulates the reversible and 
unitary halting protocol the state- locking pulse field Psl is applied contin- 
uously to the quantum control system from the beginning to the end of the 
quantum program. If the quantum system of the quantum computer which 
includes the quantum control system now is acted on by a unitary operation 
such as one of the unitary operations U^, U^, Uf, and the trigger pulse P t of 
the quantum program, then actually it is acted on by both the state-locking 
pulse field and the unitary operation simultaneously. Then the state-locking 
pulse field could be designed in such a way that the effect of the state-locking 
pulse field on the quantum system can be negligible during the period of the 
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unitary operation applying to the quantum system. For example, the unitary 
transformation on the quantum system is approximately equal to the single 
unitary transformation of the functional operation Uj when both the func- 
tional operation Uj and the state-locking pulse field Psl are applied to the 
quantum system simultaneously in the quantum program. This is because in 
this case the state-locking pulse field has a negligible effect on the quantum 
system. However, as shown in next section, there may also be another case 
that the unitary transformation during the trigger pulse Pt could be really 
generated by both the state-locking pulse field and the trigger pulse. Then 
in this case the contribution of the state-locking pulse field is not negligible. 
These general properties of an ideal state-locking pulse field could be used 
to measure the performance of a real state-locking pulse field used in the 
reversible and unitary halting protocol. 

The quantum control process that simulates the reversible and unitary 
halting protocol could be a single time-dependent unitary evolution process, 
but generally it may be a time- and space-dependent unitary evolution pro- 
cess of the quantum control system. However, if the quantum control process 
is restricted to be dependent only upon a single time variable, there could 
be a large drawback for such a quantum control process with a two-state 
control state subspace S(C). This can be explained in detail below. Suppose 
that the state |ci) of the control state subspace is generated completely from 
the initial halting state \n h ) = |0) by the trigger pulse P t at the instant of 
time toi in the %— th cycle of the quantum program (see statement P7). The 
instant of time toi is special in that the state-locking pulse field really starts 
to act on the control state subspace S(C) at the instant of time toi in the 
quantum program Q c . Here suppose that the state-locking pulse field has a 
negligible effect on the state \ci) during the trigger pulse P t . Evidently, a 
different initial functional state \f r (xo)) corresponds to a different time t 0i , 
and there are m r possible and different times {toi, i — 1, 2, m r } at most in 
the quantum program because there are m r different initial functional states 
{\fr(xo)), %o = 0, l,...,m r — 1} of the cyclic-group state subspace S(m r ). 
Suppose that ATj (1 < i < m r ) is the time period of the i— th cycle of the 
quantum program. Obviously, = to(i-i) + Alj for % = 1, 2, m r , where too 
may be defined as too — toi — ATi. For convenience, here the time period ATj 
is set to the same one AT C for every cycle i = 1,2, m r . Suppose that in the 
quantum program Q c the duration of the trigger pulse P t is 5t r and the du- 
ration is denoted as At during which the state \c\) is converted completely 
into the state |c 2 ) in the control state subspace by the state- locking pulse 
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field. In order to show that a quantum control process that is restricted to 
be dependent only upon a single time variable is not a good one to simulate 
efficiently the reversible and unitary halting operation there are two possible 
situations to be investigated in the quantum control process. Consider the 
first situation. Note that in the i— th cycle of the quantum program the 
initial halting state \n h ) = |0) is changed completely to the state |ci) by the 
trigger pulse P t in the time interval from the time — 8t r to the time t^. 
This means that the functional state \f r (x)) is changed to the desired state 
\f r (xf)) at the end of the (i — 1)— th cycle, and during the period from the 
time toi to the time t Qi + At the state \ci) is converted into the state \c2), 
and from the time t 0i + At on, the state |c 2 ) is kept unchanged to the end of 
the quantum program by the state-locking pulse field Psl- Evidently, before 
the instant of time io(j+i) — St r i n the (i + 1)— th cycle the state \c\) must be 
completely converted into the state \c2) and since then the state \c2) is kept 
unchanged by the state-locking pulse field, otherwise it is possible that the 
residual state |ci) could be changed back to the initial halting state \n h ) = |0) 
by the trigger pulse P t during the period from the time £o(i+i) — dt r to the 
time to(i+i) i n the (i + 1)— th cycle. Therefore, for the first situation the quan- 
tum control process requires that the quantum control system be completely 
in the state \c2) in the time interval between to(«+i) — 8t r and io(i+i) in the 
(i + 1) — th cycle of the quantum program. Now consider the second situation. 
There is also another possibility that unlike the first situation the functional 
state \f r (x)) is converted into the desired state \f r (xf)) in the i — th cycle 
instead of the (i — 1)— th cycle, because there are different initial functional 
states {\f r (xo))} in the quantum program. Then during the period from the 
time to(j+i) — $t r to the time to(«+i) in the {i + 1)— th cycle the initial halting 
state \rih) = |0) is changed completely to the state |ci) by the trigger pulse 
P t . This shows that for the second situation the quantum control system is 
in the state \c\) at the time t (i+i) in the {i + 1)— th cycle. Then it can be 
seen from the two possible situations that at the instant of time £o(i+i) in 
the {i + 1)— th cycle of the quantum program the quantum control system is 
either completely in the state |ci) for the second situation or completely in 
the state \c2) for the first situation. Note that these two states |ci) and \c2) 
of the control state subspace are orthogonal to one another. Now one con- 
siders the {i + 2)— th cycle of the quantum program. The quantum program 
requires for any one of the two possible situations that the quantum control 
system be completely in the state \c2) so as to avoid any real effect of the 
trigger pulse on the quantum system during the period of the trigger pulse 
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between the time io(i+2) _ St r and the time £0(1+2) in the (i + 2)— th cycle. 
Evidently, given a state-locking pulse field and a trigger pulse as well as other 
unitary operations in the quantum program there is the same time evolution 
propagator U(t ( i+2 ) — <5t r , t (i+i)) of the quantum control system during the 
period from the time t (i+i) m the (i + 1)— th cycle to the time io(»+2) — dt r 
in the (i + 2)— th cycle no matter that the quantum control system is in the 
state \ci) for the second situation or in the state |c 2 ) for the first situation at 
the instant of time io(f-t-i)- Thus, there are two possibilities to be considered. 
The first one is that the state of the control state subspace is the state |c 2 ) 
at the instant of time to(«+i) m the (i+ 1)— th cycle, which corresponds to the 
first situation above. Since during the period between the time to(«+2) ~ St r 
and the time £o(i+2) m the (i + 2)— th cycle the quantum control system must 
be in the state \c 2 ) as required by the quantum program and circuit, one has 
the unitary transformation for the state \nh)\bh)\f r (x)) of the whole quantum 
system: 

U(t 0(i+2) - St r ,t 0{i+1) )\c 2 )\l)\0) = |c 2 )|l)|0) 

where the state of the whole quantum system is |c 2 )|l)|0) at the time t (i+i) 
in the (i + 1)— th cycle and also |c 2 )|l)|0) at the time io(i+2) _ 8t r in the 
(i+2)— th cycle. The second one is that the state of the control state subspace 
is the state \ci) at the instant of time £o(i+i) i n the (i + 1)— th cycle, which 
corresponds to the second situation above. Then for this situation the unitary 
state transformation is given by 

U(t 0(i+2) -ft r ,* o(i+ i))|ci)|l)|0) = |c 2 )|l)|0) 

where the state of the whole quantum system is also |c 2 )|l)|0) at the time 
to(i+2) ~ St r , as required by the quantum program and circuit. One sees that 
these two orthogonal states \ci) and |c 2 ) are changed completely to the same 
state |c 2 ) in the control state subspace by the same unitary transformation 
U(to(i +2 ) — 6to,to(i + i)) during the period from the time to(j+i) to the time 
to(i+2) ~~ Obviously, this is impossible and there is a conflict between the 
unitarity of the propagator U (£o(»+2) — ^o(«+i)) an d the requirement of the 
quantum program and circuit that the quantum control system be in the state 
|c 2 ) at the time to(«+2) — St r in the (i + 2)— th cycle, because the requirement 
leads to non-unitarity of the propagator U(t ( i+2 ) — 5t , i (i+i))- Therefore, 
the quantum control process with the two-state control state subspace S(C) 
could not simulate faithfully and efficiently the reversible and unitary halting 
protocol and could fail to control the quantum computational process in 
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the quantum program Q c if it is constrained to be a single time-dependent 
evolution process. This conflict could be related to the square speedup limit 
of the quantum search algorithm if the quantum program Q c is used to 
construct the quantum search algorithm. Here it should be pointed out 
that this conflict could be avoided in a larger control state subspace rather 
than the two-state control state subspace, but this could lead to that the 
output state of the quantum program is still dependent sensitively upon 
initial states so that the quantum program Q c becomes unvalued for building 
up an efficient quantum search process. Of course, it is usually better to use a 
simpler control state subspace to control the quantum computational process. 

It is well known that a quantum system with a time-independent Hamil- 
tonian satisfies the time-displacement symmetry (or invariance) [5a]. The 
time evolution process of such a quantum system is independent of any ini- 
tial time but depends upon the time difference between the end time and 
the initial time of the process. Therefore, one possible scheme to make the 
output state of the quantum program Q c independent of any initial state 
could be that the Hamiltonian that governs the quantum control process is 
restricted to be time-independent. As shown before, the times {toi} could be 
thought of as the starting times for the state-locking pulse field to really act 
on the control state subspace in the quantum control process. Actually, the 
quantum control process starts to work after the trigger pulse is applied at 
the instant of time — 8t r , and it may be stated simply as that the initial 
halting state \nh) = |0) is changed completely to the state \ci) by the trigger 
pulse Pt in the time interval from the time toi — St r to the time toi, then 
the state \ci) is changed to the state \c?) by the state-locking pulse field in 
the time interval Ato from the time ten to the time toi + Ato, and from the 
time toi + At on, the state |c 2 ) is locked by the state-locking pulse field. 
According to this picture the quantum control process could be expressed 
conveniently in terms of a sequence of unitary transformations on the state 
\ n h)\bh)\fr{x)) of the whole quantum system, 

{P SL ({ ¥ » fc },toi,t O i-^r), J P t }|0)|l)|0) = |ci)|l)|0), (i) 
PsL(W k },t,toi)\ci)\l)\0) 

= {e(t,t Qi )\ Cl ) + e -^*°Vl - |e(Moi)| 2 |c2)}|l)|0), t w < t. (2) 

Here PsL{Wk\-> ^ ten) i s the unitary propagator of the state-locking pulse field 
Psl applying separately to the quantum system during the period from the 
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initial time t 0i to the time t and {Psl toi, toi ~ dt r ), P t } represents the 

unitary propagator when the trigger pulse P t and the state-locking pulse field 
Psl are applied to the quantum system simultaneously during the pulse du- 
ration 5t r of the trigger pulse from the time t 0i — 5t r to the time toi- The 
parameters {y? fc } in the unitary propagator Psi,({</? fe }, t, to) are the control 
parameters of the state-locking pulse field which may be generally dependent 
on the time variable, the spatial variables, or even the quantum field vari- 
ables. Here the unitary propagator {PsL({<^ fc }, toi, toi—St r ), Pt} is really equal 
to the unitary transformation U t during the trigger pulse P t , indicating that 
the unitary transformation U t is really generated by both the state-locking 
pulse field and the trigger pulse. As supposed before, the state-locking pulse 
field is negligible during the trigger pulse in the quantum program Q c , and 
hence the unitary transformation Ut is really generated approximately by the 
single trigger pulse. In the unitary transformation (2) j(t, t 0i ) is the phase 
angle of the amplitude of the state \c2) and e(t, toi) the residual amplitude of 
the state \ci) at the time t after the unitary transformation P SL ({ip k },t,t Qi ) 
acts on the state \ci). As required by the quantum program Q c , the state |ci) 
must be converted completely into the state in the control state subspace 
by the unitary transformation Psx({</? fc }, t, tot) within the period At from 
the time tot to the time t = tot + At - Here the time interval At is shorter 
than the cyclic period AT C minus the duration 5t r of the trigger pulse, that 
is, At < AT C — 5t r . Evidently, when the time t > t 0i + At the absolute 
amplitude value \e(t,toi)\ should approach infinitely zero in theory but can 
not equal exactly zero for every time toi for i = 1, 2, m r even for an ideal 
state-locking pulse field, and generally it should be close to zero for a real 
state-locking pulse field. One may say that the amplitude e(t, toi) equals zero 
in theory if it approaches infinitely zero. The amplitude value \e(t,t 0i ) \ mea- 
sures the performance of a real state-locking pulse field, that is, the smaller 
the amplitude value \e(t,t 0i )\ is for t > t 0i + At , the better the performance 
is for a real state- locking pulse field Psl- If now the quantum control pro- 
cess is governed by a time-independent Hamiltonian during the state-locking 
pulse field, then the time evolution process of Eq. (2) from the state |ci) to 
the state \c2) does not depend directly upon any initial time to* or the end 
time t but it is dependent upon the time difference Ati = t — toi, that is, 
PsLiiVk}-, t, toi) — PsL(Wk}i Ati). This also means that both the amplitude 
e(Moi) of the state \ci) and e~ tl( - t,toi ^ a/1 — \e(t,t 0i )\ 2 of the state |c 2 ) in Eq. 
(2) are dependent upon the time difference Ati, that is, e(t,t 0i ) = e(Ati) 
and 7(t, t 0i ) = j(Ati). It was pointed out before that there are m r differ- 
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ent times {toi, i = 1,2, ...,m r } at most for the quantum control process in 
the quantum program. If the state-locking pulse field PgL is designed in 
such a way that the Hamiltonian that governs the quantum control process 
of Eq. (2) is time-independent, then the time difference Ati = t — t 0i for 
i = 1, 2, m r replaces the time variable t to become a new time variable of 
the quantum control process of Eq. (2), and consequently it is not necessary 
to deal with directly the time variable t in the quantum control process of 
Eq. (2). The unitary transformations (1) and (2) indicate that by the quan- 
tum program Q c with the quantum control process of Eq. (2) governed by 
a time-independent Hamiltonian these m r different initial functional states 
{\f r (x ))} (as well as the initial halting state and branch-control state) can 
be transferred one-by-one to the m r states on the right-hand side of Eq. (2) 
whose phases and amplitudes are dependent upon the time differences {Ati} 
for % = 1,2, ...,m r , respectively. That the Hamiltonian that governs the 
quantum control process of Eq. (2) under the state-locking pulse field is con- 
strained to be time-independent could be helpful for designing a state-locking 
pulse field with a good performance to control the quantum computational 
process in the quantum program Q c . 

However, the time-displacement symmetry is not sufficient to solve thor- 
oughly the conflict mentioned above and it can not yet figure out completely 
a state-locking pulse field. This is because although the amplitude e(Ati) 
of the state \ci) in Eq. (2) is independent of any single initial time t i, h 
is still dependent upon the time difference Ati. Obviously, a different time 
difference Ati generally leads to a different residual amplitude e(Ati) of the 
state |ci) in Eq. (2). Then there could be still a problem that the amplitude 
\e(Ati)\ could fail to be close to zero during the trigger pulse P t from the 
time tok — 5t r to time tok, here the time difference (io& — St r — t 0i ) > At for 
k = i + i + 2, m r . As a result, the residual state \ci) with a large ampli- 
tude |e(Atj)| may be changed back to the initial halting state \rih) = |0) by 
the trigger pulse P t again. Obviously, this result is not consistent with the 
requirement of quantum program that any residual amplitude |e(Atj)| of the 
state \ci) in Eq. (2) equal zero in theory when Ati > At for % = 1,2, ...,m r . 
Now suppose that these two states \ci) and \c2) of the control state subspace 
are degenerative eigenstates of the time-independent Hamiltonian H and the 
unitary propagator Psx({<^ fc }, t, toi) = exp(-iHAti) (H = 1). Then these two 
eigenstates have the same eigenvalue and the eigen-equations are written as 
H\ci) = A|ci) and H\c 2 ) = A|c 2 ) with the common eigenvalue A, respectively. 
Thus, the time evolution process similar to Eq. (2) with the initial superpo- 
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sition state a\c\) + 6|c 2 ) of the control state subspace S(C), which is drived 
by the time-independent Hamiltonian H, may be generally expressed by 

PsdWk}, A*i)(a|ci) + 6|c2»|l)|0) = exp(-iAA*0(a|ci) + &|c 2 ))|l)|0). 

Evidently, only the global phase factor of the state a\ci) + b\c 2 ) is changed by 
the unitary propagator -Psz,({<^ fc }, Atj). The global phase factor is dependent 
upon the time difference At i? but the absolute amplitude of the state does not 
change as the time difference. Therefore, the time-independent Hamiltonian 
that drives the time evolution process of Eq. (2) in the two-state control state 
subspace is very useful for keeping the amplitude of any state of the control 
state subspace unchanged for a long time. However, such a Hamiltonian 
could not be suitable for transferring the state \ci) to the state \oi) in the 
control state subspace, while the state transfer is necessary for a quantum 
control process such as Eq. (2) to simulate the reversible and unitary halting 
protocol. 

It is shown above that it is not sufficient to build up the state-locking 
pulse field with a good performance by constraining the Hamiltonian of the 
quantum control system with the two-state subspace S(C) under the state- 
locking pulse field to be time-independent. A more suitable Hamiltonian 
that governs the quantum control process of Eq. (2) may be dependent upon 
the spatial variables and/or the quantum field variables but independent of 
the time variable so that the propagator P SL ({ip k }, At) is still dependent 
upon the time difference At. In a quantum computer architecture differ- 
ent quantum bits of the quantum system of the quantum computer must be 
addressed spatially or distinguished from each other by some properties of 
the quantum system such as the spectroscopic properties so that they can 
be manipulated at will. However, such spatial-dependent properties of the 
quantum system are static and different from a space-dependent evolution 
process. A quantum computational process generally is considered as a uni- 
tary time evolution process of a quantum system [3] which may be generally 
dependent upon both time and space. According as quantum mechanics [5], 
time- and space-dependent evolution processes of a quantum system such as 
the conventional quantum scattering process, the quantum tunneling process, 
the quantum collision process, the molecular chemical dissociation process, 
and so on, obey the Schrodinger equation as well and hence they are also 
governed by the unitary quantum dynamics in time and space. The force 
to drive a time- and space-dependent evolution process such as a quantum 
scattering process usually could be the motional momentum of a particle 
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or an electromagnetic field and so on. The time- and space-dependent uni- 
tary evolution processes could also be used to build up quantum compu- 
tational processes just like the conventional quantum gates [25], although 
space-dependent unitary evolution processes usually are more complicated 
and difficult to be manipulated at will than those space-independent ones. 
A large advantage of a time- and space-dependent unitary evolution process 
over a space-independent one for building up a quantum computational pro- 
cess could be that a time- and space-dependent unitary evolution process may 
be manipulated separately either in the time dimension or in the space di- 
mensions or in both the time and space dimensions. While a time-dependent 
and space-independent unitary evolution process could be inadequate as the 
quantum control process of Eq. (2) for the two-state control system, a time- 
and space-dependent unitary evolution process could be better to act as the 
quantum control process. Thus, a time- and space-dependent unitary evo- 
lution process could be very useful for some specific purposes in quantum 
computation, although a quantum computational process usually is simply 
designed to be a space-independent unitary evolution process of a quantum 
system and any space-dependent evolution processes of the quantum system 
are suppressed so that in algorithm the quantum computational process may 
be constructed as simply as possible. As shown before, if the quantum con- 
trol process is purely time- dependent in the two-state subspace S(C), there 
is a conflict between the unitarity of the quantum control process and the 
performance of the quantum control process that the state \oi) in the control 
state subspace S{C) is kept unchanged for a long time by the state-locking 
pulse field. Then it could be better to use a time- and space-dependent uni- 
tary evolution process such as a quantum scattering process to realize the 
quantum control process of Eq. (2), meanwhile the quantum computational 
process may be set to a single time-dependent unitary evolution process of 
the quantum system of the quantum computer. This is because both the 
single time-dependent quantum computational process and the time- and 
space-dependent quantum control process may be manipulated separately in 
time and space and hence they become almost independent upon each other. 
If the Hamiltonian to drive the time- and space-dependent quantum con- 
trol process of Eq. (2) is space-dependent and time-independent, then the 
energy of the quantum control system is conservative during the quantum 
control process of Eq. (2) and hence both these two states \c±) and |c 2 ) of 
the control state subspace have the same energy. Then the quantum control 
process of Eq. (2) in the quantum control system is directly dependent upon 
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the time interval Atj = t — t 0i rather than the initial time £ i or the time 
variable t separately. Since now these two states \ci) and \c2) are degenerate 
in energy only their global phase is dependent upon the time difference Ati 
but their amplitudes are not during the quantum control process of Eq. (2), 
and hence the state \c2) may be kept unchanged for a long time under the 
state-locking pulse field. However, here the state transfer from the state \c\) 
to the state \c2) in the control state subspace S(C) could be achieved by the 
time- and space-dependent unitary evolution process such as the quantum 
scattering process if the time-independent and space-dependent Hamiltonian 
of the quantum control system is chosen suitably. Obviously, the quantum 
scattering process should be designed suitably according to the properties of 
an ideal state-locking pulse field. 

3. An atomic physical model to simulate efficiently the reversible 
and unitary halting protocol 

The detailed analysis in the former sections for the quantum control pro- 
cess shows that a quantum control process that simulates faithfully and ef- 
ficiently the reversible and unitary halting protocol should contain a time- 
and space-dependent unitary evolution process such as a quantum scattering 
process and the control state subspace S(C) should not be restricted to be 
only the smallest two-state subspace. Generally, the quantum control process 
to simulate efficiently the reversible and unitary halting protocol in the quan- 
tum program Q c may be implemented in a real quantum physical system. 
A trapped atomic-ion system has been proposed as a real physical system 
to implement quantum computation [33]. Here a simple quantum physical 
system of an atomic ion or a neutral atom in one-dimensional double-well 
potential field is proposed as the quantum control system of the quantum 
program and circuit Q c with the control state subspace S(C) which is larger 
than the two-state control state subspace. Now the simple atomic physical 
system will be used to illustrate how the quantum control process simu- 
lates really the reversible and unitary halting protocol and how to construct 
explicitly the state-locking pulse field. In this simple physical system the 
atomic ion or the neutral atom in the double-well potential field is called 
the halting-quantum-bit atom or the halting-qubit atom briefly. Hereafter 
the halting-qubit atom is referred to as the atomic ion or the neutral atom 
in the double-well potential field unless otherwise stated. In the double-well 
potential field the left-hand potential well could be approximately a con- 
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ventional harmonic potential well, while the right-hand potential well could 
be simply a square potential well. Here also suppose that the right-hand 
square potential well is sufficiently wide such that the halting-qubit atom 
motions almost freely in one- dimensional space in the square potential well. 
The intermediate part between these two potential wells is a square potential 
barrier which is used to block free transportation of the halting-qubit atom 
from one potential well to another in the double-well potential field. Here 
the maximum height of the right potential wall of the left-hand potential 
well is just equal to the height of the intermediate square potential barrier. 
The left and right potential walls of the double-well potential field may be in- 
finitely high in theory, but the intermediate square potential barrier is finitely 
high and wide. Actually the left-hand potential well should be an asymmet- 
ric harmonic potential well with an infinitely high left potential wall and 
a finitely high right potential wall, respectively. If the intermediate square 
potential barrier is infinitely high and finitely wide then any two states of 
the halting-qubit atom in the left- and right-hand potential wells respectively 
are completely orthogonal to one another, and if the square potential barrier 
is high enough and finitely wide then the two states are also considered to 
be orthogonal to one another approximately. The time-independent double- 
well potential field in the atomic physical system could be generated by the 
external electromagnetic field [30]. More generally the double- well potential 
field could be thought of as an effective potential field of the halting-qubit 
atom, that is, this potential field could be generated effectively by the in- 
teraction between the halting-qubit atom and the external electromagnetic 
field, the interactions between the halting-qubit atom and those atoms of 
the computational state subspace in the quantum system of the quantum 
computer, and those interactions between the halting-qubit atom and its en- 
vironment. The halting-qubit atom in the double- well potential field could 
be coupled to those quantum-bit atoms of the computational state subspace 
either through the Coulomb interactions between charged atomic ions [31] 
or through the atomic dipole-dipole interactions between atomic ions in an 
atomic- ion physical system [32b], while the halting-qubit atom could also be 
coupled to other atoms by the dipole-dipole interactions of neutral atoms 
in a neutral atomic physical system [32a, 32c]. These interactions may be 
used to set up two-qubit quantum gates between the halting-qubit atom and 
those quantum-bit atoms of the computational state subspace in the quan- 
tum system. With these two-qubit quantum gates and one-qubit quantum 
gates one can build up those efficient unitary operations that act on both the 
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halting qubit and those qubits of the computational state subspace such as 
the unitary operation in the quantum program Q c and V A C in the quantum 
control unit Qh (see below). It is required by the quantum control process 
that these interactions be available only when the halting-qubit atom is in the 
left-hand harmonic potential well, while they are negligible when the halting- 
qubit atom is in the right-hand potential well due to that the halting-qubit 
atom in the right-hand potential well is much farther from those atoms of 
the computational state subspace, and both the Coulomb interaction and 
the dipole-dipole interaction may become very weak as the distance between 
the interacting atoms become large [31, 32]. Therefore, when the halting- 
qubit atom enters the right-hand potential well from the left-hand one these 
interactions between the halting-qubit atom and those atoms of the compu- 
tational state subspace should be decoupled and can be negligible so that the 
two-qubit quantum gates between the halting-qubit atom in the right-hand 
potential well and those atoms of the computational state subspace can not 
be built up effectively. More generally any quantum operations involved in 
the halting-qubit atom in the left-hand potential well could be really hung up 
when the halting-qubit atom leaves the left-hand potential well. Therefore, 
the halting operation could be achieved due to that these quantum opera- 
tions are hung up when the halting-qubit atom enters into the right-hand 
potential well from the left-hand one. 

Generally, an atom has both the internal electronic states and the center- 
of-mass motional states. Here a center-of-mass motional state of an atom 
may be a wave-packet motional state, and for a heavy particle the wave- 
packet picture in quantum mechanics is close to the classical particle picture 
[5a] . The internal electronic states of an atom are generally quantized bound 
states, but the center-of-mass motional states may be either the quantized 
bound states in a potential field or the continuous states in free space [5a, 5c] . 
A quantum bit of an atom generally may be chosen suitably as a pair of the 
internal electronic ground states of the atom, but sometime the quantized 
motional states of an atom in a harmonic potential field are also taken as 
the quantum bits in quantum computation [30, 33]. Thus, the halting-qubit 
atom has the internal electronic states and also the center-of-mass motional 
states in the double-well potential field. The halting quantum bit generally 
may be chosen as a pair of the specific internal electronic ground states 
of the halting-qubit atom. The time-independent double-well potential field 
generally affects the center-of-mass motional states of the halting-qubit atom 
[30], but it could have a negligible effect on the internal electronic states of 
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the halting-qubit atom so that these internal electronic states could keep 
unchanged when the halting-qubit atom moves from one potential well to 
another [34] . Actually, the internal electronic states of the halting-qubit atom 
are determined mainly by the internal interactions of the halting-qubit atom 
itself, although the complete quantum structure of the halting-qubit atom 
in the double-well potential field is determined by both the external double- 
well potential field and the internal interactions of the halting-qubit atom 
itself. As shown in the previous section, the implementation of the reversible 
and unitary halting protocol in the atomic physical system is involved in the 
time- and space-dependent unitary evolution process that the halting-qubit 
atom moves from one potential well to another in the double-well potential 
field. Then the center-of-mass motional states and especially the wave-packet 
motional states of the halting-qubit atom in the double-well potential field 
will play an important role in implementing the reversible and unitary halting 
protocol in the atomic physical system. 

If now the atomic physical system of the halting-qubit atom in the one- 
dimensional double-well potential field is considered to be a quantum control 
system, then one must define explicitly the initial halting state \rih) and 
the states \ci) and \c2) of the control state subspace S{C) in the quantum 
program Q c . First of all, the total wavefunction of the halting-qubit atom in 
the one-dimensional double-well potential field may be generally written as 

\n' h ,CM,R) = \n' h )\CM,R). 

Here the states \n' h ) and \CM,R) represent the internal electronic state and 
the center-of-mass motional state of the halting-qubit atom, respectively. 
The integer n' h and the index CM are the quantum numbers of the internal 
state \n' h ) and the motional state \CM, R), respectively, and R is the spatial 
coordinate of the center of mass of the halting-qubit atom with the motional 
state \CM,R) in the double-well potential field. The center-of-mass spatial 
position R is generally time-dependent, i.e., R = R(t). Actually, CM is also 
used to represent the expectation value or eigenvalue of the motional energy 
(or momentum) of the halting-qubit atom in the double-well potential field 
particularly when the halting-qubit atom is in a wave-packet motional state. 
Before the quantum program Q c is performed the halting-qubit atom is lo- 
cated in the left-hand harmonic potential well of the double-well potential 
field by the ionic or atomic trapping techniques [30, 35]. For convenience 
the left-hand potential well may be prepared temporarily as a conventional 
harmonic potential well before the quantum program starts to work. This 
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can be achieved easily by setting the height of the right-hand wall of the 
left-hand potential well to be sufficiently large, since the left-hand potential 
well can be thought of approximately as a conventional harmonic poten- 
tial well when the right-hand potential wall is sufficiently high (note that 
the left-hand potential wall is infinitely high in theory). Thus, before the 
quantum program starts the internal and motional states of the halting- 
qubit atom may be really prepared to be the ground internal state and the 
ground motional state of the conventional harmonic oscillator by the laser 
cooling techniques [30, 36], respectively. Now the global ground state of the 
halting-qubit atom in the left-hand harmonic potential well may be written 
as \0,CM0,Rq) = \0)\CM0,Rq), which is the product state of the ground 
internal state \n' h ) = |0) and the ground motional state \CM0, i? ) of the 
atom in the harmonic potential well. Note that the ground motional state 
\CM0, Rq) of the atom in the conventional one- dimensional harmonic poten- 
tial well is a Gaussian wave-packet motional state [5a] . After the total ground 
state \0,CM0,Rq) is prepared the left-hand harmonic potential well is sud- 
denly changed back to the original double-well potential field at the starting 
time of the quantum program. Actually, this process changes merely the 
sufficiently high right-hand wall of the harmonic potential well to the finitely 
high one of the left-hand potential well of the double-well potential field. Ac- 
cording as quantum mechanics [5a] that a wavefunction of a quantum system 
must be continuous in time, the state \0,CM0, R Q ) still remains unchanged 
at the starting time of the quantum program when the harmonic potential 
field is suddenly changed back to the double-well potential field. Then at 
the starting time of the quantum program the motional state for the halting- 
qubit atom in the left-hand potential well of the double-well potential field 
is just \CM0, Rq) and hence still a one-dimensional Gaussian wave-packet 
motional state. The wave-packet state |0, CMO, Rq) could not be exactly 
the global ground state of the halting-qubit atom in the double-well poten- 
tial field. However, the energy of the wave-packet state |0, CMO, Ro) may 
be very close to that one of the global ground state of the halting-qubit 
atom in the double-well potential field if the double-well potential field is 
designed suitably. Therefore, the initial halting state \rih) = |0) in the quan- 
tum program Q c may be set to the wave-packet state |0,CM0, Ro) . Then 
the halting quantum bit may be chosen as the two ground internal states 
{\n' h ), n' h = 0,1} of the wave-packet states {\n' h , CMO, R )}. The state |ci) 
of the control state subspace in the quantum program Q c could be taken as 
the wave-packet state |1, CMl(i i)> Ri(toi)) with the internal state \n' h ) = |1) 
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and the wave-packet motional state \CMl(t 0i ), Ri(t 0i )) of the halting-qubit 
atom in the left-hand potential well. Here the internal state \n' h ) = |1) 
could be chosen as another hyperfine ground internal electronic state of the 
halting-qubit atom different from the hyperfine ground internal electronic 
state \n' h ) = |0), and the wave-packet spatial position Ri(t 0i ) is within the 
left-hand potential well. The wave-packet motional state |CMl(t j)> Ri(toi)) 
is generated from the wave-packet motional state \CM0, Ro) by the trigger 
pulse P t (see below). The mean motional energy (CM1) of the wave-packet 
motional state \CMl(t 0i ), Ri(t 0i )) is much higher than that one (CMO) of 
the motional state \CM0, R ) and also the height of the intermediate poten- 
tial barrier in the double-well potential field. Thus, the wave-packet state 
|ci) = |1, CMl(t 0i ), Ri(t 0i )) is an unstable state. When the halting-qubit 
atom in the left-hand potential well is in the unstable wave-packet state 
|ci), its motional energy (CM1) is much higher than the height of the in- 
termediate potential barrier so that the halting-qubit atom is driven by the 
high motional momentum (CM1) of the atom to pass the intermediate po- 
tential barrier to enter into the right-hand potential well. This is a time- 
and space-dependent quantum scattering process for the halting-qubit atom 
in the double-well potential field. This quantum scattering process will be 
taken as the quantum control process of Eq. (2). Note that the quantum 
scattering process starts at the initial time t 0i . Now suppose that the halting- 
qubit atom enters completely into the right-hand potential well at the time 
t m i = toi + At in the quantum scattering process and at the time t m i the 
wave-packet state of the halting-qubit atom in the right-hand potential well 
is denoted as |1, CM2(t m i), R2{t m i)). Here the wave-packet spatial position 
R2(t m i) is within the right-hand potential well. Since the quantum scatter- 
ing process for the halting-qubit atom from one potential well to another 
in the double-well potential field does not change the ground internal states 
of the atom [34] both the wave-packet states |ci) = |1, CMl(t 0i ), Ri(t 0i )) 
and |c' 2 ) = \l,CM2(t mi ), R2(t mi )) have the same internal state with n' h = 1. 
Now another state |c 2 ) of the control state subspace in the quantum pro- 
gram Q c could be temporarily set to the wave-packet state |c' 2 ). Actually, 
the state |c2) in the quantum program Q c will correspond to any wave-packet 
state of the halting-qubit atom in the right-hand potential well, as can be 
seen later, for the control state subspace of the atomic physical system is 
not a two-state subspace. Due to the motional energy conservation during 
the quantum scattering process the motional energy (CM1) of the wave- 
packet state |1, CMl(t 0i ), Ri(t 0i )) of the halting-qubit atom in the left-hand 
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potential well is really equal to that one (CM2) of the wave-packet state 
\l,CM2(t m i), Ri(t m i)) of the atom in the right-hand potential well. Thus, 
the wave-packet state |1, CM2(t m i), Rzitmij) is also an unstable state just 
like the wave-packet state |1, CMl(t j), Ri(toi))- The spatial spreads of these 
two wave-packet motional states \CMl(t 0i ) : Ri(t 0i )) and \CM2(t mi ), i^^mi)) 
could be generally different from each other, but these two motional states 
could be approximately Gaussian in coordinate space [5a, 5c] just like the mo- 
tional state \CM0, Rq). Though these two wave-packet motional states could 
not be exactly orthogonal to one another, they could be almost completely or- 
thogonal to one another if the wave-packet spatial distance \R2(t m i)— Ri(toi)\ 
between their spatial positions Ri(t 0i ) and -R2 (^mi) is large enough, because 
the overlapping integral between these two wave-packet motional states de- 
cays exponentially as the wave-packet spatial distance increases. Therefore, 
the width of the intermediate potential barrier of the double-well potential 
field must be large enough such that any wave-packet motional state of the 
halting-qubit atom within one potential well is almost completely orthogonal 
to that one within another potential well in the double-well potential field. 

In the atomic physical system the time-independent double-well poten- 
tial field could be considered as one of the components of the state-locking 
pulse field Psl- More generally, a complete state-locking pulse field in the 
atomic physical system could consist of three parts below. The first part is 
the double- well potential field itself. The second is the sequences of the time- 
and space-dependent electromagnetic pulse fields which are applied only to 
the right-hand potential well. As can be seen below, these sequences in- 
clude mainly the unitary decelerating sequence and the unitary accelerating 
sequence. The unitary decelerating sequence (the unitary accelerating se- 
quence) is used to decelerate (accelerate) the halting-qubit atom in motion 
in the right-hand potential well. In this part the sequences of the electromag- 
netic pulse fields are used to manipulate coherently the wave-packet motional 
state of the halting-qubit atom in the right-hand potential well so that the 
halting-qubit atom can stay in the right-hand potential well for a long time 
required by the quantum computational process in the quantum program. 
The third part is the sequences of the time- and space-dependent electro- 
magnetic pulse fields and includes also the time-dependent potential fields 
which are applied mainly to the left-hand potential well after the quantum 
program Q c terminates. In this part the electromagnetic pulse fields asso- 
ciated with the unitary accelerating sequence in the second part are used 
to drive the halting-qubit atom in the right-hand potential well to go back 
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to the left-hand potential well after the quantum computational process fin- 
ished, and at the same time transfer each of possible wave-packet states of 
the halting-qubit atom in the right-hand potential well to some given wave- 
packet state of the atom in the left-hand potential well so that the final state 
of the halting-qubit atom in the left-hand potential well is not dependent 
sensitively upon these possible wave-packet states. Therefore, in both the 
second and third parts these electromagnetic pulse fields should have a negli- 
gible effect on any motional state of the halting-qubit atom if the atom locates 
within the left-hand potential well during the quantum computational pro- 
cess. In this and next sections it is discussed how the double-well potential 
field affects the wave-packet state of the halting-qubit atom before the atom 
leaves the left-hand potential well in the quantum computational process, 
while in the section 5 it will be discussed how the sequences of the time- and 
space-dependent electromagnetic pulse fields acts on the halting-qubit atom 
to keep the atom in the right-hand potential well for a long time and how by 
the electromagnetic pulse fields and the time-dependent potential fields each 
of the possible wave-packet states of the halting-qubit atom in the right-hand 
potential well is changed to some given wave-packet state of the atom in the 
left-hand potential well after the quantum computational process finished. 
The time-independent double-well potential field must be designed properly. 
Before the halting-qubit atom leaves the left-hand potential well during the 
quantum computational process the atom may have its zero-point oscillatory 
motion in the left-hand harmonic potential well, and it could also penetrate 
through the intermediate potential barrier and enter into the right-hand po- 
tential well due to the quantum tunneling effect [5a] even if the motional 
energy of the halting-qubit atom in the left-hand potential well is lower than 
the height of the intermediate potential barrier. The zero-point oscillatory 
motion is allowed normally in quantum computation [30], but the quantum 
control process that simulates efficiently the reversible and unitary halting 
protocol could become degraded if the halting-qubit atom enters into the 
right-hand potential well in a non-negligible probability due to the quan- 
tum tunneling effect before the halting operation is performed according as 
the quantum program. Thus, the double-well potential field should be de- 
signed in such a way that the probability for the halting-qubit atom going 
from the left-hand potential well to the right-hand one due to the quantum 
tunneling effect should be minimized and can be neglected for the quantum 
computational process. The height and width of the intermediate potential 
barrier in the double-well potential field may control the quantum tunneling 
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effect, that is, the higher and wider the intermediate potential barrier, the 
lower the penetrating probability for the halting-qubit atom. For example, 
consider a simple physical model that a free particle with motional energy 
Eh and mass rrih hits a square potential barrier with height Vo and width a 
[5a]. The particle will be reflected and/or transmitted by the square poten- 
tial barrier. If the motional energy Eh of the particle is much less than the 
potential barrier height Vq such that (3a » 1 with (3H = \/2m/^ — Eh), 
then the transmission coefficient of the particle is approximately proportional 
to the exponential factor exp(— 2 (3 a) [5a]. Therefore, the probability for the 
free particle to penetrate through the potential barrier is also proportional 
to the exponential factor exp(— 2(3a). When the potential barrier height Vq 
and width a are large enough, this probability falls off rapidly. Now a bound 
particle with the same motional energy Eh « Vq like the halting-qubit atom 
in the left-hand potential well is more difficult to penetrate through the same 
square potential barrier than a free particle. An atom has a much heavier 
mass rrih than an electron. The probability for the halting-qubit atom with 
the motional energy Eh « Vq in the left-hand potential well to penetrate 
through the intermediate square potential barrier with the height Vq and 
width a decays rapidly in an exponential form as the height Vq, the width 
a, and the atomic mass rrih increase. Evidently, if the intermediate potential 
barrier is high and wide enough such that the ground-state motional energy 
(CMO) of the halting-qubit atom is much lower than the barrier height, then 
the quantum tunneling effect could have a negligible effect on the initial 
halting state \rih) = |0, CMO, Ro) and also the state \rih) = \1,CM0, Rq) of 
the halting quantum bit, and hence the wave-packet states {\n' h ,CM0, R )} 
with n' h = 0, 1 of the halting qubit atom in the left-hand potential well may 
keep almost unchanged during the quantum computational process before 
the halting-qubit atom leaves the left-hand potential well due to the halting 
operation in the quantum program. 

The coherent stimulated Raman adiabatic passage (STIRAP) method 
has been used to prepare the 'Schrodinger Cat' superposition state of a 
trapped atom [37]. Here the coherent STIRAP method [37-41] may also 
be used to transfer selectively the halting state \rih) = \1,CM0, R ) of 
the halting-qubit atom in the left-hand potential well to the unstable state 
\ci) = \l,CMl(toi), Ri(t 0i )) of the control state subspace. Therefore, the 
state-dependent trigger pulse P t in the quantum program Q c could be cho- 
sen as the coherent Raman adiabatic laser pulse. The spatial action zone of 
the coherent Raman adiabatic laser pulse P t is confined within the left-hand 
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potential well. A coherent Raman adiabatic laser pulse P t consists of a pair 
of the coherent adiabatic laser beams. Here denote these two adiabatic laser 
beams as A and B, respectively. When the coherent Raman adiabatic laser 
pulse P t is applied to the halting-qubit atom in the left-hand potential well, 
one of these two adiabatic laser beams (e.g., the beam A) first excites se- 
lectively the transition of the halting-qubit atom from the wave-packet state 
\n' h , CMO, Rq) with the internal state \n' h ) = |1) to some specific excited 
state \n e ,CM e ,R e ), meanwhile the halting-qubit atom in the excited state 
\n e , CM e , R e ) is stimulated by another adiabatic laser beam B to jump to 
the state \c±) = |1, CMl(^), Ri(t 0i )). Here the excited internal state \n e ) 
(n e 7^ 0, 1) of the excited state \n e , CM e , R e ) has a higher internal energy 
than the ground internal states {\n' h ), n' h = 0, 1} and the halting-qubit atom 
with the wave-packet motional state \CM e ,R e ) is still within the left-hand 
potential well. For example, if the halting-qubit atom is chosen as a single 
9 B+ ion, then these two internal states {\n' h ), n' h = 0,1} could be taken 
as the ionic hyperfine ground states 2 Si/ 2 (F = 2,m F = —2) and 2 S 1 / 2 
(F = I, trip = —1), respectively, while the excited internal state \n e ) could 
be the excited electronic state 2 P\/2 {F = 2, rnp = —2) of the ion [37a]. This 
state-dependent excitation process under the coherent Raman adiabatic laser 
trigger pulse P t may be simply expressed as 

\1,CM0,R ) A \n e ,CM e ,R e ) £ \l,CMl(t 0i ),R 1 {t 0i )). 

In order that only the two desired transitions: |1, CMO, R ) <-> \n e , CM e , R e ) 
and \n e , CM e , R e ) <-> |1, CMl(i i)> Ri(toi)) are excited selectively by the co- 
herent Raman adiabatic laser pulse Pt the frequencies of both the laser beams 
A and B must be close to the resonance frequencies of these two desired tran- 
sitions, respectively, and they are also much far from resonance frequencies of 
any other transitions including the transition: |0,CM0,i? ) <-> \n e ,CM e: R e ). 
Any wave-packet state \n' h) CM, R) of the halting-qubit atom with the inter- 
nal state \n' h ) ^ |1) will not be affected effectively by any one of these two 
adiabatic laser beams. Thus, the coherent Raman adiabatic laser pulse Pt 
does not act on any wave-packet state with the internal state \n' h ) ^ |1) 
such as the initial halting state |0, CMO, R ). On the other hand, in order 
to suppress irreversible spontaneous emission processes both the laser beams 
A and B are detuned properly from the excitation state \n e ,CM e ,R e ). It 
might be better that the wave vector difference of these two laser beams 
A and B is set to point to the left-hand potential wall of the double-well 



47 



potential field as the left-hand potential wall may be sufficiently high in 
practice. This means that the halting-qubit atom moves along — x axis to- 
ward the left-hand potential wall under the action of the coherent Raman 
adiabatic laser pulse P t . Obviously in the state-dependent excitation pro- 
cess the halting state |1,CM0, R ) is excited to the higher motional energy 
state |ci) = \l,CMl(t 0i ), Ri(t 0i )) by the coherent Raman adiabatic laser 
pulse P t . In particular, the mean motional energy of the motional state 
\CMl(toi), Ri(toi)) of the halting-qubit atom must be much higher than the 
height of the intermediate potential barrier so that the unitary quantum scat- 
tering process can take place automatically for the halting-qubit atom from 
the left-hand potential well to the right-hand one. 

The quantum control process in the quantum program Q c usually should 
be modified properly when it is implemented in a real quantum physical sys- 
tem such as the simple atomic physical system mentioned above. For the 
quantum control system of the atomic physical system the initial halting 
state \rih) and these two states |d) and |c 2 ) of the control state subspace 
S(C) in the quantum program Q c are defined as the wave-packet states of 
the halting-qubit atom in the left-hand potential well: \nh) = |0, CMO, Rq), 
|ci) = \l,CMl(t 0i ),Ri(toi)), and |c 2 ) = |c' 2 ) = |1, CM2(t mi ), R 2 (t mi )), re- 
spectively. The state-dependent trigger pulse P t is taken as a suitable co- 
herent Raman adiabatic laser pulse that is applied to within the left-hand 
potential well in space. The state-locking pulse field in the atomic physi- 
cal system consists of the double-well potential field and the sequences of 
the time- and space-dependent electromagnetic pulse fields and the time- 
dependent potential fields, as mentioned in the previous paragraphs. Now a 
quantum control process (or unit) Qh of the atomic physical system which 
replaces the statement P7 of the quantum program Q c is designed to simulate 
efficiently the reversible and unitary halting protocol. The quantum control 
process (or unit) Qh is written as 

While \f r (x)) = \f r (x f )), 

Do U c h : \n h )\f r (x f )) = \0,CM0, R )\l) - \0,CM0, Ro)\0) , 
V h c : \n h )\f r (x)) = \0,CM0,Ro)\0) - |1, CMO, Rq)\0) 
State-dependent excitation process (P t ) : 

K> = \1,CM0,R ) -> |d) = |l,CMl(^),#i(^)> 
Quantum scattering process in time and space (Psl) '■ 
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|d) = |l,CMl(^),#i(*o*)) -> |c' 2 ) = \l,CM2(t mi ),R 2 (t mt )) 

Here the desired functional state is \f r (%f)) — |1) if the quantum pro- 
gram Q c works in the multiplicative-cyclic-group state subspace S(m r ) and 
\f r {xf)) = |0) if the quantum program Q c works in the additive-cyclic-group 
state subspace S(Z mr ) and in later case the unitary operation U£ may be 
omitted from the quantum control unit Qh- In the quantum control unit 
Qh the unitary operation U% is the same as the original one in the quantum 
program Q c : 

U c h : \0,CM0,R )\b h )\l) <- |0,C7MO,i2o)|6 ft )|0), 

and the new conditional unitary operation Vfi is defined as 

V£ : |0,CMO,i2o)|6 ft )|0) <- |1, CMO, R )\b h )\0) 

with bh = 0,1, and the conditional unitary operation U^ r during the state- 
dependent trigger pulse P t is simply defined by 

U c tr : \l,CM0,Ro)\b h )\f r (x)) <- \l,CMl(t i),Ri(t 0t ))\b h )\f r (x)). 

Here the unitary operation U£ r corresponds to the original one U t in the quan- 
tum program Q c . The conditional unitary operations U% and V h c generally 
are dependent upon both the functional states and the internal states of the 
halting-qubit atom but may be independent of any motional states of the 
atom. They may be built up efficiently out of the interactions between the 
halting-qubit atom in the left-hand potential well and those atoms of the com- 
putational state subspace. It follows from the quantum control unit Qh that 
only when the functional state \f r (x)) becomes the desired state \f r (xf)) can 
the unitary operations and V£ take a real action on the quantum system of 
the quantum computer. One of the important processes in the quantum con- 
trol unit Qh is the state-dependent excitation process involved in the trigger 
pulse P t with the pulse duration 8t r . Though the coherent Raman adiabatic 
laser trigger pulse P t is applied to the halting-qubit atom in the left-hand 
potential well at the starting time toi — 5t r for every cycle of the quantum pro- 
gram Q c with the cyclic index i — 1, 2, m r , it can only take a real action on 
those wave-packet states of the atom with the internal state \n' h ) = |1) such as 
the halting state |1, CMO, Rq). Therefore, only when the initial halting state 
\n h ) = |0,CM0,i? ) is changed to the halting state \n h ) = \1,CM0, R ) by 
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the conditional unitary operation V£ in the quantum program can the trigger 
pulse P t excite the state \nh) = \1,CM0,R ) with the ground motional en- 
ergy (CMO) to the unstable wave-packet state \c\) = |1, CMl(£ i), Ri(toi)) 
with a much higher motional energy (CM1) and a different spatial position 
Ri(toi) 7^ Ro in the left-hand potential well. Since the mean motional energy 
(CM1) of the wave-packet state |1, CMl(t 0i ), Ri(t 0i )) is much higher than 
the height of the intermediate potential barrier the halting-qubit atom in the 
left-hand potential well with the wave-packet state |1, CMl(toi), Ri{toi)) can 
easily pass the intermediate potential barrier to enter into the right-hand po- 
tential well. This process is a quantum scattering process in space and time: 
|1, CMl(t 0i ), Ri(t 0i )} — > \l,CM2(t mi ), R 2 (t m i)), that is, a time- and space- 
dependent unitary evolution process which is driven by the motional momen- 
tum of the halting-qubit atom in the unstable state |1, CMl(toi), -Ri(^oi))- 
This is a key process for the quantum control process Qh to achieve the re- 
versible and unitary halting operation. The halting operation will take place 
when the halting-qubit atom leaves the left-hand potential well because the 
coherent Raman adiabatic laser trigger pulse P t does not have a real action 
on any wave-packet state of the halting-qubit atom if the atom is not in the 
left-hand potential well. On the other hand, if the halting-qubit atom enters 
into the right-hand potential well, then any two-qubit quantum gates become 
unavailable between the halting-qubit atom and those atoms of the compu- 
tational state subspace due to that the effective interactions vanish between 
the halting-qubit atom and those atoms. Then any one of the unitary oper- 
ations and Vfi in the quantum control unit Qh becomes yet unavailable 
and consequently the halting operation is achieved too. 

4. The unitary evolution processes for the quantum program 
and circuit in the atomic physical system 

Now the original quantum control unit (the statement P7) of the quan- 
tum program Q c is replaced with the quantum control unit Qh of the atomic 
physical system. The time evolution process of the atomic physical system of 
the quantum computer under the quantum program and circuit Q c is inves- 
tigated in detail below. This investigation will be helpful for understanding 
more clearly and deeply the general properties of a state-locking pulse field 
and especially the properties of the unitary transformations related to the 
state-locking pulse field. First of all, the quantum program and circuit Q c 
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with the quantum control unit Qh may be written in the simple form 
Q c = {P SL : OFF}{U}P SL (Wih At )U^U c h U^{P SL : ON}. 

Here the initial state of the quantum circuit Q c is set to the basis state 
|n h )|6 h )|/ r (a:)) = |0, CMO, R ) |0) |/ r (x )>. The state-locking pulse field P SL 
is applied continuously to the quantum system from the beginning to the 
end of the quantum circuit. During the periods of these unitary operations 
Uj, U£ r , Uf V£, and U£ the quantum system is really acted upon simulta- 
neously by both the state-locking pulse field and these unitary operations. 
The unitary transformation PslH^i}, A^o) represents the quantum scatter- 
ing process during the period from the time to« to the time t m i = toi + Ato- 
Obviously, if a unitary operation is applied only to the computational state 
subspace it commutes with the unitary propagator PsL{{Vk\-> ^ ^o) of the 
state-locking pulse field as the state-locking pulse field is applied only to 
the halting-qubit atom in the double-well potential field. Now the unitary 
functional operation and the unitary operation U§ are applied only to 
the computational state subspace in the quantum circuit. Thus, the unitary 
propagator Psi,({</? fe }, t, to) always commutes with these conditional unitary 
operations: 

P SL {{ Vk },t,t )U c b = U c b P SL {{ Vk },t,t ), (3) 

P SL (W k },t,t )U c f ee U c f P SL (W k },t,t ). (4) 

Though the conditional unitary operations and V£ may be independent of 
any atomic motional states according as their definitions, these two unitary 
operations may require that the wave-packet motional state \CM0, Rq) of the 
halting states {\n' h , CMO, R ),n' h = 0, 1} of the halting-qubit atom in the left- 
hand potential well keep unchanged up to a global phase factor when these 
two unitary operations are applied to the quantum system in the period from 
the initial time to of the quantum circuit to the time toi — St r before the state 
|1, CMO, Rq) is changed to the state \ci) = |1, CMl(toi), i?i(to«)) by the state- 
dependent trigger pulse P t . Here the quantum circuit Q c starts to run the first 
cycle at the initial time t . The requirement may be necessary when both the 
unitary operations are built up out of the dipole-dipole interactions or the 
Coulomb interactions which are dependent on the interdistances of atoms. 
Since the double-well potential field is considered as one of the components of 
the state-locking pulse field Psl in the atomic physical system, any motional 
state of the halting-qubit atom in the double-well potential field is affected 
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inevitably by the state-locking pulse field. However, as pointed out before, 
the intermediate potential barrier is so high and wide that up to a global 
phase factor the wave-packet motional state \CM0, Rq) of the halting states 
{\n' h , CMO, Ro)} is almost unchanged in the period from the initial time to to 
the time t 0i —5t r before the state |1, CMO, R ) is changed to the unstable state 
|ci). This property of the state- locking pulse field could be simply expressed 
by the unitary transformation: 

P SL (W k },t,t )\n' h ,CMO,R ) 

= exp[-iE (t - t )/h]\n' h , CMO, R ), t <t< t 0i - St r . (5) 

where n' h = 0, 1 and E is the motional energy of the ground motional 
state \CM0, R ) . Hereafter for convenience the global phase factor such as 
exp[— iE (t — t )/H] in Eq. (5) is omitted without confusion. Therefore, both 
the unitary operations and Vfi commute approximately with the unitary 
propagator Psx({<^ fc }, t, t ) (to < t < toi — St r ) of the state-locking pulse field, 

P SL (W k },t,to)U c h = U^P SL (W k },t,to), (6) 

P SL (W k },t,t )V£ = V^P SL {{ Vk },t,t ). (7) 

Moreover, the commutation relations (6) and (7) still hold when the halting- 
qubit atom enters into the right-hand potential well from the left-hand one, 
because in this case these two unitary operations U% and V£ have not any 
real effect on the halting-qubit atom and are reduced theoretically to the 
unity operation. Therefore, the commutation relations (6) and (7) hold for 
the whole quantum circuit. The commutation relations (5), (6), and (7) lead 
to that there hold the state unitary transformations, 

{P SL (W k },t,t ),U c h }\n h ,CMO,Ro)\b h )\f r (x)) 

= U c h \n' h ,CM0,R )\b h )\f r (x)), t < t < t 0l - St r , (8) 

{P SL (W k },t,t ),V^}\n' h ,CMO,R )\b h )\f r (x)) 

= V^\n' h ,CM0,R )\b h )\f r (x)), t <t<t 0i -St r . (9) 

Evidently, the unitary propagator Psz,({</? fe }, t, to) generally could not com- 
mute with the unitary operation Ujr r of the trigger pulse P t in the atomic 
physical system, 

p SL (W k },t,t )u? r ? uf r p SL (W k },t,t ). 



52 



This is because in the atomic physical system the state-dependent exci- 
tation process of the trigger pulse P t from the state \1,CM0, Rq) to the 
unstable state |1, CMl(ta), Ri(toi)) is affected inevitably by the left-hand 
harmonic potential field. Then the coherent Raman adiabatic laser trigger 
pulse P t must be designed suitably such that when both the trigger pulse 
P t and the state-locking pulse field Psl are applied simultaneously to the 
halting-qubit atom in the left-hand potential well the unitary propagator 
{Pshii^k}) toi, toi — St r ), P t } satisfies the relation: 

{P SL (W k },t 0l ,t 0t -5t r ),P t }\l,CM0,R }\b h }\f r (x)} 

= U^\l,CMO,Ro)\b h )\f r (x)) = \l,CMl(toi),Ri(t 0i ))\b h )\f r (x)). (10) 
The approximate calculation in theory for the unitary propagator U£ r = 
{PsL{Wk\-> toi-* ~~ fitr)>Pt} could be carried out on the simple physical 
model of the time- dependent forced harmonic oscillator [5c, 37, 50g, 50h, 
50i, 51b]. 

Suppose again that the six unitary operations Uj, Psz,({<^ fc }, Aio), U£ r , 
Vfi, Ufl, and U§ in the quantum circuit Q c have the durations Stf, At , 5t r , 
St' h , 5th, an d Stb, respectively, and the period of each cycle of the quantum 
circuit is AT = 5tb + 5th + 5t' h + 5t r + At^ + Stf. As stated before, the unstable 
state \ci) = \l,CMl(toi), Ri(toi)) of the control state subspace is generated 
completely at the instant of time toi in the i— th cycle of the quantum circuit. 
Then it follows that the functional state \f r (x)) is converted into the desired 
functional state \f r (%f)) during the period Stf from the time (t 0i — 8t r — St' h — 
5th — 5tb) — Stf to the time (t 0i — 5t r — 5t' h — 5th — Stb) in the (i — 1)— th cycle. 
Note that the initial state of the quantum circuit is the wave-packet state 
|0, CMO, R o )\0)\f r (x o )) and at the initial time t the halting-qubit atom is in 
the state |0, CMO, R ) and in the left-hand potential well. Then the desired 
functional state \f r {xf)) takes the integer Xf = (xq + i — 1) modm r . Before 
the i— th cycle the initial halting state \nh) = |0, CMO, Rq) and the initial 
branch-control state \bh) = |0) keep unchanged but only the initial functional 
state \ f r (x )) is consecutively changed to other functional state \ f r (x + j)) 
for j — 0, 1, 2, i — 1 in the quantum circuit, where the last functional state 
\f r (xo + % — 1)) is the desired functional state \f r (xf)}. Obviously, before 
the i— th cycle the time evolution process of the whole quantum system of 
the quantum computer with the initial state |0, CMO, Rq) |0) \f r {xo)) can be 
expressed as 

{Ps L ({ ¥ » fc },* j ,to),(C/r) J '}|0,CM0, J Ro)|0)|/ r (a;o)) 
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= (C/ r y|0,CMO, J R )|0)|/ r (a;o)) 

= \0,CM0,R o )\0)\f r (x +j)), t J= t +jAT, < j < i - 1, (11) 

where the unitary operator Ut is denoted as the unitary operation sequence 
UfPsL(Wk}i ^o)U^ r V h c U^U^ of the quantum circuit Q c . In the unitary trans- 
formation (11) the first equation shows that up to a global phase factor the 
unitary propagator of the quantum system under both the state-locking pulse 
field Psl and the unitary operation sequence (UtY for < j < i — 1 acting 
on the initial state is really equal to the single unitary operation sequence 
{UtY acting on the same initial state before the i— th cycle of the quantum 
circuit. Here the integer j is also used as the cyclic index of the quantum 
circuit. 

For example, there are some cases to be considered for the unitary trans- 
formation (11). For the first case (i) the initial functional state \f r (xo)) hap- 
pens to be just the desired functional state \f r (xf)). In this case the index 
% — 1. Then Xf = xq and the cyclic index takes j = in the unitary transfor- 
mation (11). Therefore, tj = to for j — 0. Here toi — + + + St' h + 5t r . 
Both the unitary propagator -Psz,({<£>fc}, tj, t ) and the unitary operation se- 
quence (UtY become the unity operator in effect as the cyclic index j = and 
the unitary transformation (11) is a state identity. For the second case (ii) 
the initial functional state |/ r (^o)) could be changed to the desired functional 
state \f r (xf)) at the end of the first cycle of the quantum circuit. In this case 
the index % = 2. Then Xf — (x + 1) modm r and the cyclic index takes j = 
and 1 in the unitary transformation (11). Note that the end time of the first 
cycle is just equal to the beginning time of the second cycle in the quantum 
circuit. Thus, to = to, t\ = t +AT = t i+At -|-5t/ = to2—5t r —5t' h —5th—Stb. 
Here toi = to + 5t& + 5th + St' h + 8t r and to2 = toi + AT. Generally, for 
the third case (Hi) the initial functional state \f r (xo)) could be changed 
to the desired functional state \f r (xf)) at the end of the (i — 1)— th cy- 
cle of the quantum circuit. In this case the cyclic index i > 1. Then 
Xf = (xq + % — 1) modm r and the cyclic index takes j = 0, 1, (i — 1) 
in the unitary transformation (11). Here the end time of the (i — 1)— th cycle 
is just the beginning time of the %— th cycle in the quantum circuit. Thus, 
t = t and tj = t + j AT = t j + At + Stf = to(j+i) — St r — 5t' h — 5t h — Stf, for 
j = 1, (i - 1). Here t 01 = t + 5t b + St h + St' h + 5t r and t o+i) = t 0i + j'AT 
for j = I, ■■■,{%- 1). 

The i— th cycle (i > 1) of the quantum circuit will be analyzed separately 
as follows. In the i— th cycle of the quantum circuit the starting time is 



54 



t(j-i) — t + (i — 1)AT = t 0i — 5t r — 8t' h — 5th — 8h an d the starting state 
|0, CMO, R o )\0)\f r (xf)) with Xf = (x + i — l)modm r . Suppose that now 
the quantum circuit starts to execute the %— th cycle. First, the unitary 
transformation U£ changes the initial branch-control state \bh) = |0) to the 
state |1), 

{P 5L (M},V 1) + 5t 6 ,V 1) ),f/ 6 c }|0,^M0,i? )|0)|/ r ( a ; / )) 

= P SL (W k }, t (i -i) + St b , t (i _i))^|0, CMO, R )\0)\f r (x f )) 

= \0,CM0,R o )\l)\f r (x f )) . (12) 

Here the unitary propagator {-Psl({ + <^&> i s identical 

to the unitary operation Ps L ({</? fe }, fy-i) + 5t b ,t^i))U^ due to the fact that 
both the unitary transformation -Psz,({<^fc}, + of the state- 

locking pulse field and the unitary operation U£ commute with each other, 
as shown in Eq. (3). The second equation holds due to the state unitary 
transformation P SL ({ip k }, t (i _i } + 5t b , fy-i))|0, CMO, R ) = |0, CMO, R ), as 
shown in the unitary transformation (5). Then, the unitary operation 
converts the desired functional state \f r (xf)) = |1) into the state |0), 

{PsdWk}^h-i) + ttb)M}\^CMO,R )\l)\f r {x f )) 

= PsL(Wk}^^(i-i) + 5t b )Uh\0,CM0,Ro)\l)\f r (xf)) 

= \0,CM0,R )\1)\0), t = t {i - 1) + 6t b + 6t h , (13) 

and the unitary operation V h c further changes the initial halting state \nh) = 
\0,CM0,R ) to the state \1,CM0,Rq), 

{P SL (Wk}, t, + St b + 5t h ), V h c }\0, CMO, R )\l)\0) 

= P SL (W k },t,t {l „ 1) + 5t b + 5t h )V h c \0,CMO,Ro)\l)\0) 

= |l,CM0,i2o>|l>|0>, (14) 

where t = i(i-i) +5t b +5th + St' h = toi — 5t r . These two unitary transformations 
(13) and (14) may be obtained by the relations (8) and (9), respectively, and 
here equation (5) has also been used. Now the state-dependent trigger pulse 
P t starts to act on the halting state \nh) = |1, CMO, R ) at the time t 0i — St r 
due to that the halting state has the internal state with n' h — 1 and the 
halting-qubit atom now is in the left-hand potential well. With the help of 
the relation (10) the state-dependent excitation process of the trigger pulse P t 
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from the halting state to the unstable state \ci) of the control state subspace 
may be expressed as 

8t r ),U? r }\l,CM0,R o )\l)\0) 

= \l,CMl(toi),Ri(toi))\l)\0) (15) 

where t 0i = t(i^i)+5t b + 5t h + 5t' h + 5t r and {P 5L ({y? fc }, t i, t 0i -5t r ), U£ r } of the 
quantum circuit Q c is just defined as the unitary operation U^ r (see Eq. (10)). 
This excitation process increases the motional energy of the halting-qubit 
atom in the left-hand potential well so that the following unitary quantum 
scattering process for the atom can take place automatically. 

From the instant of time toi on, the quantum circuit Q c really starts to ex- 
ecute simultaneously its own two almost independent processes: the quantum 
control process and the quantum computational process. The quantum con- 
trol process could be really thought to start at the time t 0i and at the initial 
state |ci) = |1, CMl(i 0i ), Ri(t 0i )) for the physical system of the halting-qubit 
atom. Here for convenience the state-dependent excitation process of the 
trigger pulse is not included in the quantum control process. In the quantum 
control process the halting-qubit atom with the state |1, CMl(toi), Ri(toi)) 
first carries out a unitary quantum scattering process in space and time. Dur- 
ing the period At of the quantum scattering process the halting-qubit atom 
is driven by its own motional momentum to leave the left-hand potential 
well and pass the intermediate potential barrier to enter into the right-hand 
potential well. This quantum scattering process may be expressed formally 
by the unitary transformation: 

{P5L({^J,^,^), J P5L({^},Ato)}|l,C'Ml(^), J Rl(^))|l)|0) 

= P SL ({ Vk },At )\l,CMl(t 0t ),R 1 (t 0t )}\l}\0} 

= \l,CM2(t mi ),R 2 (t mi ))\l)\0), t mi = t 0i + At , (16) 

where {P 5L ({y? fc }, t mi , t 0i ), P 5L ({y9 fe }, At )} in the quantum circuit Q c is de- 
fined as Psx({<£> fe }, Ato), and the duration At must be long enough so that by 
the quantum scattering process the wave-packet state |1, CMl(t 0i ), Ri(t i)) 
of the halting-qubit atom in the left-hand potential well can be almost com- 
pletely transferred to the wave-packet state |1, CM2(t mi ), P 2 (t m j)) of the 
atom in the right-hand potential well. Actually, the duration At must en- 
sure that the wave-packet spatial distance \R2(t m i) — Pol is large enough such 
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that both the wave-packet state |1, CMO, R ) and |1, CM2(t mi ), R-2it mi )) are 
almost orthogonal to each other. Note that the distance \R2{t m i) ~ Ri(toi) \ is 
longer than \R 2 {t m i) — i?o I • This means that such a duration Ato ensures that 
both the wave-packet states |1, CMl(t j), -Ri(^oi)) and |1, CM2(t mi ), -fi^mi)) 
are also almost orthogonal to each other. The quantum scattering process 
(16) is dependent only upon the time difference At rather than the in- 
stant of time t mi or t 0i , that is, P SL (Wk}, tmi, t i) = PsUWh), Af ), and 
it is also an energy conservative process that the motional energy of the fi- 
nal state \c' 2 ) = \l,CM2(t m i), R2(t mi )) is really equal to that of the initial 
state |ci) = \l,CMl(t i), Ri(toi)). After the quantum scattering process the 
quantum computational process continues to execute the functional unitary 
operation Vj as usual, but from the time t m i on, the quantum computa- 
tional process is really halted in effect and meanwhile the halting-qubit atom 
also motions continuously in the right-hand potential well under the state- 
locking pulse field. This process may be expressed in terms of the unitary 
transformation: 

{P 5L (M},t mi + %,t mj ),f/ / c }|l,C'M2(t mi ), J R 2 (t™))|l)|0) 

= P 5L ({ V 9j,t mi + 5t / ,t mi )f/ / c |l,CM2(t mi ), J R 2 (t mi ))|l)|0) 

= \n' h (t mi + 5t f ), CM2(t mt + 8t f ), R 2 (t mi + %))|1)|0), (17) 

where \n' h (t mi +5tf)) is the internal state of the halting-qubit atom at the time 
t m i + Stf. The functional operation Wj does not have a real effect on the state 
\f r (x)) = |0) also due to the branch-control state \bh) = |1). Here it must be 
ensured that the halting-qubit atom is in the right-hand potential well during 
the period from the time t mi to the time t mi + 5tf so that the wave-packet 
state \n' h (t),CM2(t), R 2 (t)) with t mi < t < t mi + 5tf is almost orthogonal 
to any one of the three wave-packet states \n' h , CMO, R ) (n' h = 0, 1) and 
|1, CMl(toi), Ri(t 0i )). Obviously, here the width of the intermediate potential 
barrier should be large enough so that a wave-packet state of the halting- 
qubit atom in one potential well is almost orthogonal to any wave-packet 
state of the atom in another potential well in the double-well potential field. 

From the (i + 1)— th cycle to the end of the quantum circuit the condi- 
tional unitary operations U£, Vfi, U^ r , and Uj do not really affect the 
quantum system, although according as the quantum circuit these unitary 
operations are still applied continuously to the quantum system of the quan- 
tum computer, and only the state-locking pulse field takes a real action on 
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the halting-qubit atom in the right-hand potential well. Then the time evo- 
lution process of the quantum system from the time ti = t m i + 8tf to the end 
of the quantum computational process could be generally written as 

{PsL({<Pk},ti+j,U), (UT) j }\n' h (t i ),CM2(t i ),R 2 (t i ))\l)\0) 

= PsLiWkhti+jMrihiU), CM2(t i ),R 2 (t i ))\l)\0) 

= \n' h (t i+j ), CM2(t i+j ), R 2 (t i+j ))\l)\0). (18) 

Here the time ti+j = t« + jAT for < j < m r — i and the end time of the 
computational process is given by t mr = t + m r AT. From Eq. (11) and Eq. 
(18) one sees that before the initial halting state \0,CM0,Rq) is changed 
to the state |1,CM0, R ) the unitary propagator {P S L(Wk}i tj, t ), (UtY} 
acting on the initial state |0, CMO, R )\0)\f r (x )) is really equal to the single 
unitary operation sequence (UtY acting on the same initial state (see Eq. 
(11)), but after the halting-qubit atom enters into the right-hand potential 
well the unitary propagator {PsL({<£k}i U+j, ti), (UtY} acting on the starting 
state \n' h (ti), CM2(tj), R 2 (ti))\l)\0) is really equal to the single unitary prop- 
agator PsL{{(fk\-> U+j, U) °f the state-locking pulse field acting on the same 
state (see Eq. (18)). Here any wave-packet state \n' h (t), CM2(t), R 2 (t)} for 
ti <t < t nw of the halting-qubit atom in the right-hand potential well must 
be orthogonal or almost orthogonal to any one of the three wave-packet states 
\n' h ,CM0,R ) (n' h = 0, 1) and |1, CMl(i w ), Ri(t 0i )) of the atom in the left- 
hand potential well. Therefore, the quantum program and circuit requires 
that the halting-qubit atom be in the right-hand potential well from the time 
t mi to the end of the quantum computational process. 

Though the halting-qubit atom must be in the right-hand potential well 
at the end of the computational process no matter what the initial functional 
state \f r (xo)) is with xo = 0, 1, ...,m r — 1, each possible wave-packet state of 
the halting-qubit atom \n' h (t mr ), CM2(t mr ), R 2 (t mr )) with t mr = U + (m r - 
i)AT for i — 0,1, ...,m r — 1 could be different in spatial position in the right- 
hand potential well at the end of the computational process. This is because 
a different initial functional state \f r (xo)) corresponds to a different wave- 
packet state \n' h (t mr ), CM2(t mr ), R 2 (U + (m r — i)AT)) which is located at a 
different spatial position R 2 (ti + (m r — i)AT) in the right-hand potential well. 
Here n' h (t mr ) and CM2(t mr ) with t mr = ti + (m r — i)AT each takes a same 
value for all different index values i — 0,1, ...,m r — 1, respectively. Actually, 
similar to the situations in the conventional halting protocol [11, 54], all these 
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m r possible wave-packet states {\n' h (t mr ) , C M2(t mr ) , R 2 (U + (m r — i)AT))} 
for i — 0,1, m r — l should be almost orthogonal to one another. Therefore, 
the output wave-packet states of the quantum circuit at the end of the com- 
putational process are dependent sensitively upon the initial functional states 
{ | /r (#()))}■ However, as pointed out before, the quantum circuit Q c will not 
be a suitable component of the quantum search processes based on the uni- 
tary quantum dynamics [24] if its output state is dependent sensitively upon 
any initial functional state \f r (xo))- Evidently, if each of these m r possible 
wave-packet states {\n' h (t mr ), C M2(t mr ) , i? 2 (^i + ijn r — i)AT))} at the end of 
the computational process can be further transferred to some desired state 
in a high probability close to 100% by a given unitary transformation, then 
the output state of the quantum circuit could be considered to be almost 
independent of any initial functional state. Of course, it is impossible that 
the unitary transformation can change all these m r wave-packet states to the 
same desired state in the probability 100%. It could be better to choose the 
desired state as the wave-packet state \nh) = \n' h , CMO, R ) (n' h = or 1) 
of the halting-qubit atom in the left-hand potential well as the wave-packet 
state is stable in the double-well potential field. Thus, it is necessary to ma- 
nipulate and control coherently the halting-qubit atom by the state-locking 
pulse field after the halting-qubit atom enters into the right-hand potential 
well. The coherent manipulation process have two purposes in the quantum 
control process to simulate efficiently the reversible and unitary halting pro- 
tocol. The first one is that the halting-qubit atom can stay in the right-hand 
potential well for a long time till the computational process finished after 
the halting-qubit atom enters into the right-hand potential well. The second 
is that after the computational process finished the halting-qubit atom can 
return the left-hand potential well from the right-hand one and the returning 
halting-qubit atom in the left-hand potential well is in the wave-packet state 
\n' h , CMO, Rq) (n' h = or 1) with a high probability close to 100% no matter 
what the wave-packet state \n' h (t mr ),CM2(t mr ), R 2 (ti + (m r — i)AT)) is for 
i — 0,1, ...,m r — 1. The coherent manipulation of the halting-qubit atom in 
the right-hand potential well generally starts after the computational process 
finished, but actually this manipulating process may start at a much earlier 
time t mi when the halting-qubit atom enters into the right-hand potential 
well rather than at the end of the computational process. 

5. Coherently manipulating the halting-qubit atom in time and 
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space 

In this section the issues to discuss are focused on how the state-locking 
pulse field manipulates coherently the wave-packet states of the halting-qubit 
atom in the right-hand potential well of the double- well potential field so that 
the halting-qubit atom can stay in the right-hand potential well for a long 
time till the computational process finished and how the quantum control 
process that simulates efficiently the reversible and unitary halting proto- 
col in the atomic physical system makes its output state insensitive to any 
initial functional state of the quantum circuit Q c . As shown in the previ- 
ous section, when the halting-qubit atom leaves the left-hand potential well 
and enters into the right-hand one at the time t mi = t 0i + At in the i— th 
cycle of the quantum circuit, the conditional unitary operations Ufa and Vfi 
and the state-dependent trigger pulse Pt of the quantum circuit become un- 
available. Thus, from the time t m i on, the quantum computational process 
really terminates in effect, although it runs continuously to the end of the 
quantum circuit. Meanwhile, the halting-qubit atom evolves continuously in 
the right-hand potential well under the state-locking pulse field. Note that 
the state-locking pulse field consists of the double-well potential field itself 
and the sequences of the time- and space-dependent electromagnetic field 
pulses which are applied only to the right-hand potential well during the 
computational process and then to the double-well potential field after the 
computational process finished. Obviously, these sequences of the time- and 
space-dependent electromagnetic field pulses of the state-locking pulse field 
have a negligible effect on any state of the halting-qubit atom when the atom 
is in the left-hand potential well during the computational process. Thus, it 
follows from the unitary transformations (17) and (18) that, from the time 
t mi on, the time evolution process of the whole quantum system of the quan- 
tum computer may be reduced to the simpler quantum control process of the 
halting-qubit atom in the right-hand potential well under the state-locking 
pulse field, 

PsdWk}, t, t mi )\l, CM2(t mi ), R 2 (t mi )) 

= \n' h (t),CM2(t),R 2 (t)), t mi <t. (19) 

Here the wave-packet state \n' h (t),CM2(t), R 2 (t)} (t mi < t) generally could 
be expanded in terms of the motional basis states of the halting-qubit atom 
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in the right-hand potential well [38-43] , 

\n' h (t),CM2(t),R 2 (t)) = E « CM2,t)\n' h )\CM2,R 2 ). 

n' h CM2 

Now the first purpose of the quantum control process is to design a unitary 
sequence of the time- and space-dependent electromagnetic pulses and/or 
the shaped potential fields of the state-locking pulse field PsL({<fk}i t, t m i) 
to manipulate the halting-qubit atom so that the atom is able to stay in the 
right-hand potential well till the end of the computational process. For con- 
venience, suppose that the lowest point of the left-hand harmonic potential 
well is equal to the bottom of the right-hand square potential well and both 
are set to zero. As stated before, when the halting-qubit atom is in the unsta- 
ble wave-packet state |1, CMl(t i), Ri(toi)) the total motional energy of the 
atom which includes the kinetic and potential energies in the left-hand har- 
monic potential well but not the atomic internal energy is much higher than 
the height of the intermediate potential barrier in the double-well potential 
field. Since the quantum scattering process is an energy-conservative pro- 
cess the total motional energy of the halting-qubit atom remains unchanged 
when the halting-qubit atom enters into the right-hand potential well from 
the left-hand one, but it is completely converted into the kinetic energy as the 
potential energy of the halting-qubit atom is zero in the right-hand square po- 
tential well. Suppose that the relativistic effect is negligible for the motional 
halting-qubit atom in the double-well potential field. Then the motional ve- 
locity of the halting-qubit atom is given by Vh = \/2E h /m h at the time t mi 
after the atom enters into the right-hand potential well, where and 
are the total motional energy and mass of the halting-qubit atom, respec- 
tively. Therefore, the motional velocity of the halting-qubit atom could 
become very large when the atom enters into the right-hand potential well 
from the left-hand one. Since the geometric length of the right-hand square 
potential well is limited it is impossible to keep the halting-qubit atom in 
the right-hand potential well for a long time required by the computational 
process if the motional velocity is very large. Thus, the motional velocity 
of the halting-qubit atom must be decelerated greatly by the state-locking 
pulse field so that the halting-qubit atom does not leave in a short time the 
right-hand potential well. This decelerating process could be achieved by the 
unitary sequence of the time- and space-dependent laser light pulses of the 
state-locking pulse field. The decelerating process of the halting-qubit atom 
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is really very similar to the conventional atomic laser cooling processes [36]. 
The interactions between the halting-qubit atom in motion and the laser light 
pulse field become important in the decelerating process, while it is known 
that the dipole force of the laser light pulse field exerting a motional atom 
plays an important role in the atomic laser cooling processes [36b] . There are 
a variety of the atomic laser cooling methods and techniques which have been 
discovered and developed in the past decades and used extensively to cool 
an atomic ensemble to an extremely low temperature [36] , but most of these 
atomic laser cooling methods are non-unitary. Thus, not all these atomic 
laser cooling methods are suitable for building up the unitary decelerating 
sequence because the decelerating process for the halting-qubit atom must 
be a unitary process. This is quite different from a conventional atomic laser 
cooling process in an atomic ensemble. Only when an atomic laser cooling 
method is unitary or can be made unitary can it be exploited to decelerate the 
halting-qubit atom in the quantum control process. Thus, only the coherent 
atomic laser cooling methods could be used to build up the unitary deceler- 
ating sequence. Another difference between an atomic laser cooling process 
and the unitary decelerating process is that the unitary decelerating process 
is simpler as the current atomic physical system is the pure quantum-state 
system of an atomic ion or a neutral atom in the double-well potential field, 
while it is usually more complex to cool an atomic ensemble by an atomic 
laser cooling technique. The third difference is that it could not be necessary 
to decelerate the halting-qubit atom to zero velocity, while the target of a 
conventional laser cooling technique is to cool an atomic ensemble to an ex- 
tremely low temperature as close to zero degree as possible. Therefore, the 
unitary decelerating process for the halting-qubit atom is relatively easy to 
be achieved by a coherent atomic laser cooling technique. 

The mechanisms for the laser cooling in an atomic ensemble have been 
studied extensively and thoroughly in the past years [36]. However, in or- 
der to investigate the mechanism of the unitary decelerating process for the 
halting-qubit atom here the basic atomic-laser-cooling mechanism is intro- 
duced briefly. A conventional atomic laser cooling method usually consists of 
both the atomic optical-pumping process (or the atomic optical-absorption 
process) and the atomic optical-emission process [36c]. The optical pump- 
ing process is that an atom under cooling is excited from the ground state 
to the excited state by absorbing photons from the laser light field, while 
the optical emission process is that the atom in the excited state returns 
the ground state by emitting photons to the laser light field or environment. 
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The optical pumping process generally may be made unitary easily, but the 
optical emission process usually is simply chosen as a spontaneous and ran- 
dom process in a conventional laser cooling technique partly due to that the 
lifetime of the atomic excited state usually is very short. Hence the optical 
emission process usually is a non-unitary process in an atomic ensemble for a 
conventional laser cooling technique. Suppose that a free atomic ion or neu- 
tral atom with the mass m in the ground state is irradiated by a laser light 
field and makes a transition from the ground state to the excited state by 
absorbing a photon from the laser light field. For convenience, here only non- 
relativistic limit is considered for the atomic motion. The optical pumping 
process generally obeys the energy-, momentum-, and angular momentum- 
conservative laws. Denote that before the transition the atom in the ground 
state has the internal-state energy E a , the kinetic energy p a /(2m), the mo- 
mentum p a = mv„, and the angular momentum J a , respectively, and the 
photon that will be absorbed by the atom has the photonic energy E c = hu, 
the momentum p c = Kk (\p c \ = fku/c), and the angular momentum J c , re- 
spectively. After the transition the photon is absorbed by the atom and 
hence the photonic energy, momentum, and angular momentum are trans- 
ferred to the atom. The atom now is in the excited state. Suppose that after 
the transition the atom in the excited state has the internal-state energy E b , 
the kinetic energy E b = p\j (2m), the momentum p b = rnvj, and the angular 
momentum J b , respectively. Then the energy conservation before and after 
the transition shows that there holds the relation: 



where on is the photonic frequency. The momentum conservation leads to 
that 



where k (|k| —oj/c) is the wave vector of the photon before the transition 
and v a and v b are the motional velocity vectors of the atom before and after 
the transition, respectively. The angular momentum conservation between 
the angular momentum J a of the atom and the photonic angular momentum 
J c before the transition and the angular momentum J b of the atom after the 
transition will not be discussed in detail here. If the optical pumping process 
(or the photon absorption process) for the atom is one-dimensional, then the 
wave vector k of the photon is either co-direction to the motional velocity v a 
of the atom or opposite to the velocity v a before the transition. For the first 



E a + p 2 J(2m) + fku = E b + p 2 b /(2m) 



(20) 



mv„ + Kk = mv b , 



(21) 



63 



case that both the wave propagating of the laser light field with the vector 
k and the motion of the atom with the velocity v a are co-direction to each 
other the motional velocity of the atom after the transition is given by 

Vb = v a + Hk/m > v a , (22) 

hence the motional atom is accelerated by the copropagating laser light field, 
while for the second case that the laser light field propagates in the oppo- 
site direction to the atomic motion the atomic motional velocity after the 
transition is written as 

v b = v a - Hk/m < v a , (23) 

hence the motional atom is decelerated by the opposite propagating laser 
light field. Obviously, here the atom is slowed down by Hk/m when the atom 
absorbs a photon from the opposite propagating laser light field. When the 
atom is in the excited state it no longer absorbs any photons from the oppo- 
site propagating laser light field and hence can not be further slowed down. 
One of the schemes to decelerate further the atom is that the atom in the 
excited state first jumps back to the ground state without changing signifi- 
cantly its total motional momentum, and then it absorbs a photon from the 
opposite propagating laser light field again and hence is decelerated further. 
Note that the process that the excited-state atom jumps back to the ground 
state by emitting photons to the laser light field or its environment is just 
the atomic optical-emission process. Thus, the atomic laser cooling process 
consists of a number of the atomic optical absorption-emission cycles that the 
atom in the ground state absorbs a photon to make a transition to the excited 
state and then jumps back to the ground state from the excited state by emit- 
ting photons. The optical emission process is either a spontaneous emission 
process in a random form or the stimulated emission process in a coherent 
form. In the spontaneous emission process photons are emitted in a random 
form by the atoms in the excited state so that the total momentum of the 
emitting photons is zero and hence the atomic motional momentum does not 
change significantly after the emission process. Therefore, the atom is slowed 
down basically by the optical pumping process in every optical absorption- 
emission cycle if the optical emission process is spontaneous and random. 
Most of the conventional atomic laser cooling methods and techniques use 
the spontaneous optical-emission process as their key component to cool an 
atom ensemble. Every optical absorption-emission cycle can make the atom 
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to be slowed down by Hk/m and the atom can be slowed down continuously 
by a sequence of the optical absorption-emission cycles. However, the spon- 
taneous optical-emission process of the conventional laser cooling methods 
is not allowed due to its own non-unitarity if these laser cooling methods 
are used to slow down the halting-qubit atom in the quantum control pro- 
cess. On the other hand, the coherent optical-emission process is different 
from the spontaneous optical-emission process in that the coherent optical- 
emission process may be a unitary process. The coherent optical-emission 
process generally may be stimulated by an external laser light field [5a]. The 
momentum of the emitting photons from the atom in the excited state in the 
coherent optical-emission process is not zero on average and hence can make 
a significant contribution to the motional momentum of the atom after the 
emission process. The coherent atomic optical-emission process still obeys 
the energy-, momentum-, and angular momentum-conservative laws. There- 
fore, if the emitting photon travels along the same direction to the atomic 
motion, then the atom will lose part of its motional momentum after the 
emission process and hence is slowed down. In the optical-emission process 
not only the atomic internal energy (E^ — E a ) of the excited state is trans- 
ferred to the photonic energy but also part of the motional energy of the 
atom in the excited state is converted into the photonic energy. However, 
the atom will receive a recoil momentum from the emitting photon and hence 
is accelerated after the emission process if the emitting photon travels along 
the opposite direction to the atomic motion. In this atomic optical-emission 
process the atomic internal energy of the excited state is transferred partly 
to the photonic energy and partly to the atomic kinetic energy at the same 
time. In the quantum control process the halting-qubit atom in the right- 
hand potential well must be first slowed down greatly so that it is able to 
stay in the right-hand potential well for a long time till the end of the compu- 
tational process in the quantum circuit Q c , and then the atom is sped up in 
a unitary form after the computational process finished such that the atom 
can return to the left-hand potential well. The coherent atomic laser cooling 
methods and techniques therefore provide a possible way to generate both 
the decelerating and accelerating processes for the halting-qubit atom in the 
quantum control process. 

A conventional laser cooling method based on the optical absorption- 
emission (spontaneous) cycles usually is realized more easily than a coherent 
one in an atomic ensemble. In general, the spontaneous optical-emission 
process from the atomic excited state to the ground states in the atomic en- 
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semble occurs easily in nature as the atomic excited state usually has a much 
shorter lifetime than those ground states of the atom. This in turn implies 
that a coherent atomic laser cooling method could be more complex as the 
non-unitary spontaneous optical emission must be avoided in the coherent 
laser cooling process. The stimulated Raman adiabatic passage (STIRAP) 
laser cooling method is one of the important atomic laser cooling methods. 
The STIRAP method has been used extensively to cool an atomic ensem- 
ble to an extremely low temperature [42, 43], to cool a trapped atomic ion 
to the ground state for quantum computation [30, 44], to manipulate a co- 
herent atomic beam in the atomic interferometry [38, 39, 40, 41], and to 
prepare and manipulate a nonclassical motional state in a trapped-ion phys- 
ical system [37, 45]. In particular, the STIRAP laser cooling method could 
be used conveniently to cool a multi-level atomic ensemble with many inter- 
nal states to an extremely low temperature. The coherent STIRAP laser 
cooling (or decelerating) method could be a better candidate to avoid the 
non-unitary spontaneous optical emission of the atom from the excited state 
to the ground states. This is because the cooling (or decelerating) atom does 
not stay in the excited state at all or could stay in the excited state in a 
much shorter time than the lifetime of the excited state during the coher- 
ent STIRAP laser cooling process if the Raman adiabatic laser pulses are 
detuned properly from the excited state. Since an adiabatic laser beam usu- 
ally has a much wider frequency bandwidth than a conventional CW laser 
beam the STIRAP laser cooling method is able to take the Doppler effect 
into account conveniently during the atomic laser cooling process. It is well 
known that a general Raman adiabatic laser pulse consists of a pair of the 
adiabatic laser beams with the specific characteristic parameters. Generally, 
the characteristic parameters for the adiabatic laser beams of a Raman adia- 
batic laser pulse include the carrier frequencies and detunings, the frequency 
bandwidths, the amplitudes and phases of the adiabatic laser light fields, the 
laser-beam durations, the propagation directions and polarizations (e.g., <t+ 
or <t_), and the spatial action positions and zones. Suppose now that the 
states \gi), \g2), and |n e ) are three different internal states of the cooling 
atom and their corresponding wave-packet states of the atom are written 
as \gi,CMi, Ri), \g2,CM 2 , R 2 ), and \n e , CM e , R e ), respectively. These two 
internal states \g±) and \g 2 ) usually may be chosen as a pair of ground inter- 
nal states or two lowest energy-level internal states of the atom, while the 
internal state \n e ) may be an excited state whose energy level is much higher 
than those of the internal states \g±) and \g 2 ). In the coherent STIRAP laser 
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cooling method an adiabatic laser beam A may be used to pump the atom 
from the ground internal state \gi) to the excited state \n e ) and at the same 
time another adiabatic laser beam B is applied to stimulate the atom in the 
excited state \n e ) to jump back to the internal state \g 2 ) [42, 43, 44]. The 
coherent ST I RAP process may be formally expressed in terms of the unitary 
transition process: 

IguCM^Ri) A \n e ,CM e ,R e ) £ \g 2 ,CM 2 ,R 2 ). 

Obviously, the carrier frequency of the adiabatic laser beam A should be 
close to the resonance frequency between the ground state \gi) and the ex- 
cited state \n e ), while the carrier frequency for the adiabatic laser beam B 
is close to the resonance frequency between the ground state \g 2 ) and the 
excited state \n e ). In order to avoid occurring the non-unitary spontaneous 
optical emission for the atom from the excited state \n e ) to the ground states 
both the adiabatic laser beams A and B are detuned properly from the ex- 
cited state \n e ). For example, if the unitary state transfer \gi) <-> \g 2 ) is 
achieved by the conventional CW laser light irradiation, that is, the internal 
state \gi) is first transferred completely to the excited state \n e ) and then 
to \g 2 ) by the CW irradiation method, then the decoherence effect usually 
affects largely the state transfer since lifetime of the excited state \n e ) usu- 
ally is very short and also much shorter than those of the ground states \gi) 
and \g 2 ), whereas the ST I RAP method can avoid such decoherence effect 
on the state transfer. While the direct transition from the ground state \gi) 
(\g 2 )) to another ground state \g 2 ) (\g 2 )) is prohibited under the CW laser 
light irradiation, the coherent STIRAP method is a better scheme to ex- 
cite indirectly the transition between these two ground states. Thus, the 
coherent STIRAP method has some advantages over the conventional CW 
irradiation method to transfer the ground state \gi) (\g 2 )) to another ground 
state \g 2 ) (Ifi^)) in a unitary form. Obviously, the effective spatial bandwidth 
(ESB) for the Raman adiabatic laser pulse must be greater than the spatial 
displacement (SD) of the atom during the Raman adiabatic laser pulse. The 
spatial displacement (SD) is not more than the pulse duration {t p ) of the 
Raman adiabatic laser pulse times the maximum velocity (%) of the atom 
during the Raman adiabatic laser pulse, that is, SD < v M x t p < ESB. 

While the ground internal state \gi) (\g 2 )) is completely transferred to 
the state \g 2 ) (\g\)) by the Raman adiabatic laser pulse, the corresponding 
motional state \CMi, Ri) (\CM 2 , R 2 )) of the atom is also changed to another 
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motional state \CM 2 ,R 2 ) (\CM 1 ,R 1 )). During the coherent STIRAP pro- 
cess the atom could be either sped up or slowed down and this is mainly 
dependent upon the characteristic parameter settings for these two adiabatic 
laser beams A and B of the Raman adiabatic laser pulse and also the initial 
atomic motional velocity and direction, as mentioned earlier. An example is 
given below. Suppose that the atom is in the wave-packet state |gi, CMi, Ri) 
at the initial time, the propagating directions of these two adiabatic laser 
beams A and B are opposite to each other, and the beam A propagates in 
the opposite direction to the atomic motion. Then the atom will be slowed 
down by HkA/nT> + Hks/m when the wave-packet state \g±, CMi, R\) is trans- 
ferred to \g 2 ,CM 2 , R 2 ) by the Raman adiabatic laser pulse [43]. Here sup- 
pose that the initial atomic velocity is much greater than Hh^/m + fik B jm. 
This decelerating process could be understood intuitively: {%) when the state 
| #1, CMi, Ri) is induced a transition to the excited state \n e , CM e , R e ) by the 
laser beam A the atom is slowed down by hkA/m because the atom absorbs 
the photonic momentum Hk^ from the laser light field of the beam A which 
travels along the opposite direction to the atomic motion; (ii) when the atom 
is stimulated by the laser beam B to jump to the state \g 2 ,CM 2 , R 2 ) from 
the excited state \n e , CM e , R e ) it releases the momentum hk B to the laser 
light field of the beam B and the atom therefore is slowed down further 
by hk B /m as the atomic motional direction is the same as the propagating 
direction of the beam B. Evidently, the atom can also be sped up when 
the atomic wave-packet state \g\, CMi, Ri) is transferred to \g 2 , CM 2 , R 2 ) by 
the Raman adiabatic laser pulse with the proper characteristic parameter 
settings. Furthermore, the atom may be slowed down or sped up continu- 
ously by many Raman adiabatic laser pulses with the proper characteristic 
parameter settings. For example, suppose that one wants the atom to be 
decelerated further after the atom is slowed down by hk a/ 'm + hks / 'm by the 
Raman adiabatic laser pulse R(A, B) with the beams A and B. Then one 
may apply another Raman adiabatic laser pulse R(A±, B\) with the beams A\ 
and B\ to the state \g 2 , CM 2 , R 2 ) to decelerate further the atom. Since both 
the spatial positions R\ and R 2 are different for these two wave-packet states 
|<7i, CMi, Ri) and \g 2 , CM 2 , R 2 ) the applying spatial position (R 2 ) of the Ra- 
man adiabatic laser pulse R(Ai, B\) is different from that one (i?i) of the Ra- 
man adiabatic laser pulse R(A, B). Here the adiabatic laser beam Ai should 
travel along the opposite direction to the atomic motion, while the beam Bi 
propagates in the opposite direction to the beam A\. Then the atom is decel- 
erated by HkA 1 /m+ Hk Bl /m again after the wave-packet state \g 2 , CM 2 , R 2 ) 
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is transferred to another state \g±, CM 3 , R 3 ) by the Raman adiabatic laser 

pulse : \g 2 ,CM 2 ,R 2 ) <H \n e ,CM' e , R' e ) h CM 3 , i? 3 ). The uni- 

tary decelerating (or accelerating) sequence of the state-locking pulse field 
used to decelerate (or accelerate) the halting-qubit atom in the right-hand 
potential well is built up out of these coherent Raman adiabatic laser pulses 
with the proper characteristic parameter settings. Obviously, both the uni- 
tary decelerating and accelerating sequences are time- and space-dependent. 
The conversion efficiency from one internal state (|<7i)) to another internal 
state {\g 2 )) measures the performance of a coherent Raman adiabatic laser 
pulse. A good-performance Raman adiabatic laser pulse should be able to 
convert completely the internal state \gi) (\g 2 )) into the state \g 2 ) (l^i))- It 
has been shown theoretically [46, 38] that in a three-state atomic system 
a ground internal state \gi) {\g 2 )) can be transferred completely to another 
ground state \g 2 ) {\g\)) in a unitary form through the excited state \n e ) by 
a Raman adiabatic laser pulse with the proper characteristic parameter set- 
tings. 

The halting-qubit atom generally may be chosen as a multi-level atom 
with many internal states in addition to the two internal states {|^), rih = 
0, 1} of the halting quantum bit. Now the coherent STIRAP method is 
used to manipulate the halting-qubit atom after the atom enters into the 
right-hand potential well. Here in the STIRAP method these two internal 
states {\n' h ), n' h = 0, 1} could be conveniently set to the internal states \gi) 
and \g 2 ), respectively, and the excited state \n e ) to some specific excited 
electronic state of the halting-qubit atom. A unitary decelerating sequence 
UD{Wk\->tmi + Ti>,t m i) consisting of the Raman adiabatic laser pulses with 
the proper parameter settings then is constructed to decelerate the 

halting-qubit atom when the atom enters into the right-hand potential well 
at the time t m i, where To is the total duration of the unitary decelerating 
sequence. The total duration Tp must be much shorter than the period AT 
of each cycle of the quantum circuit. Note that there are m r possible different 
times {t mi , i = 1,2, ...,m r } for the quantum circuit Q c . The halting-qubit 
atom may enter into the right-hand potential well at any time t m k of these 
m r possible times {t m i} for k = 1, 2, m T . In order to decelerate the halting- 
qubit atom the unitary decelerating sequence must be applied at every time 
t mi of these m r possible times {t mi } in the quantum circuit. As known 
in the previous sections, the halting-qubit atom completely enters into the 
right-hand potential well at the time t m i = toi + At in the i— th cycle of 
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the quantum circuit and at the time t m i the halting-qubit atom is in the 
wave-packet state |c' 2 ) = \l,CM2(t m i), R 2 (t m i)). Then the unstable wave- 
packet state |c' 2 ) will be changed to the stable wave-packet state \c2{t m i + 
To)) of the control state subspace after the unitary decelerating sequence 
[/^({^fc}? t m i + Td, t m i) acts on the halting-qubit atom at the time t m i in the 
right-hand potential well, 

\c 2 (t mt + T D )} = U D ({<p k },t mi + T D ,t mi )\l,CM2(t mi ),R 2 (t mi )) 
= \0,CM2(t mt + T D ),R 2 (t mt + T D )), 

meanwhile the initial motional velocity Vh of the halting-qubit atom at the 
time t m i is slowed down to the velocity Vq << Vh by the unitary decelerating 
sequence and the initial internal state |1) is also changed to |0). Here it 
is important that after the halting-qubit atom is acted on by the unitary 
decelerating sequence it is no longer acted on by next unitary decelerating 
sequences. Then the motional velocity vq of the halting-qubit atom must be 
greater than zero so that the halting-qubit atom itself can leave in the velocity 
Vq the effective spatial action zone of the unitary decelerating sequence before 
next decelerating sequence starts to apply at the time t m>i+ i = t mi + AT. The 
spatial displacement of the halting-qubit atom is given by SD = v x (AT — 
Td) during the period (AT — To) from the time t m i + To after the atom 
is acted on by the unitary decelerating sequence to the time i m ,j+i before 
next decelerating sequence is applied. The spatial displacement SD must be 
large enough to ensure that the entire wave-packet state \c 2 (t mji+ i)) of the 
halting-qubit atom at the time t m ^ + \ is outside the effective spatial action 
zone of the unitary decelerating sequence, and hence it is also much greater 
than the effective wave-packet spread of the wave-packet state \c 2 (t mt i + i)) . 
Now both the atomic motional velocity Vq and energy E = itlVq/2 are much 
less than the initial velocity Vh = \j2Eh/rrih and energy Eh, respectively. 
The wave-packet state \c 2 (t mi + T D )) is stable in the sense that the atomic 
motional energy E$ of the wave-packet state is much lower than the height of 
the intermediate potential barrier in the double-well potential field. If now 
the halting-qubit atom goes in the velocity vq through a fixed distance AR 
in the right-hand potential well it spends the time equal to AR/vo. Because 
Vh >> vq this time interval AR/vq is much longer than the time period 
AR/vh during which the atom passes the same distance AR in the initial 
velocity Vh- Thus, one may imagine that the wave-packet state of the halting- 
qubit atom is locked in the right-hand potential well for a long time (AR/v ) 
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after the atom is slowed down greatly by the unitary decelerating sequence. 
As required by the quantum program and circuit, the halting-qubit atom 
should stay in the right-hand potential well until the end time (t mr ) of the 
computational process. Then the time interval AR/v > t mr — (t m i + Td), 
where the time t m \ is the earliest one among all these m T possible times {t m i}. 
For convenience, suppose 5tf > T D . Then at the end time t mr = t + m r AT 
of the computational process the center-of-mass spatial position i?2(^m r ) °f 
the wave-packet state |c2(t mr )) of the halting-qubit atom is given by 

R2,i( t m r ) = R2(t mi + T D ) +v x (tr^ -t mi -T D ), l<i< m r , (24) 

when the halting-qubit atom enters into the right-hand potential well at the 
time t mi in the i— th cycle of the quantum circuit. Here the spatial position 
R2(t mr ) is denoted as #2,1 (W) so as to show explicitly that the position is 
dependent of the cyclic index value % = 1,2, ...,m r . Though each possible 
spatial position R2(t m i + To) is the same for the index value % — 1,2, ...,m r 
just like R2(t mi ), the wave-packet spatial position i?2,j(W) a ^ the end time 
t mr of the computational process is different for a different index value i. 
This is because the halting-qubit atom stays in the right-hand potential well 
for a longer time and hence passes a longer spatial distance before the end 
time t mr of the computational process if it enters into the right-hand poten- 
tial well at an earlier time t mi . The maximum and minimum wave-packet 
spatial positions for the halting-qubit atom at the end time t mr correspond 
to the halting-qubit atom entering into the right-side potential well in the 
first and the last cycle of the quantum circuit, respectively. Evidently, these 
m r possible different wave-packet spatial positions of the halting-qubit atom 
at the end time t mr satisfy the following inequality: 

R2,m r (tm r ) < ••• < R2,2(tm r ) < ^2,1 (W)" (25) 

Here as usual the +x coordinate direction is defined as from the left-hand po- 
tential well toward the right-hand one in the double-well potential field. The 
inequality (25) is also correct for the case 8tf < Td. Of course, in this case 
the atomic wave-packet state and its spatial position R2,i{t mr + — ^/) 
at the time t mr + To — 5tf correspond to the wave-packet state and its 
spatial position R2,i(t mr ) a t the end time t nir in the case 5tf > To, re- 
spectively. Actually, the time t mr + T D — Stf is the end time of the total 
quantum circuit consisting of the quantum circuit Q c and the unitary de- 
celerating sequence for the case Stf < T D , while if <5ty > T D the end time 
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of the total quantum circuit is really just the end time t mr of the single 
quantum circuit Q c . Hereafter only the situation 8tf > Tp is considered. 
Obviously, the halting-qubit atom at the time t mr is always in the spatial 
region [R 2 ,m r (t mr ) - 5R(t mr )/2, R 2 ,i(t mr ) + 5R(t mr )/2] of the right-hand po- 
tential well, that is, any center-of-mass spatial position R 2 ^{t mr ) of the wave 
packet state \c 2 (t mr )) for i = 1,2, ...,m r satisfies: R 2 ,i(t mr ) G [R 2 ,m r (t mr ), 
R 2 ,i(tm r )]- Here for convenience suppose that the spatial shape of the wave- 
packet state \c 2 (t mr )) is symmetrical and 5R(t mr ) is the effective spatial 
spread of the wave-packet state at the time t mr . Since the spatial dis- 
tance between any two nearest wave-packet states takes the same value: 
AR 2 ,i, i + 1 (t mr ) = R 2 ,i(t mr ) ~ R2,i+i(t mr ) = w x AT for % = 1,2, ...,m r - 1, 
as shown in Eq. (24), these m r possible wave-packet states distribute uni- 
formly in the spatial region [R 2 ,m r (t mr ) - 8R(t mr )/2, R 2 ,i{t mr ) + 5R{t mr )/2] 
of the right-hand potential well. In general, it follows from Eq. (24) that 
the spatial distance between the spatial positions R 2 ,i(t mr ) and R 2 j(t mr ) 
(1 < i < j < m r ) of the halting-qubit atom at the time t mr can be calculated 
by 

AR 2:iJ (t mr ) = R 2 ,i(t mr ) - R 2:j (t mr ) = v (j - i)AT. (26) 

Here the time difference (j — i)AT is the duration between the halting- 
qubit atom entering into the right-hand potential well in the j— th cycle 
and the i— th cycle (j > i) in the quantum circuit. Obviously, the maximum 
spatial distance which is the dimensional size of the spatial region [R 2 ,m r (t mr ), 
R2,i{t mr )} is g iven by AR 2Atrrir (t mr ) = v (m r - 1)AT. From the end time t mr 
on, there are no longer any unitary operation of the quantum circuit and any 
unitary decelerating sequence applying to the whole quantum system of the 
quantum computer. However, in order that the output state of the reversible 
and unitary halting protocol is not dependent sensitively upon any initial 
functional state in the quantum circuit the wave-packet state \c 2 (t mr )) = 
|0, CM2(t mr ), R 2t i(t mr )) of the halting-qubit atom must be changed back to 
the stable halting state such as the state |1, CMO, Rq) in a high probability. 
Thus, the halting-qubit atom must ultimately return the left-hand potential 
well from the right-hand one after the computational process finished. Here 
the control state subspace S(C) in the atomic physical system consists of 
a series of wave-packet states and is not a two-state subspace. Then the 
state \c 2 ) of the control state subspace S(C) in the quantum program Q c 
really corresponds to these wave-packet states of the halting-qubit atom in 
the right-hand potential well and also the stable halting state |1,CM0, R ) 
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finally. 

One possible scheme to force the halting-qubit atom in the right-hand 
potential well to return to the left-hand potential well is to increase the 
atomic motional energy and invert the motional direction of the atom. Now 
a unitary accelerating sequence which consists of the Raman adiabatic laser 
pulses with the proper characteristic parameter settings is constructed to 
speed up in a unitary form the halting-qubit atom in the right-hand po- 
tential well after the computational process and the unitary decelerating 
sequence finished. The characteristic parameter settings for the Raman adi- 
abatic laser pulses of the unitary accelerating sequence are clearly different 
from those of the unitary decelerating sequence. A unitary accelerating pro- 
cess could be thought of as the inverse process of a unitary decelerating 
process except the atomic motional direction and the spatial action zone. 
There are two purposes for the unitary accelerating process to speed up the 
halting-qubit atom. The first purpose is simply that after the halting-qubit 
atom is sped up by the unitary accelerating sequence it hits the right-hand 
wall of the right-hand potential well to change its motional direction and 
then returns to the left-hand potential well in a higher velocity so that the 
halting-qubit atom can pass the intermediate potential barrier to arrive in the 
left-hand potential well in a shorter time. Here define the arriving time (Tj, 
% = 1, 2, m r ; see below) as the instant of time at which the entire effective 
wave-packet state of the halting-qubit atom enters into the left-hand poten- 
tial well from the right-hand one and moreover the center-of-mass position 
of the wave-packet state is some given spatial position (e.g., Ri(toi)) within 
the left-hand potential well. For a heavy atom the wave-packet picture in 
quantum mechanics is very similar to the classical particle picture [5a]. Then 
from viewpoint of the particle picture it could be better to choose the given 
spatial position such that the mean motional speed and kinetic energy of 
the halting-qubit atom is zero at the given spatial position within the left- 
hand potential well, that is, at the given spatial position the total motional 
energy of the halting-qubit atom is pure potential energy. Evidently, when 
the halting-qubit atom arrives in the left-hand potential well there are m r 
possible different arriving times for the halting-qubit atom with the m r pos- 
sible wave-packet states {|0, CM2(t mr ), R2,i(t mr ))} (i = 1,2, ...,m r ) at the 
end time t mr of the computational process. Each arriving time corresponds 
one-to-one to a possible wave-packet state (|0, CM2(t mr ), i?2,i(W))) which 
locates at a different spatial position (i?2,j(t mr )) i n the right-hand poten- 
tial well. Generally, the first wave-packet state |0, CM2(t mr ), i?2,i(W)) wm 
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arrive in the left-hand potential well at the earliest time, while the last one 
|0, C M2(t mr ) , R2,m r (tmr)) enters into the left-hand potential well at the latest 
time. Then the second purpose is particularly important for the quantum 
control process in that the unitary accelerating sequence is really used to 
shorten greatly any time differences among the m r possible different arriving 
times for the halting-qubit atom. As pointed out earlier, all these m r possi- 
ble wave-packet states {|0, CM2(t mr ), R2,i(tm r ))} °f the halting-qubit atom 
at the end time t nw are within the spatial region [R2,m r (t mr ) — SR(t mr )/2, 
R2,i(t mr ) + 8R(tm r )/2] m the right-hand potential well. In order that any 
one of these m r possible wave-packet states can be changed back to the sta- 
ble halting state \n' h , CMO, R ) (n' h = or 1) in a high probability each of 
these m r possible wave-packet states may be acted on by the same unitary 
accelerating sequence such that any time differences among these m r possible 
arriving times can be shorten greatly. Here the effective width of the spatial 
action zone of every Raman adiabatic laser pulse of the unitary accelerat- 
ing sequence must be greater than the dimensional size of the spatial region 
[R2, mr (t mr )-SR(t mr )/2, R 2y i(t mr ) + 5R(t mr )/2]. Since the halting-qubit atom 
moves also a spatial displacement during the Raman adiabatic laser pulse the 
effective spatial-action-zone width of the Raman adiabatic laser pulse must 
also take the spatial displacement into account in addition to the dimensional 
size of the spatial region. In technique it could be better to choose those spa- 
tially uniform ultra-broadband adiabatic laser pulses [47] as the adiabatic 
laser beams of the Raman adiabatic laser pulses of the unitary accelerating 
sequence. Denote such an ultra-broadband unitary accelerating sequence as 
UA({<Pk}it + TA,t), where Ta is the total duration of the accelerating se- 
quence. When the ultra-broadband unitary accelerating sequence acts on 
the halting-qubit atom at the end time t mr of the computational process any 
one of these m r possible wave-packet states {|0, C M2(t mr ) , R 2 ,i(t mr ))} of the 
halting-qubit atom is transferred to the corresponding unstable wave-packet 
state: 

\c'2,i(t mr + T A )) = UA(W k },t mr +T A ,t mr )\0,CM2{t mr ),R2,i(tm r )) 

= \l,CM2(t mr + T A ),R2,i(tm r + T A )) 

where the internal state |0) of the halting-qubit atom is changed to the state 
|1) after the unitary accelerating sequence, meanwhile the halting-qubit atom 
is sped up from the initial motional velocity vq to a great velocity v >> vq. 
The motional velocity v usually may be greater than the motional velocity 
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Vh = ^/2E h /m h , that is, v > Vh >> v , so that the halting-qubit atom has 
an enough high motional energy to pass the intermediate potential barrier 
to enter into the left-hand potential well. Note that the motional velocity v 
of the halting-qubit atom has an upper-bound value c, where c is the light 
speed in vacuum, and usually v « c. The important point for the unitary 
accelerating process is that the center-of-mass spatial distance between any 
pair of the wave-packet states among these m r possible wave-packet states 
{|0, CM2(t mr ), R2,i(tm r ))} °f the halting-qubit atom is kept unchanged be- 
fore and after the unitary accelerating process, although the halting-qubit 
atom is accelerated under the action of the unitary accelerating sequence and 
its wave-packet spatial position has been changing along the +x direction. 
Therefore, after the unitary accelerating sequence these m r possible wave- 
packet spatial positions {R2,i{tm r +TA)} of the halting-qubit atom still satisfy 
the inequality (25) and their possible spatial distances remain also unchanged 
and are still given by Eq. (26), that is, AR 2 ,ij(t mr + Ta) = v (j — 2) AT 
for 1 < % < j < m r . After the action of the accelerating sequence the 
halting-qubit atom moves in the velocity v along the direction +x toward 
the right-hand potential wall of the double-well potential field and hits ulti- 
mately the potential wall in an elastic form. Then the halting-qubit atom is 
bounced off the right-hand potential wall and its motional direction therefore 
is reversed and hence changed to the direction — x. Evidently, this elastic 
bouncing process is unitary [5a, 48]. Now the halting-qubit atom moves in 
the velocity v (v > Vh) along the direction —x toward the left-hand poten- 
tial well. It first goes across the right-hand potential well, then passes the 
intermediate potential barrier, and finally enters into the left-hand harmonic 
potential well. Note that the spatial position R2,i(t mr + Ta) is nearest the 
right-hand potential wall among all these m r possible wave-packet spatial 
positions {R2,i{t mr +T4)}. Evidently, the shortest spatial distance between 
the spatial position #2,1 (W + Ta) and the right-hand potential wall must 
be much greater than half the wave-packet spread: 5R(t mr + Ta)/2. The 
halting-qubit atom needs to spend a short time period when the atom moves 
from the spatial position R2,i{t mr + Ta) to the right-hand potential wall, 
bounces off the potential wall, and then returns the original spatial position 
R2,i(t mr +Ta)- Denote this short period as the atomic bouncing dead time td- 
Suppose that the time period is denoted as t a when the halting-qubit atom 
arrives in the left-hand potential well from the spatial position R2,i(t mr +Ta) 
after the atom bounces off the potential wall. Then the longest time period 
of the quantum control process from the starting time (t 0i ) of the quantum 
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scattering process to the time when the halting-qubit atom arrives in the left- 
hand potential well is not longer than (t mr —t m ) + T A + td+t a . It follows from 
the inequality (25) that if the halting-qubit atom enters into the right-hand 
potential well from the left-hand one in the first cycle of the quantum cir- 
cuit, then it will first return the left-hand potential well from the right-hand 
one after the unitary decelerating and accelerating sequences, whereas the 
halting-qubit atom returns the left-hand potential well at the latest time if 
it enters into the right-hand potential well in the latest cycle of the quantum 
circuit. Suppose that the halting-qubit atom returns to the left-hand poten- 
tial well from the right-hand one and arrives at some given spatial position 
within the left-hand potential well at the arriving time T (Tj > t„ lr + T A ) 
for i — 1,2, m r if the atom enters into the right-hand potential well from 
the left-hand one in the %— th cycle of the quantum circuit. It follows from 
the inequality (25) that these m r possible arriving times {Tj, i — 1,2, ...,m r } 
satisfy the following inequality: 

T x < T 2 < ... < T mr (27) 

and the equation (26) shows that the arriving-time difference AT^j = Tj — Tj 
for 1 < i < j < m r is given by 

AT,- j = AR 2 , tJ (t mr + T A )/v = (j - i)ATv /v, (28) 

and the maximum arriving-time difference equals 

AT mril = (m r - l)ATv /v. (29) 

It is known that the time difference between the halting-qubit atom entering 
into the right-hand potential well in the j— th cycle and the i — th cycle (j > i) 
of the quantum circuit is given by (j — i)AT. But after the halting-qubit atom 
is acted on by the unitary decelerating and accelerating sequences in the 
right-hand potential well the corresponding arriving-time difference becomes 
AT,- ; j = (j — i)ATv /v. Since the motional velocity v is much greater than 
the velocity v , that is, the time-compressing factor v /v « 1, the arriving- 
time difference AT^j is much shorter than the original time difference (j — 
i)AT, indicating that the original time difference is greatly compressed after 
the time- and space-dependent quantum control process which contains the 
unitary decelerating and accelerating processes. 

During the quantum control process the halting-qubit atom carries out 
consecutively the quantum scattering process in which the atom goes from 
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the left-hand potential well to the right-hand one, decelerating process, ac- 
celerating process, elastic bouncing process, the second decelerating process 
(see below), and finally the second quantum scattering process in which the 
atom returns from the right-hand potential well to the left-hand one. But 
the internal state of the halting-qubit atom could be changed only in the 
unitary decelerating and accelerating processes among these processes. For 
simplicity, here the second decelerating process is not considered temporar- 
ily. For convenience, for the case that the halting-qubit atom enters into the 
right-hand potential well in the %— th cycle of the quantum circuit and then 
returns and arrives in the left-hand potential well at the arriving time 7] the 
wave-packet state of the halting-qubit atom in the left-hand potential well at 
the arriving time Tj is denoted as |1, CMl(Tj), i? 1; j(Tj)) for i = 1,2, ...,m r . 
Here the wave-packet state |1, CMl{Ti), Ri^(Ti)) has the same internal state 
as the state \1,C M2(t mr + Ta), R2,i(t mr + Ta)) of the halting-qubit atom 
in the right-hand potential well after the unitary accelerating process. Evi- 
dently, in an ideal case all these wave-packet states {|1, CMl(Tj), Ri^T))} 
are really the same: 

11,6^1(71), ^(71)) = \1,CM1(T 2 ),R 1>2 (T 2 )) 

= ... = \l,CMl{T mr ),R 1>mr {T mr )), 

although the arriving time T for the halting-qubit atom is different for a 
different cycle index value i = 1,2, ...,m r , as can be seen in (27). How- 
ever, the wave-packet motional state \CMl(Ti), R l! i(T i )) for i = l,2,...,m r 
is generally different from the desired motional state \CM0,Rq). Since 
the motional energy of the halting-qubit atom with the wave-packet state 
|1, CMl(Ti), Rij(Ti)) is much higher than the height of the intermediate 
potential barrier in the double-well potential field the wave-packet state 
\l,CMl(Ti), Ri :i (Ti)) is unstable and hence it must be transferred to the 
stable halting state \1,CM0,Rq). 111 general, there is a unitary transforma- 
tion that transfers completely the wave-packet state |1, CMl(Ti), i?ii(71)) 
to the halting state |1, CMO, Rq). This unitary transformation is defined as 

U({<p k }, T t + 71, 71)11, CM\{T X ), #1,1(71)) = |1, CMO, R ). (30) 

Here the duration T t of the unitary operation U({ip k },T t + 71, 71) usually 
could be much longer than the maximum arriving-time difference AT mr i = 
T mr — T\ (see Eq. (29)) and {<^} are the control parameters of the uni- 
tary operation. The unitary operation U({ip k },T t + 71,71) starts to act on 
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the wave-packet state |1, CMl(Ti), Ri^iTi)) of the halting-qubit atom in the 
left-hand potential well at the arriving time Ti and it does not change any 
internal state of the halting-qubit atom. It could be generated by applying 
the coherent Raman adiabatic laser pulse and the time- dependent potential 
field to the left-hand harmonic potential well. The coherent Raman adiabatic 
laser trigger pulse P t transfers the lower motional-energy state |1, CMO, Rq) 
to the higher motional-energy wave-packet state |1, CMl(t 0i ), i?i(ioi))> while 
here the unitary operation U({<p k }, T t + Ti, Ti) converts the higher motional- 
energy wave-packet state |1, CMl(Ti), .Ri^Ti)) into the lower motional- 
energy state |l,CM0,i? ) in the left-hand potential well. However, the 
unitary operation U({(f k },T t + Ti,Ti) is more complex than the coherent 
Raman adiabatic laser trigger pulse P t . There is a difference between these 
two motional states \CMl(T-\),R\j_(T\)) and \CM0, Rq) . The difference is in 
the atomic motional energy, momentum, spatial position, wave-packet shape 
(e.g., the effective spread) and phase, and so forth. Thus, in order to convert 
the wave-packet motional state |CMl(Ti), i2 lj:L (Ti)) into the ground motional 
state \CM0, R ) one could apply the coherent Raman adiabatic laser pulses 
and also an external time- dependent potential field to the halting-qubit atom 
in the left-hand potential well. Here the external time- dependent potential 
field may be used to modulate the left-hand harmonic potential field or even 
the whole double-well potential field. Both the external time-dependent po- 
tential field and the coherent Raman adiabatic laser pulses could also be 
thought of as the components of the state-locking pulse field. 

As shown in the inequality (27), the arriving time Tj is different for % = 
1, 2, m r for the halting-qubit atom returning and arriving in the left-hand 
potential well after the atom enters into the right-hand potential well in a 
different cycle {i— th cycle) of the quantum circuit. The earliest and latest 
arriving times are Ti and T mr , respectively. It is certain that the halting- 
qubit atom is within the left-hand potential well at the arriving time Tj 
if the atom enters into the right-hand potential well in the %— th cycle of 
the quantum circuit. However, the halting-qubit atom could not be in the 
left-hand potential well at the arriving time Tj if the atom enters into the 
right-hand potential well in the j— th cycle of the quantum circuit with the 
cyclic index j ^ % rather than in the %— th cycle. If the halting-qubit atom 
enters into the right-hand potential well in the j— th cycle (m r > j > 1) 
rather than in the first cycle of the quantum circuit, then at the arriving 
time Ti the wave-packet state of the halting-qubit atom could not be the 
state |l,CMl(Ti),i2i i i(Ti)) or the state |1, CMl(Tj), Rij(Tj)), but it could 
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be another wave-packet state |1, CMl(Ti), Rij(T\)) different from the state 
|1, CMl(Ti), -Ri,i(Ti)). Here the spatial position Rij{T{) for m r > j > 1 
is also different from Ri^(T\) or Rij(Tj) and will not be constrained to be 
within the left-hand potential well. Then the unitary operation U ({y?fc}, T t + 
Ti,Ti) of Eq. (30) transfers the wave-packet state |1, CMl(Ti), i2ij(Ti)) 
with the cyclic index m r > j > 1 to the state |1,CM0, i? ) in a probability 
less than 100%. The unitary state transformation may be generally written 
as 

|1, CMO(j), R (j)) = U({tp k }, T t + T u Ti)|l, CM\(T X ), Ri,j(Ti)) 

= Aj(CM0, Rq, T t , CMO, i? ) 
+ ^-(CM, ^(Tx), T t , T0I1, CM, i?>, (31) 

CM 

where the first term on the right-hand side is the desired state |1, CMO, Rq) 
and the second term a superposition state which is orthogonal to the de- 
sired state |1, CMO, Rq). The absolute amplitude \Aj(CM0, Rq, T t , T x )\ mea- 
sures the conversion efficiency from the state |1, CM1(7\), iZ^Ti)) (m r > 
j > 1) to the state |1,CM0, Rq) under the action of the unitary operation 
U({ip k },T t + Ti,Ti). By comparing Eq. (30) with Eq. (31) one sees that the 
amplitude A^CMO, i? ,T t ,Ti) = 1 and A^CM, i2i,i(7i), T t , 7\) = for any 
index value CM. A theoretical calculation for the amplitude Aj(CM0, R , T t , 
Ti) (m r > j > 1) usually could be more complex. 

Denote that Ho and Uo(t,to) are the Hamiltonian and time evolution 
propagator of the halting-qubit atom in the time-independent double-well po- 
tential field without the Raman adiabatic laser pulses and the time-dependent 
external potential field, respectively. The propagator U (t, t ) does not change 
any internal state of the halting-qubit atom. Then there is the relation be- 
tween both the wave-packet states |1, CMl(Ti), and |1,CM1(7}), 
R i,j( T j)) f° r J = l,2,...,m r , 

1/0(3), T0I1, C7M1(T 1 ), J R 1J -(T 1 )> = |1, CMl(Tj), Rij(Tj)). (32) 

This is because the halting-qubit atom in the wave-packet state |1, CM1(7\), 
Rij(Ti)) at the time T\ will arrive in the left-hand potential well at the ar- 
riving time Tj > Ti and moreover it is in the wave-packet state |1, CMl(Tj), 
Rij(Tj)) at the arriving time Tj. By the equation (32) and the relation 
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\l,CMl(Tj), Rij(Tj)) = |l,CAfl(Ti), the unitary state transfor- 

mation (31) may be reduced to the form 

U({cp k },T t + T 1 ,T 1 )\l,CMl(T 1 ),R 1J (T 1 )) 

= U({<p k }, T t + TuTjUoiTj, 7\) + |l, CMl{T x ), Ri^T,)). (33) 

The wave-packet motional state |CM1(T 1 ), i? 11 (T 1 )) generally is not an exact 
eigenstate of the Hamiltonian Hq of the halting-qubit atom in the double- 
well potential field, although the motional state \CM0, R ) could be approx- 
imately an eigenstate of the Hamiltonian Hq with the motional energy eigen- 
value E = Huo/2 since it is approximately the ground motional state of the 
halting-qubit atom in the left-hand harmonic potential well. Thus, equation 
(33) shows that not every state |1, CMl(Ti), i?ij(Ti)) for the cyclic index 
j = 1,2,...., m r can be converted completely into the same state |1, CMO, Rq) 
by the same unitary operation U({ip k },T t + 71, Ti). It follows from Eqs. 
(30), (31), and (33) that the amplitude Aj(CM0, R , T t , T ± ) of the state 
|1, CMO, R ) on the right-hand side of Eq. (31) can be written as 

AjiCM^Ro^Ti) = 

(1, CMliT,), R ltl (T 1 )\U (T j , T0 + I1, CM1(7\), Ri^T,)) (34) 

This equation may be used to calculate the amplitude Aj(CM0, Rq, T t , Ti) 
if the wave-packet state \l,CMl(Ti), R lt i(Ti)) and the unitary operation 
Uo(Tj,Ti) are explicitly given for j = 1,2, ...,m r . It follows from Eqs. (32) 
and (34) that the probability \ Aj(CM0, R , T t , Ti) | 2 is really the project prob- 
ability of the wave-packet state |1, CMl(Ti), Rij(Ti)) for j = 1,2, ...,m r to 
the wave-packet state |1, CMl(Ti), Ri^(Ti)). 

The quantum control process really starts at the time to« (here for con- 
venience the excitation process of the trigger pulse P t is not considered) 
and its initial state is the unstable wave-packet state |1, CMl(t j), -Ri(toi)) 
of the halting-qubit atom in the left-hand potential well. For convenience 
the mean motional velocity and kinetic energy of the halting-qubit atom 
may be prepared to be zero when the atom is prepared to be in the wave- 
packet state |l,CMl(t i), -Ri(toi)) & t the time toi. This can be achieved by 
the suitable state-dependent coherent Raman adiabatic laser trigger pulse 
P t , as shown in the quantum control unit Qh- From the viewpoint of the 
particle picture the total motional energy of the halting-qubit atom in 
the left-hand harmonic potential well is really pure potential energy at the 
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time toi due to zero motional velocity of the atom when the atom is in the 
wave-packet state |1, CMl(i i), Ri(toi))- Thus, the center-of-mass position 
Ri{toi) of the halting-qubit atom in the left-hand harmonic potential well 
is really determined uniquely by the total motional energy Eh- The wave- 
packet motional state |CMl(i j)> Ri(toi)) may be considered as a coherent 
state [49] for the halting-qubit atom in the left-hand harmonic potential 
well as it is generated from the ground motional state \CM0, Rq) by the 
coherent Raman adiabatic laser trigger pulse P t , as shown in Ref. [37]. 
The wave-packet state |1, CMl(t i), -Ri(^oi)) is transferred to the wave-packet 
state |1, CMl(Tj), i?i 5 i(Tj)) (i = 1,2, ...,m r ) by the quantum control process 
that the halting-qubit atom starts at the position Ri(t 0i ) and at the time t 0i 
from the left-hand potential well to enter into the right-hand potential well 
by the quantum scattering process, then it is manipulated by the decelerat- 
ing process, the accelerating process, and the elastic bouncing process in the 
right-hand potential well, and finally it returns to the left-hand potential well 
by another quantum scattering process and arrives in the left-hand potential 
well at the arriving time Tj. For convenience, consider the specific case that 
the halting-qubit atom arrives at the original position Ri(toi) in the left-hand 
potential well at the arriving time Tj after it goes a cycle along the double- well 
potential field in the quantum control process. Then in an ideal condition 
the wave-packet state |1, CMl(Tj), R lti (Ti)) is just the original wave-packet 
state |1, CMl(t 0i ), Ri(t 0i )) and hence it is also a coherent state. In order to 
achieve this point the total motional energy of the halting-qubit atom in the 
wave-packet state |1, CMl(Tj), R\^(Ti)) must be equal to that one Eh of the 
atom in the wave-packet state |1, CMl(t 0i ), Ri(t 0i )). Then this requires that 
the motional velocity v of the halting-qubit atom after the unitary accelerat- 
ing process be equal to the motional velocity Vh of the atom just before the 
unitary decelerating process. Actually, this means that for the halting-qubit 
atom the quantum scattering process from the left-hand potential well to the 
right-hand one is the inverse process of the second quantum scattering process 
from the right-hand potential well to the left-hand one. On the other hand, 
both the wave-packet shapes for the motional states |CMl(t j), Ri{toi)) and 
\C Ml(Ti) , i?! j(Tj)) could be different from each other in practice. The wave- 
packet spread of the motional state \CMl(Ti), Ri ti (Ti)) usually is greater 
than that of the motional state \CMl(t 0i ) : Ri(t 0i )) [5a, 48, 53] even if the 
halting-qubit atom returns the original position Ri(t 0i ) at the arriving time 
Ti in the left-hand potential well after it goes a cycle along the double-well 
potential field. This is because those processes including the quantum scat- 
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tering processes, the decelerating and accelerating processes, and the elastic 
bouncing process as well as the freely atomic motional processes in the quan- 
tum control process can change the wave-packet shape of motional state of 
the halting-qubit atom and usually could broaden the wave-packet spread of 
the motional state. In order to minimize the effect of these processes on the 
wave-packet spread of the motional state in the quantum control process the 
unitary accelerating process may be constructed as the inverse process of the 
unitary decelerating process except the atomic motional direction and spatial 
position. This means that the motional velocity Vh of the halting-qubit atom 
before the unitary decelerating process equals the velocity v of the atom after 
the unitary accelerating process. Then in the case v — v h the effect of the 
decelerating process could cancel that effect of the accelerating process on 
the wave-packet spread of the motional state and likely the effect of the first 
quantum scattering process could cancel that of the last quantum scatter- 
ing process in the quantum control process. As a result, the net effect for 
both the decelerating and accelerating processes and both the first and the 
last quantum scattering processes on the wave-packet spread of the motional 
state could become so small that it can be neglected. Thus, in the quantum 
control process the effect on the wave-packet spread of the motional state 
could mainly come from the elastic bouncing process and the freely atomic 
motional processes. However, this effect could also be small and negligible if 
the halting-qubit atom has a large mass m h and there is a short period for the 
quantum control process [53]. Therefore, in the case v = v h the wave-packet 
motional state |CMl(Tj), Ri^(Ti)) (i = 1,2, ...,m r ) still could be considered 
approximately as the original motional state |CMl(t j), -Ri(^oi)) in practice, 
although the wave-packet spread of a motional state of the halting-qubit 
atom could change slightly after the quantum control process. Since the mo- 
tional state \C Ml(toi) , Ri(toi)) (i = 1) is a coherent state the wave-packet 
motional state |CMl(Ti), R^i(Ti)) may be approximately expressed in a 
coherent-state form [49, 37, 5c], 




where the relevant global phase factor is omitted and the state \k) is the en- 
ergy eigenstate of the Hamiltonian of the halting-qubit atom in the left-hand 
harmonic potential well and a a complex parameter. The absolute value \a\ 2 
is the mean motional energy of the halting-qubit atom in the wave-packet 



82 



motional state \C Ml(Ti) , Ri ^Tx)) in the left-hand harmonic potential well 
[5c]. Now consider other wave-packet motional states {\CMl(Ti), Rij(Ti))} 
(m r > j > 1). If the entire effective wave-packet motional state |CMl(Ti), 
.Rij(Ti)) is within the left-hand harmonic potential well at the time Ti, it is 
also a coherent state approximately just like the coherent state |CMl(Ti), 
or \CMl{Tj),R lt j{Tj)). Here these two coherent states |CMl(Ti), 
Rij(Ti)) and \CMl(Ti), Ri^Ti)) are connected by the unitary transforma- 
tion (32). If now the right-hand potential wall of the left-hand potential 
well is sufficiently high, then that the time evolution process of the halting- 
qubit atom in the left-hand potential well is governed by the Hamiltonian 
of the double-well potential field is really reduced to the simple one that 
the time evolution process is governed by the Hamiltonian of the single left- 
hand harmonic potential well. Here the Hamiltonian Ho and the propagator 
Uo(t,to) = exp[— iH (t — to)/H] of the double-well potential field are also re- 
duced to those of the single left-hand harmonic potential field, respectively. 
Since the state \k) in the coherent state of Eq. (35) is an eigenstate of the 
Hamiltonian of the single left-hand harmonic potential well, there holds the 
eigen-equation: Uq(Tj, Ti) + \k) = exp[i(k + l/2)uj (Tj —Ti)]\k). It follows 
from Eqs. (32) and (35) that the coherent state |CMl(Ti), Ri,^)) for 
j = 1, 2, m r can be written as 



|CMl(Ti), J Ri, J -(Ti)) = \aeM^o(Tj - Ti)]) 
= exp(-^| 2 )f; {aeXpl ^- Tl)]}fc |A:) (35a) 



2 to 

where the global phase factor is also omitted and c^o — a/ K/rrih is the basic 
oscillatory frequency of the harmonic oscillator, i.e., the halting-qubit atom in 
the left-hand harmonic potential well, and K a force constant of the harmonic 
oscillator. 

There are two different situations to be considered below. For the first 
situation the maximum arriving-time difference T mr — 7\ for the halting- 
qubit atom is much shorter than the basic oscillatory period 2ti/uoq of the 
halting-qubit atom in the left-hand harmonic potential well, that is, T mr — 
Ti << tv/uq. In the first situation there is not the second unitary decelerating 
process in the quantum control process. From the viewpoint of the particle 
picture the halting-qubit atom at the position i?i(toi) leaving the left-hand 
potential well needs to spend roughly the time ti/uq. Then this really means 
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that in the case T nir — 7\ << ti/ujq each of these m r possible wave-packet 
states {\l,CMl(Ti), Ri t j(Ti))} for the halting-qubit atom is able to enter 
completely into the left-hand potential well at the time T\. Once each of 
these m r possible wave-packet states {|1, CMl(Ti), Rij(Ti))} of the halting- 
qubit atom enters completely into the left-hand potential well at the time T±, 
one may suddenly change the double-well potential field to the single left- 
hand harmonic potential field. Then in this case it is much easier to calculate 
theoretically the amplitude Aj(CM0, Ro, T t , T\) in Eq. (31) as the quantum 
dynamics can be exactly solved in a single harmonic oscillator even when the 
Hamiltonian of the harmonic oscillator is time-dependent. Now all these m r 
possible wave-packet states {|1, CMl(Ti), R 1 j(T i ))} are the coherent states 
which are given explicitly by Eqs. (35) and (35a). With the help of Eq. (32) 
and by inserting the coherent states of Eqs. (35) and (35a) into Eq. (34) 
the absolute amplitude \Aj(CM0, Ro, T t , 7\)| of Eq. (34) is reduced to the 
simple form 

^■(CMO, #0,^,71)1 = \(1,CM1(T 1 ),R 1A (T 1 )\1,CM1(T 1 ),R 1J (T 1 ))\ 

= exp{-|a| 2 (l - cos[u (Tj - 7\)])} (36) 

where the relation for a pair of coherent states \a) and \/3) : (cx\f3) = 
exp{-\a\ 2 /2-\(3\ 2 /2+a*(3} [49, 51b, 52] is used. When the arriving time Tj = 
7i it follows from Eq. (36) that the absolute amplitude (A^CMO, R , T t , T x )\ 
= 1. This result is consistent with Eq. (30). However, if the arriving time 
Tj > T\ (m r > j > 1), then equation (36) shows that the absolute ampli- 
tude \Aj(CM0, Ro, T t , Ti) \ decays exponentially as the mean motional energy 
| a | 2 of the halting-qubit atom and the term (1 — cos[uJo(Tj — ^i)]). In or- 
der that the absolute amplitude \Aj(CM0, R ,T t ,Ti)\ is close to unity one 
should decrease greatly either the mean motional energy \a\ 2 and the basic 
oscillatory frequency ujq or the arriving-time differences {(Tj — T\)}. Since 
for the current situation the basic oscillatory frequency u>o is fixed and the 
mean motional energy \a\ 2 of the halting-qubit atom must be much greater 
than the height of the intermediate potential barrier one can only decrease 
greatly the arriving-time differences {(Tj — 7\)} to make the absolute am- 
plitude \Aj(CM0, Ro, T t , Ti)| close to unity. Because the arriving-time dif- 
ference Tj -Ti = (j - i)ATv /v for 1 < % < j < TTlf . clS shown in Eq. (28), 
and the atomic motional velocity v — Vh is determined by the mean motional 
energy \a\ 2 in the current situation, one has to slow down greatly the atomic 
motional velocity v h to a much smaller velocity v which could be close to 
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zero by the unitary decelerating process, so that the arriving-time differences 
can be shortened greatly and thus the probability \Aj(CM0, Rq, T t , T\) | 2 of 
Eq. (36) can be close to 100%. Though a larger atomic motional veloc- 
ity v leads to a smaller arriving-time difference (Tj — Ti), it also leads to a 
larger mean motional energy | ct | 2 cc v 2 , and hence the absolute amplitude 
\Aj(CM0, R , T t) Ti)| does not become larger and closer to unity significantly 
for a larger atomic motional velocity v — v^. 

Now consider the second situation that the arriving-time differences are 
not always much shorter than the basic oscillatory period 2tt/lj of the 
halting-qubit atom in the left-hand harmonic potential well, that is, the 
condition Tj — Ti < ti/uj may not hold for some index values j = 1, 2, m r . 
In this case not every of these m r possible wave-packet states {|1, CMl(Ti), 
.Rij(Ti))} for the halting-qubit atom is able to enter completely into the left- 
hand potential well at the time T\. Since in this case some arriving-time dif- 
ferences {(Tj — Ti)} are larger the corresponding amplitudes {\Aj(CM0, R , 
T u T\) | } of Eq. (36) are not close to unity. One scheme to solve the problem 
is that the basic oscillatory frequency Uq is switched to a smaller value such 
that the new basic oscillatory period can be much longer than the maximum 
arriving-time difference and meanwhile the mean motional energy \a\ 2 of the 
halting-qubit atom is lowed down accordingly before the halting-qubit atom 
enters into the left-hand harmonic potential well. In this case the amplitudes 
{\Aj(CM0, R , T t , Ti)|} become larger and may be close to unity, as can be 
seen from Eq. (36). This scheme may be described as follows. Since in 
the current situation the shortest arriving-time difference ATj+i^ = ATvq/v 
(1 < i < m r — 1) generally is long and it is also much longer than the pulse 
duration of the Raman adiabatic laser pulse, one may use a unitary deceler- 
ating sequence consisting of the Raman adiabatic laser pulses to decelerate 
the halting-qubit atom so that the atom has a low motional energy and speed 
before the atom enters into the left-hand potential well from the right-hand 
one by the second quantum scattering process. Here every arriving-time 
difference AT^ for 1 < i < j < m T must be kept unchanged before and 
after the unitary decelerating process. This decelerating process is the sec- 
ond unitary decelerating process of the quantum control process. It is really 
a space-compressing process for the m r possible wave-packet states of the 
halting-qubit atom. The unitary decelerating sequence may be applied after 
the atom bounces off the right-hand potential wall of the double-well poten- 
tial field and its spatial action zone could be the same as that one of the 
first unitary decelerating sequence of the quantum control process. After the 
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quantum computational process finished and long before the halting-qubit 
atom returns the left-hand potential well the original oscillatory frequency 
ujq of the halting-qubit atom in the left-hand potential well is switched to 
a much smaller oscillatory frequency and the height of the intermediate po- 
tential barrier is changed to be much lower than the motional energy of the 
halting-qubit atom such that the atom with a low motional energy can also 
enter completely into the left-hand potential well even after it is decelerated 
by the second unitary decelerating sequence. Of course, the intermediate 
potential barrier may also be cancelled (i.e., switched off) and in the case 
the height of the right-hand potential wall of the left-hand potential well 
has to be lowed down correspondingly. These operations can be achieved 
by applying external potential fields to modulate the left-hand harmonic po- 
tential well and the intermediate potential barrier. Denote that u c is the 
basic oscillatory frequency for the halting-qubit atom in the new left-hand 
harmonic potential well and |a c | 2 the mean motional energy of the halting- 
qubit atom after the atom is decelerated by the second unitary decelerating 
sequence. Then u c « u and |o; c | 2 << \a\ 2 . Now the basic oscillatory pe- 
riod 2it/uj c is longer for a smaller basic oscillatory frequency uj c . Since the 
time-compressing factor vq/v = Vo/vh << 1 and 2tt/u « 27i/u c now the 
condition: T mr — T\ — (m r — l)AT(f /f/ l ) << n/u c may be met easily by 
setting an enough small basic oscillatory frequency u c . Then in this case 
each of these m r possible wave-packet states {|1, CM1(7\), Rij(Ti))} is able 
to enter completely into the new left-hand potential well at the time Ti af- 
ter the second decelerating process. After each of these m r possible states 
{|1, CMl(Ti), Rij(Ti))} enters completely into the new left-hand potential 
well one may change suddenly the double- well potential field to the single left- 
hand harmonic potential field by increasing greatly and quickly the height of 
the right-hand potential wall of the left-hand potential well. Here the spatial 
position Rl ^Ti) of the wave-packet state |1, CMl(Ti), i^ ^Ti)) in the new 
left-hand potential well is determined by the mean motional energy \a c \ 2 in- 
stead of \a\ 2 . Obviously, the wave-packet state \a c ) = |1, CMl(Ti), R\^(T{)) 
of the new left-hand potential well is fully different from the original state 
| a) = |l,CMl(Ti),i2i ) i(7 , i)) and the state |1, CMl(t 0i ), Ri{hi)) of the origi- 
nal left-hand potential well, but both the mean motional velocity and kinetic 
energy are zero for the halting-qubit atom with any one of the three wave- 
packet states. Evidently, these wave-packet states {|1, CMl(Ti), R^j(Ti))} 
are still approximately coherent states of the new left-hand harmonic po- 
tential well and the equation (36) holds also for them, where the oscillatory 
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frequency uj and the mean motional energy \a\ 2 are replaced with u c and 
\a c \ 2 , respectively. Now the time-compressing factor Vo/vh << 1 leads to that 
w c {Tj — Ti) « 1 and |a c | 2 (l — cos[u c (Tj — Ti)]) << 1 for j = 1, 2, ...,m r . 
Thus, it follows from Eq. (36) that the project probability of the wave-packet 
state ll^Ml^i),^.^)) to the wave-packet state |1, CMl(Ti), i^ ;1 (Ti)) 
for the halting-qubit atom in the new left-hand potential well may be given 
by 

= exp{-2|a c | 2 (l - co8[u c (Tj - Ti)])}. (37) 

Here one must pay attention to that unlike in Eq. (36) the project probability 
in Eq. (37) generally is not equal to the probability \Aj(CM0, Rq, T t , T\)\ 2 
of the desired state \1,CM0,Rq) in Eq. (31), because the current left-hand 
potential well is not the original one. 

In the definition (30) of the unitary transformation U({(f k },T t + Ti,Ti) 
the wave-packet state |1, CM1(7\), Ri^(Ti)) now is replaced with the state 
|1, CMl(Ti), R{ i{T\j) of the halting-qubit atom in the new left-hand po- 
tential well, while the state \1,CM0, Rq) is still of the atom in the original 
left-hand potential well. If now there is a unitary transformation Vi({ip k }, r+ 
7i,7i) such that the motional state \a c ) = |CM1(T 1 ), R^^Ti)) of the new 
left-hand potential well is completely converted into another coherent state 
\(3) of the original left-hand potential well and there is also another unitary 
transformation ^({v^fc}; T t +Ti, r+Ti) such that the coherent state \/3) is fur- 
ther converted completely into the desired motional state \CM0, Rq), then 
it can turn out that the project probability of Eq. (37) is really just the 
probability \Aj(CM0, R , T t , T^] 2 of the desired state \l,CM0,R o ) in Eq. 
(31), and the probability \Aj(CM0, R , T t , T{)\ 2 now can be expanded as the 
power series of the arriving-time difference (Tj — Ti), up to the second order 
approximation, 

\Aj(CM0, Ro^TJl 2 = 1 - |a c | 2 (^ c ) 2 (AT) 2 (j - l) 2 (^oK) 2 + .... (38) 

Therefore, the motional state |CM1(T 1 ), i?£ -(Ti)) of the new left-hand po- 
tential well for j = l,2,...,m r is converted into the desired motional state 
\CM0,R ) in the probability \Aj(CM0, R , T t , T x ) | 2 of Eq. (38) by the uni- 
tary transformation U ({<^,t}, T t + Ti, Ti) of Eq. (30) which now is a product 
of the unitary transformations V\ and V2. The equation (38) shows that the 
lower bound of the probability \Aj(CM0, Rq,^,^)] 2 of the desired state 
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|1,CM0, R ) in Eq. (31) is dependent upon the time-compressing factor 
Vo/vh in a quadratic form when the time-compressing factor v^/vh « 1. 
One sees from Eq. (38) that the probability \Aj(CM0, Ro,T t ,Ti)\ 2 becomes 
closer to unity when the initial motional velocity (vt = v) increases. This 
point is quite different from that one in the first situation (see Eq. (36)) and 
it may make the scheme to increase the probability \Aj(CM0, R ,T t ,Ti)\ 2 
in the second situation better than that one in the first situation. 

Obviously, the unitary transformation U({ip k },T t + 71,71) of Eq. (30) 
now may be explicitly expressed as 



U(W k }, T t + 71, 71) = V 2 (W k }, T t + T 1 ,r + TJV^W,}, r + 71, 71). 



Here the unitary transformations Vi({ip k }, r+Ti, 71) and V2({<p k }, T t +71, r+ 
71) are respectively defined as 



The unitary transformation V 2 ({(p k }, T t + 71, r + 71) can be generated sim- 
ply by the coherent Raman adiabatic laser pulse applying to the original 
left-hand harmonic potential well [37] as both the motional states \/3) and 
\CM0, Ro) are coherent states of the halting-qubit atom in the original left- 
hand potential well. Of course, the unitary operation V2({<p k }, T t +71, r + 71) 
is the unity operation if the motional state \f3) is just \CM0, Rq). Now one 
needs to construct the unitary transformation Vi({<p k },T + 71,71). Firstly, 
one must pay attention to that both the desired motional state \CM0, Rq) 
and the coherent motional state \/3) belong to the halting-qubit atom in the 
original left-hand harmonic potential well whose basic oscillatory frequency 
is lu , while the coherent motional state \a c ) = |CM1(71), 72^(71)) is of the 
atom in the new left-hand harmonic potential well whose basic oscillatory 
frequency is u c . Therefore, the unitary transformation Vi({ip k }, r + 71, 71) is 
involved in the time-dependent oscillatory-frequency-varying evolution pro- 
cess of the halting-qubit atom in the left-hand potential well and in the 
process the initial and final oscillatory frequencies are u c and uq, respec- 
tively. In order to construct the unitary transformation Vi({(p k }, r + 71, 71) 
one may apply the time- dependent potential field to modulate the left-hand 
harmonic potential well with the initial and final basic oscillatory frequen- 
cies ijj c and uq, respectively, after the halting-qubit atom enters into the new 



V 1 ({i Pk },T + T 1 ,T 1 )\a c ) = \{3), 



(39) 



V 2 ({<p k }, T t + T 1 ,r + 71)|/3) = \CM0, Ro). 



(40) 
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left-hand harmonic potential well. Note that the time-dependent potential 
field does not change any internal state of the halting-qubit atom in the 
left-hand potential well. Denote that H(t) and U(t,to) are the Hamiltonian 
and evolution propagator of the halting-qubit atom in the time-dependent 
frequency-varying left-hand harmonic potential well, respectively. Then the 
time evolution propagator U (t, t ) of the atomic physical system which can 
be considered as a conventional harmonic oscillator may be written as 

t 

[/(Mo) = Texpi-'-J H(t')dt'} (41) 

to 

where the operator T is the Dyson time-ordering operator and the time- 
dependent frequency- modulation Hamiltonian H(t) of the harmonic oscilla- 
tor is simply written as 

where the frequency- modulation function ui{t) satisfies ou(t ) = u c for t = T\ 
and u){t) = ujq for t = t + T 1 . The unitary transformation Vi({<p k }, r + Ti, Ti) 
is just the propagator U(t + Ti,Ti). The unitary propagator U(t + Ti,Ti) 
may be generally expressed as [50g, 50h, 50i, 51, 52] 

Uir + T,, TO = exp[^(r,T 1 )a 2 -^(r,T 1 )*(a+) 2 ] 

xexp[-i^(r,7i)a + a] (43) 

where the operators a and a + are creation and destruction operators of 
the conventional harmonic oscillator, respectively, the complex parameter 
z{t,Ti) = \z(t, Tl)| exp[i(p(r, Ti)], and = 0(r,Ti) is the real angular 
frequency. The time-dependent frequency-modulation function u(t) of the 
Hamiltonian (42) must be designed suitably. The unitary operation U(t + 
Ti,Ti) may be generally determined from Eq. (41) and the Hamiltonian of 
Eq. (42) by the Magnus expansion method [50a, 50c, 50d, 50e] or by the Lie 
algebra approach [50b, 50c, 50f, 50g, 50h, 50j]. The most convenient method 
to calculate approximately the parameters |z(t, Ti)|, </?(r,Ti), and 0(r,Ti) 
in Eq. (43) from the Hamiltonian H(t) of Eq. (42) could be the Magnus 
expansion method [50a, 50c] or the average Hamiltonian theory [50d, 50e, 
50j]. 
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The unitary operation Vi({ip k }, r + Ti, Ti) = {7(r + Ti, Ti) of Eq. (43) can 
be determined from the Hamiltonian of Eq. (42) with the time-dependent 
frequency-modulation function u(t), but here the important thing is how 
to generate a proper time-dependent frequency- modulation function uj(t) so 
that the coherent state \a c ) = |CMl(Ti), R^ ^Ti)) can be completely trans- 
ferred to another coherent state \f3) by the unitary operation Vi({ip k },T + 
Ti,Ti) in Eq. (39). Firstly, the second term exp[— i<f>(r, Ti)a + a] of the uni- 
tary operation U(r + Ti,Ti) can transfer the coherent state \a c ) to another 
coherent state [49, 5c, 51, 52], 

exp[-i0(r,Ti)a + a]|a c ) = |o; c exp[-^(r, Ti)]). (44) 

Note that the first term of the unitary operation U(t + Ti, Ti) is a double- 
photon unitary operator S(z) = exp[|(za 2 — z*(a + ) 2 )} of the conventional 
harmonic oscillator [51, 52]. Then, the double-photon unitary operator 
S(z(t, Ti)) acting on the coherent state of Eq. (44) will transfer the coherent 
state to other coherent states. There is a general formula to calculate the 
transition amplitude between a pair of coherent states |a) and \j3) induced 
by the double-photon unitary operator S(z) [51]: 

MSWfa) = C^ 2 eM-^(\ ai \ 2 + \^\ 2 ) + 0^1(3,} 

x exp{^ r C- 1 [/3?exp(^) - (a{) 2 exp(-^)]} (45) 

where the parameters z = rexp(i^), C r = coshr, and S r = sinhr. Now in- 
serting the unitary operation Vi({</? fc }, r + Ti, Ti) — U(t + Ti, Ti) ofEq. (43) 
into the equation (39), then using the state transformation (44) and the for- 
mula (45), and setting (3 1 = a c exp[— i(j>(r, Ti)] = \a c \ exp(i7 c ) exp[— i<f>(r, Ti)], 
ai = P = \{3\exp(ij), r = \z(t,Ti)\, and (p = if(r, Ti), one obtains 

\{P\Vi({<p k },T + T 1 ,T 1 )\ ae )\ = 
C; 1 ' 2 ex p{-i(|a c | 2 + |/?| 2 ) + C^\a c (3\ cos[- 7 + l c - <Kt, Ti)]} 
x exp^C^M 2 cos[2 7c - 20(r, Ti) + (p( T , Ti)]} 

xexpi-^C-^^cos^ + ^r,^)]}. (46) 
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In order that the coherent state \a c ) = |CMl(Ti), R\ 1 (T' 1 )) can be trans- 
ferred completely to the coherent state \f3) by the unitary operation Vi({(p k }, 
t + Ti,Ti) the absolute amplitude \((3\Vi({<p k },T + Ti,Ti)\a c )\ must be equal 
to unity. Then it follows from Eq. (46) that there holds the relation when 
((/^({y^r + T^TOK)! 2 ^: 

lnC r = -M 2 (l - S r C; 1 cos[2 7c - 20(r,T 1 ) + ¥>(r,Ti)]) 

+2\a c p\C; 1 cos[- 7 + 7 C - <t>ij, Ti)] 

- |/3| 2 (1 + ^.C- 1 cos[2 7 + <p(r, Ti)]). (47) 

Equation (47) is used to construct the time-dependent frequency-modulation 
function uo(t) in the Hamiltonian H(t) of Eq. (42), because the equation (47) 
has to be met if there exists the Hamiltonian H(t) of Eq. (42) such that 
the coherent state \a c ) can be completely converted into the coherent state 
\f3) by the unitary transformation U(t + Ti,Ti). Here the time-dependent 
frequency- modulation function u(t) must satisfy the initial and final condi- 
tions: oj(T\) = uj c and u>(t + Ti) = cu . The amplitude \a c \ and phase 7 C are 
given by the coherent state \a c ) and hence there are only four independent 
variables \(3\, 7, and r = \z(t,Ti)\, 0(r, Ti) in Eq. (47). In particular, when 
r = such that C r = 1, lnC r = 0, and S r = 0, the equation (47) is reduced 
to the form 

(H - \P\) 2 + 2\a c (3\(l - cos[- 7 + 7 C - 0(r, Ti)]) = 0. 

This equation has the unique solution: \(3\ = \ct c \ and 7 = 2/c7r + 7 c — 0(r, Ti), 
where k is an integer and usually set to zero. This result is obvious because 
the unitary transformation Vi({</? fc }, t+Ti, Ti) becomes the unitary operation 
exp[— i<f>(r, Ti)a + a] if r = 0. For a general case r = |z(t, Ti)| 7^ the 
solutions to the equation (47) are more complex. For simplicity, suppose 
that the phase 7 = const. The equation (47) is reduced to the form 

ai|/3| 2 + 6i|/3|+ci = (48) 

where the three coefficients a±, bi, and c\ are given by 

ai = 1 + SrC' 1 cos[2 7 + (p( T , Ti)], 

61 = -2|a c |C T T 1 cos[-7 + 7 C - 0(r, Ti)], 
ci = In C r + |a c | 2 (l - cos[2 7c - 2^(r,Ti) + ^(r,7i)]). 
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Here the parameters r = \z(t, ip(r, Ti), and 0(r,Ti) are determined 
from the Hamiltonian //(t) of Eq. (42). The equation (48) has the solutions 
when the coefficient a\ ^ : 

101 = 2^{-6i ± V^ 7 ^}- 

The equation (48) has a real solution at least only when b\ — 4a\Ci > 0. 
However, the reasonable solutions to the equation (48) are those positive 
real solutions. Therefore, the time-dependent frequency-modulation function 
uj(t) of the Hamiltonian H(t) of Eq. (42) must be constructed such that there 
exists at least one positive real solution to the equation (48). 

After finishing both the quantum computational and quantum control 
processes the initial state |0, CMO, -Ro)|0)|/ r (^o)) with the initial functional 
state \f r {xo)) is transferred by the quantum circuit Q c to the output state 
|1, CMO, i? o )|l)|0) in the probability |A(CM0, R , T t , T\)| 2 given by Eq. (36) 
or (38), where the indices i and xo satisfy xq = (xf — i + l)modm r or 
i = (xf — xo + l)modm r for i = 1,2, ...,m r . Obviously, the halting state 
|l,CMO,-Ro) and the branch-control state \bh) = |1) in the output state can 
also be further converted completely into the halting state |0, CMO, Rq) and 
the state |0) by the unitary transformations, respectively. Now the quan- 
tum circuit Q c with the initial state \0, C MO, R )\0)\f r (x )) outputs the de- 
sired state |0, CMO, i? )|0)|0) in the probability \Ai(CM0,R , T t ,T{)\ 2 given 
by Eq. (36) or (38). As shown in Eqs. (36) and (38), the probability 
\Ai{CM0, Rq, T t , Ti)| 2 for i = (xj — xo + l)modm r = 1,2, ...,m r is close 
to 100% when the time-compressing factor Vo/vh « 1 no matter what the 
initial functional state |/ r (^o)) is in the initial state \0,CM0, Ro)\0)\f r (x )) 
of the quantum circuit. This result shows that the desired output state 
(or the state of Eq. (31)) of the quantum circuit Q c is almost independent 
upon any initial functional state l/r^o)) with x = 0, 1, ...,m r — 1. This is 
just the characteristic feature of the reversible and unitary halting protocol 
based on the state-locking pulse field. Evidently, such a quantum circuit 
whose output state is dependent insensitively upon any initial state could 
be used to build up an efficient quantum search process just like the quan- 
tum program and circuit proposed in the previous paper [24]. However, the 
lower bound for the probability \A, t (CM0, R , T t , T\)| 2 of Eq. (36) or (38) 
for any index value i = l,2,...,m r must be greater than some threshold 
value if the quantum circuit Q c is used to solve efficiently the unsorted quan- 
tum search problem in the Hilbert space of an n— qubit quantum system. 
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Here the Hilbert space should be the quantum search space. Note that the 
2™— dimensional Hilbert state space S(Z 2 n) of an n— qubit quantum system 
is really a direct product state space of n two-dimensional additive- cycle- 
group state subspaces S(Z 2 ) : S(Z 2n ) = S(Z 2 ) (g) S(Z 2 ) (g) ... (g) S(Z 2 ). It 
can turn out that the lower-bound probability of Eq. (36) or (38) must 
be greater than (1 — In p(n)/n), here p(n) is a polynomial of the qubit 
number n, if the minimum transfer probability from the entire initial state 

n n 

0{|O,CMO, J R o )|O)|/ r (a;o))}fc to the final state (g){|0, CMO, Rq) |0> |0)} fc is 

k=l k=l 

equal to l/p{n) in the 2 n — dimensional quantum search space S(Z 2 n). Thus, 
if the time-compressing factor v /v is small enough, that is, v /v « 1, such 
that the lower-bound probability of {{A^CMO, R , T t , T\)| 2 } of Eq. (36) or 
(38) is greater than (1 — lnp(n)/n), then the quantum search process based 
on the quantum circuit Q c is enough powerful to solve efficiently the unsorted 
quantum search problem in the Hilbert space of the n— qubit quantum sys- 
tem. Here the quantum parallel principle may be applied in the way that 
quantum computation is performed simultaneously in all factor state sub- 
spaces of the whole Hilbert space, and in this case the quantum parallel 
principle is compatible with both the unitary quantum dynamical principle 
and the mathematical logic principles used in the quantum search process. 

The amplitudes {A^CMO, R , T t , T\)} of the state \1,CM0,R ) in Eq. 
(31) corresponding to different initial states in the quantum circuit Q c are 
dependent upon the time-compressing factor Vq/v and the atomic mean mo- 
tional energy |a c | 2 , as shown in Eqs. (36) and (38). Of course, each of these 
amplitudes is also dependent on a global phase factor which may be differ- 
ent for a different index i. Obviously, the global phase factor has not a net 
contribution to the probability \Ai(CM0, R , T t , 7i)| 2 . Therefore, the basic 
principle to make the reversible and unitary halting protocol insensitive to 
any initial state includes two parts, one is by the time-compressing process to 
decrease the difference among the dynamical phase factors which determine 
the probabilities {|A(CM0, #0, T t , T^ 2 } of the state \1,CM0,R ) in Eq. 
(31), another is by manipulating the Hamiltonian that governs the quantum 
control process to decrease the difference among the dynamical phase fac- 
tors. Take a simple and intuitive physical model to illustrate the principle. 
Consider the simple physical model: a simple rotation pulse exp(— iu p I x At) 
acting upon the state |0) to convert it into the state |1) : 

exp(— iu p I x At)\0) = cos^jO) — «sin#|l). 
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Here the pulse field strength uj p is the control parameter of the Hamiltonian 
H p = uj p I x = u p a x /2 (a x is the Pauli operator) that governs the state transfer 
process from the state |0) to the state |1). The Hamiltonian H p is time- and 
space-independent. Note that the eigenvalue of the Hamiltonian H p is u p m x , 
here m x = 1/2 or —1/2 is the eigenvalue (h = 1) of the operator I x . Then the 
dynamical phase factor exp[— im x u p At] is determined by the dynamical phase 
angle, i.e., the rotation pulse angle 9 = 9(At) = uj p At/2. The dynamical 
phase angle is proportional to the pulse duration At and the pulse field 
strength oj p . The amplitude sin6> of the state |1) and the transfer efficiency 
(sin 2 6>) from the state |0) to the state |1) induced by the rotation pulse are 
determined by the dynamical phase angle 9 up to a global phase factor —i. If 
the rotation pulse angle 9 = u p At/2 = tt/2, then the state |0) is transferred 
completely to the state 1 1) by the rotation pulse and in this case the pulse 
duration At is equal to At p = ty/uj p . Now suppose that the pulse duration At 
is slightly different from At p , for example, At = At p — St with \St\ « At p . 
Then there is a difference between the dynamical phase angles 9 (At) and 
9{At p ) = tt/2, which is given by A9(At) = 9(At) - 9{At p ) = -u p St/2. 
Therefore, the transfer efficiency is given by 

sin 2 #(At) = cos 2 [A#(At)] = 1 - [u p St/2] 2 + .... 

One sees that the transfer efficiency sin 2 9 is determined by the dynamical- 
phase-angle difference A#(At). Now consider the case \St\ < n/u p such 
that |A#(At)| < 7r/2. Then the smaller the dynamical-phase-angle difference 
|A#(At)|, the higher the transfer efficiency sin 2 #. Note that the dynamical- 
phase-angle difference |A#(At)| is proportional to the time difference St and 
the pulse field strength lo p which controls the Hamiltonian H p . Thus, the 
transfer efficiency may be increased and is closer to unity by decreasing ei- 
ther the time difference or the pulse field strength or both. By compar- 
ing this formula with Eq. (38) one sees that there is a correspondence 
fit <-> (j — \)ATvq/v. Now suppose that this simple rotation-pulse model 
is available approximately for the quantum control process of Eq. (2). Then 
the probability 1 — \e(t, t j)\ 2 of the state |c 2 ) of the control state subspace 
in Eq. (2) could be evaluated roughly by setting St = (j — l)ATv /v for 
j = 1, 2, m r in the transfer efficiency sin 2 9 (At) above, 

1 - \e(t,t 0j )\ 2 w 1 - l -[u p (j - l)AT(v /v)} 2 . 

The probability of the state \c2) is dependent on the time-compressing factor 
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Vq/v in a quadratic form when vq/v « 1. Actually, the current atomic phys- 
ical model is much more complicated than the simple rotation-pulse model 
and especially the Hamiltonian that governs the quantum control process is 
time- and space-dependent in the atomic physical model so that this simple 
rotation-pulse model is not suited for describing the quantum control process 
in the atomic physical model, but the basic principle to make the reversible 
and unitary halting protocol insensitive to the initial states is the same for 
both the models! 

6. Discussion 

In the paper the reversible and unitary halting protocol and the state- 
locking pulse field have been investigated in detail and extensively. Though 
the halting protocol for the universal quantum computational models usu- 
ally is irreversible and non-unitary, it can be made reversible and unitary 
in the specific case that any initial state is limited to a single basis state in 
the halting protocol, and in this case it can also be simulated efficiently and 
faithfully (in a probability close to 100%) within the universal quantum com- 
putational models. From the viewpoint of the conventional halting-operation 
property that the output state of the halting operation may be completely 
independent of any initial state the reversible and unitary halting protocol 
is generally different from the classical irreversible one. In quantum com- 
putation the reversible and unitary halting protocol can only achieve such a 
unitary halting-operation property that the output state of the reversible and 
unitary halting protocol may be almost independent of any initial state (in a 
probability approaching to 100%) but it can not be completely independent of 
any initial state due to the limitation of unitarity, while the classical halting 
protocol may have completely the conventional halting-operation property 
in classical computation. The unitary halting-operation property is very im- 
portant for the reversible and unitary halting protocol because whether or 
not the quantum search process built up out of the reversible and unitary 
halting protocol can solve efficiently the quantum search problem is mainly 
dependent upon this property. The state-locking pulse field plays a key role 
in constructing such a reversible and unitary halting protocol that has the 
unitary halting-operation property in quantum computation. It is shown in 
the paper that a reversible and unitary halting protocol may be simulated 
efficiently in a real quantum physical system. A simple atomic physical sys- 
tem which consists of an atomic ion or a neutral atom in the double-well 
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potential field therefore is proposed to implement the reversible and unitary 
halting protocol. The correctness for the atomic physical model to simulate 
efficiently the reversible and unitary halting protocol is based on these facts 
that (i) for a heavy atom the atomic wave-packet picture quantum mechani- 
cally is very close to the classical particle picture, (ii) the quantum motional 
behavior of an atom is much like the classical motional behavior of a particle 
as an atom is generally much heavier than an electron, and (Hi) the time 
evolution processes for the atomic physical system in a variety of atomic 
motions still obey the Schrodinger equation and hence are governed by the 
unitary quantum dynamics. Perhaps this simple physical model could not be 
the best one, but it does provide one with a much intuitive physical picture 
to understand the mechanism of the reversible and unitary halting protocol 
and show how the state-locking pulse field works in the quantum control 
process to simulate efficiently the reversible and unitary halting protocol. 
From this simple physical model one can see clearly the reason why unitarity 
in quantum computation and hence the unitary quantum dynamics are so 
important. 

The unitarity of the quantum control process is of crucial importance 
if the quantum circuit Q c is used to build up an efficient quantum search 
process based on the unitary quantum dynamics. It follows from the quan- 
tum control process that each possible wave-packet state \l,CM0(j),R (j)) 
of Eq. (31) for j = 1,2, ...,m r corresponds one-to-one to a unique arriv- 
ing time Tj (see (27)) and a wave-packet spatial position R2,j(t mr ) at the 
end time t mr of the computational process (see Eq. (24)). On the other 
hand, it follows from the quantum circuit Q c that a different initial func- 
tional state \f r (xo)) with the index xo = 1,2, ...,m r corresponds one-to- 
one to a different wave-packet spatial position R2j{t mr ) and also a differ- 
ent arriving time Tj for the index j = (xf — Xq + l)modm r . Therefore, 
an initial functional state \f r (xo)) f° r = 1,2, ...,m r corresponds one-to- 
one to an output wave-packet state |1, CM0(j), Ro(j)) of Eq. (31) with 
the index j = (xf — xq + 1) modm r . The one-to-one correspondence is en- 
sured by the wave-packet spatial-position order such as the inequality (25): 
R2,m r {t mr ) < ■■■ < R2,2(tm r ) < R2,i(tm r ) an d also the arriving-time order such 
as the inequality (27): T\ < T 2 < ... < T mr . The one-to-one correspondence 
between the initial functional states {\f r (xo))} and the output wave-packet 
states {|1, CMO(j), R (j))} is not only necessary for the reversible and uni- 
tary halting protocol itself, but also it is of crucial importance if the quan- 
tum circuit Q c containing the quantum control unit that simulates efficiently 
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the reversible and unitary halting protocol is used to construct the quantum 
search process. The quantum search process requires that the probabilities of 
the desired state |1, CMO, R ) in the wave-packet states {|1, CMO(j), R Q (j))} 
(see Eqs. (31), (36), and (38)) become closer and closer to unity as the qubit 
number increases. This requires further that the time-compressing factor 
(v Q /v) or the atomic motional velocity v become smaller and smaller, as 
shown in Eqs. (36) and (38). Note that the atomic motional velocity v is 
limited by v << c (the speed of light in vacuum). Therefore, as the qubit 
number increases the difference among these m r possible arriving times {Tj} 
(see Eq. (28)) becomes smaller and smaller and so does the difference among 
these m r possible wave-packet spatial positions {R2,j(t mr )} (1 < j < m r ) (see 
Eq. (26)). If now the unitarity of the quantum control process is destroyed 
slightly due to some factors such as imperfect parameter settings, then the 
wave-packet spatial-position order such as the inequality (25) and/or the 
arriving-time order such as the inequality (27) could be easily destroyed in 
the quantum control process when the qubit number is large due to the small 
differences among these wave-packet spatial positions {R2,j{t mr )} and among 
these arriving times {Tj}. Then the one-to-one correspondence between the 
initial functional states {\f r (x ))} and the output states {|1, CMO(j), R (j))} 
could be destroyed easily and this will directly result in that the quantum 
search process based upon the quantum circuit Q c could not work normally. 

In particular, the importance of the unitarity of the quantum control 
process can be seen more clearly from the inverse process of the quan- 
tum control process. It is known in the quantum circuit Q c that these 
m r possible initial functional states {\f r (xo))} ar e different from each other 
and also orthogonal to one another and so are these m r possible initial 
states {\0,CM0, R )\0)\f r (x ))} of the quantum circuit. By the quantum 
circuit Q c these m r well-distinguished initial states are converted one-to- 
one into these m r output states {|l,CM0(j), Ro(j)) |1) |0)} for j = (xf — 
Xq + l)modm r = l,2,...,m r . Now for the 2 n — dimensional Hilbert state 
space S(Z 2 n) with a large qubit number n which is taken as the quantum 
search space it is required by the quantum circuit Q c that all these m r pos- 
sible probabilities {|A,(CM0, R , T t , T\)| 2 } of the state \1,CM0,R ) in the 
wave-packet states {|1, CM0(j), Ro(j))} of Eq. (31) be close to unity, that 
is, (^(CMO,^,^,^)! 2 > (1 - \np(n)/n) for j = 1,2, ...,m r , so that the 
quantum search process based on the quantum circuit can be made effi- 
cient. Thus, this means that all these m r wave-packet states {|l,CM0(j), 
R (j))} are very similar to the same state \1,CM0,Rq), indicating that all 
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these states {|1, CMO(j), R (j))} are very similar to each other and not 
yet orthogonal to one another. Thus, all these m r possible output states 
{|1, CMO(j), Ro(j)) |1) |0)} of the quantum circuit are very close to the same 
state |1, CMO, i?o)|l)|0). Hence they are very similar to each other and not 
yet orthogonal to one another. The difference among these output states be- 
comes also smaller and smaller as the qubit number increases. This property 
for these output states is completely different from that one for their corre- 
sponding initial states. When these output states {|1, CMO(j), Ro(j))\l)\0)} 
are transferred one-to-one back to the initial states {|0, CMO, Ro)\0)\f r (xo))} 
by the inverse processes of both the quantum control process and the quan- 
tum computational process of the quantum circuit the small difference among 
these output states is amplified to be a large one among those initial states, 
that is, almost indistinguishable output states are changed one-to-one back 
to well-distinguished initial states. Then some output states could not be 
converted correctly into their corresponding initial states by the inverse pro- 
cesses if the unitarity of the inverse processes is slightly destroyed. Therefore, 
the unitarity of the inverse processes is of crucial importance particularly in 
the first time stage of the whole inverse process that these output states 
are converted one-to-one back into those initial states. Here the first time 
stage of the whole inverse process is really the starting time stage (between 
T t + Ti and t mr ) of the inverse process of the quantum control process in the 
quantum circuit. 

Evidently, these m r possible initial states {\0,CM0, Ro)\0)\f r (xo))} of 
the quantum circuit have the same global phase factor and are orthogonal 
to one another. On the other hand, their corresponding m r possible out- 
put states {|1, CMO(j), R (j))\l)\0)} are very similar to each other and not 
yet orthogonal to one another, but the global phase factors for these output 
states may be very different from each other. This shows that the difference 
among these discrete and orthogonal initial states is mainly transformed to 
the global-phase difference among these similar output states by the quan- 
tum program or circuit Q c which consists of the time- and space-dependent 
unitary evolution processes. The conventional unitary operations which are 
generally space-independent are not generally suitable for constructing such 
a state unitary transformation. Only the time- and space-dependent unitary 
evolution processes can provide one with the best way to generate such a 
state unitary transformation. It is well known that there exists the square 
speedup limit for the standard quantum search algorithm whose unitary evo- 
lution process is space-independent and can be thought of as a single time- 
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dependent evolution process. Then a quantum search process could be able 
to break through the square speedup limit if it is constructed with the time- 
and space-dependent unitary evolution processes. The cost to break through 
the square speedup limit is that the quantum search space is enlarged and 
the single time-dependent unitary evolution process is extended to a time- 
and space-dependent unitary evolution process. It is also known [23] that 
an unknown quantum state can be transferred efficiently and completely to 
a large state subspace from a small subspace in the 2™— dimensional Hilbert 
state space of an n— qubit spin system whose unitary evolution process is 
space-independent, while the inverse unitary state-transfer process is rela- 
tively hard in the same Hilbert space. Indeed, when the Hilbert state space 
is added extra space dimensions and hence the unitary evolution processes 
are extended to be in a multi-dimensional Hilbert state space of time and co- 
ordinate spaces, the inverse state-transfer process that an unknown quantum 
state is transferred to a small subspace from a large one now could become 
more efficient than before, as shown in the paper. However, there still exists 
a difference between the unitary state-transfer process and its inverse pro- 
cess in the time and space multi-dimensional Hilbert state space, although 
the difference may be much smaller than that one in the 2™— dimensional 
Hilbert space in which any unitary evolution process is space-independent. 
The results in the paper show that an unknown quantum state of a large 
state subspace may be transferred to a given state of a small state subspace 
in a probability close to 100%, but it is impossible due to the limitation 
of unitarity that an unknown quantum state is transferred completely (in 
the probability 100%) into a small state subspace from a large one by any 
given unitary quantum dynamical process. On the other hand, it has been 
shown that an unknown quantum state always can be transferred completely 
(in the probability 100%) into a large state subspace from a small one [23]. 
One therefore obtains an important theorem: it is impossible to transfer com- 
pletely (in the probability 100%) every quantum state of a large state subspace 
into a small state subspace in the Hilbert space of a quantum system by any 
given unitary dynamical process quantum mechanically. This theorem could 
play an important role in understanding deeply the time evolution process of 
a quantum ensemble from a non-equilibrium state to the equilibrium state 
from viewpoint of the unitary quantum dynamics. It is well known in sta- 
tistical physics that such a non-equilibrium evolution process in a quantum 
ensemble is generally described through the stochastic probability theory. A 
direct result of the theorem is that there may exist a computational-power 
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difference [23, 24] between a time- and space-dependent unitary dynamical 
process, which transfers an unknown quantum state from one state sub- 
space to another, and its inverse process in a quantum system consisting of 
many quantum bits. The computational-power difference for the quantum 
search problem is very large for the discrete 2 n — dimensional Hilbert space 
of an n— qubit quantum system whose unitary evolution processes are space- 
independent, but it could become much smaller when the quantum system 
permits the time- and space-dependent unitary evolution processes in the 
quantum search process. Finally, it can be believed that a quantum system 
allowing to have time- and space-dependent unitary evolution processes will 
have a deep effect on the conventional quantum cloning theorem which usu- 
ally works in a quantum system with single time-dependent unitary evolution 
processes. 

The paper has not answered the question how high it is the mean motional 
energy of the unstable wave-packet state |1, CMl(t i), Ri(toi)) of the halting- 
qubit atom in the left-hand harmonic potential well so that the halting-qubit 
atom can overcome the intermediate potential barrier and enter almost com- 
pletely into the right-hand potential well by the quantum scattering process. 
One may determine the lower bound of the mean motional energy by solv- 
ing the quantum scattering problem of the atomic physical system. One 
needs also to consider the effect of the wave-packet spread on the quan- 
tum control process as the wave-packet spread of the motional state of the 
halting-qubit atom usually broadens during the atomic motional processes in 
the quantum control process. The effect may be investigated by solving the 
Schrodinger equation of the atomic physical system for the quantum control 
process. However, it can be expected that for a heavy halting-qubit atom the 
wave-packet spread should not have a significant effect on the lower-bound 
probability of Eq. (36) or (38). 
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